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should  there  be  an  accidental  surface  or  underwater  release  of  ammonia.  Although  the 
explosion  and  fire  hazard  of  such  a  spill  is  limited'  *  ''  the  health  and  pollution  haz¬ 
ard  is  not  insignificant.  The  United  States  Coast  Guard  (USCG)  recognized  the  need  to  l 
evaluate  large  accidental  spills  of  LNH-  and,  in  April  1972,  asked  Arthur  D.  Little,  Inp 
(ADL)  to  study  the  phenomena  associated  with  surface  and  underwater  ammonia  spills. 

The  purpose  of  the  study  was  to  develop  a  thorough  understanding  of  the  physio-chemical 
behavior  of  ammonia/water  reactions,  develop  analytical  model*,  and  conduct  a  scaled 
experimental  program  to  verify  and  modify  the  models  for  use  in  predicting  the  poten¬ 
tial  hazards  due  to  spills  of  one  up  to  3,000  tons  of  LNHg. 

This  report  contains  the  results  of  an  experimental  and  analytical  program  aimed 
at  evaluating  the  hazards  presented  by  a  3,000  ton  release  of  liquid  anhydrous  ammonia 
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1.  SUMMARY 


1.1  PURPOSE 

During  the  past  decade,  maritime  transport  of  liquid  anhydrous  ammonia 

(LNH^)  on  the  Inland  and  Coastal  waters  of  the  United  States  ^  has 

increased  rapidly.  Anhydrous  ammonia  is  carried  on  barges  as  a  high-pressure, 

liquefied  cargo  in  heavy-wall  steel  tanks,  and  as  a  low-temperature  cargo 

in  essentially  atmospheric-pressure  refrigerated  tanks.  The  .'frequent 

movement  of  bulk  quantities  of  LNH^  on  United  States  waters  poses  a  potential 

hazard  to  the  public,  to  aquatic  life,  and  to  other  shipping  should  there  . 

be  an  accidental  surface  or  underwater  release  of  ammonia.  Although  the 

(2  3) 

explosion  and  fire  hazard  of  such  a  spill  is  limited,  *  the  health 
and  pollution  hazard  is  not  insignificant. 


The  United  States  Coast  Guard  (USCG)  recognized  the  need  to  evaluate 
large  accidental  spills  of  LNH^  and,  in  April  1972,  asked  Arthur  D.  Little, 
Inc.  (ADL)  to  study  the  phenomena  associated  with  surface  and  underwater 
ammonia  spills.  The  purpose  of  the  study  was  to  develop  a  thorough  under¬ 
standing  of  the  physio-chemical  behavior  of  ammonia/water  reactions, 
develop  analytical  models,  and  conduct  a  scaled  experimental  program  to 
verify  and  modify  the  models  for  use  in  predicting  the  potential  hazards 
due  to  spills  of  up  to  3,000  tons  of  LNH^. 

1.2  SCOPE 


To  achieve  the  objectives  of  this  investigation,  ADL  established 
and  conducted  a  program  consisting  of  four  major  tasks. 


Task  I  — 
Task  II  — 
Task  III  - 
Task  IV  — 


Study  LNH^  Spills  on  Water 
Model  LHHj  Spills  on  Water 
Study  LNH~  Release  Underwater 
Recommend  Neutralization  Methods 


Both  Task  I  and  Task  III  were  conducted  in  three  stages.  First, 
a  series  of  laboratory  tests  was  conducted  with  spills  of  up  to  2.5  liters 
of  liquid  ammonia.  Second,  tests  were  conducted  in  the  open  in  a  small 
swimming  pool,  with  one  gallon  and  five  gallon  spills.  Finally,  quantities 
of  up  to  50  gallons  were  spilled  on  the  water  in  a  pond.  For  underwater 
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release,  there  were  0.4  gUloe  laboratory  teete,  on.  gallon  Intermediate 
scale  tests,  and  flee  g.llct  large-ecale  taots.  Discharge  depths  varied 
from  4.5  inches  to  60  inches. 

1.3  FINDINGS 

Hie  objective  of  this  work  was  to  predict  the  vapor  and  water  disper¬ 
sion  hazard,  of  up  to  3,000  ton  um3  spin,  on  water,  these  predictions 
are  given  in  Appendix  C. 

1.3.1  Laboratory  Experiments 


1.3. 1.1  Surface  Spills 

The  mean  partition  ratio,  the  fraction  of  the  total  weight  of  spilled  . 
ammonia  that  goas  into  solution  with  water,  was  about  0.73.  In  over  30%  of  the 
experiments,  the  partition  ratio  was  between  0.65  and  0.82;  was  insensitive  to 
quantity,  rate,  ana  orientation  of  spill;  and  was  independent  of  the  water 
salinity  (up  to  3.5%  NaCl)  and  air  or  water  motion.  In  addition,  restricting 
the  reaction  area  on  the  water  surface  did  not  cause  any  grei  variation 
in  the  partition  ratio.  A  slight  change  in  the  partition  ratio  was  observed 
with  changes  in  water  temperature. 

Surface  spills  reacted  with  water  in  a  small  area,  and  the  ammonium  hydroxide 
formed  spread  out  along  the  water  surface  at  about  2.5  Inches/second.  The  NH^OH 
temperature  was  about  5°  to  15°F  above  the  ambient  water  temperature. 

The  partition  ratio  results  and  the  water  temperature  data  can  be 
adequately  explained  by  a  thermodynamic  mixing  model. 

1.3 .1.2  Underwater  Release 

When  the  release  of  LNH^  was  at  shallow  depths  and  at  low  velocity, 
there  was  no  substantial  difference  in  the  partition  ratio  compared  to  the 
surface  spill  mean  value.  In  general,  when  the  release  was  at  depths  greater 
than  about  10  pipe  diameters,  the  partition  ratio  was  between  0.85  and  0.95. 

1*3.2  Intermediate-Scale  Experiments 

1.3. 2.1  Surface  Spills 

The  average  value  of  the  partition  ratio  was  about  0.56  fer  both 
one  gallon  and  five  gallon  spills,  and  increased  to  about  0.66  for  continuous 
spills.  The  water  temperature  rise  at  the  surface  was  about  5*  to  10°F,  al¬ 
though  in  some  tests  a  change  of  as  much  as  22°F  was  observed.  The  elevated 
temperature  lasted  for  about  five  to  ten  minutes. 
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Liberated  vapor  dis«*rs«d  as  c  cloud  end  passed  across  the  instrument 
rales »  15  feet  downwind.  This  was  indicated  by  a  momentary  dip  in  the  air 
temperature  recorded  at  the  rake?  which  implies  that  the  core  of  the  cloud 
was  still  cold  and  unaixed.  Distribution  of  the  vapor  mass  collected  by 
the  laplngers  on  the  instrument  nks  could  be  reasonably  fit  with  Gaussian 
profiles.  The  partition  ratios  measured  by  analyzing  the  pool  water  samples 
were  close  to  those  calculated  from  vapor  samples  collected  in  the  impingers. 

1.3. 2. 2  Underwater  Release 

At  an  average  value  of  0.85  to  0.90,  the  partition  ratio  values  were  consi¬ 
derably  higher  than  those  for  surface  spills.  No  vapor  liberation  could  be 
observed.  Water  surface  agitation  was  minimal  for  release  depths  greater 
than  about  15  inches  and  release  rates  less  than  15  ml/ sec  (about  0.25 
gal/ain) .  A  maximum  water  temperature  rise  of  about  10*F  was  recorded 
by  a  set  of  thermocouples  five  inches  off  the  axis  of  the  discharge  pipe. 

1.3.3  Large-Scale  Experiments 


1.3. 3.1  Surface  Spills 

The  partition  ratio,  calculated  from  tests  that  had  good  Implnger 
data,  was  between  0.53  and  C.63.  The  partition  ratio  date  shows  a  wide 
scatter  because  the  vagaries  of  the  wind  often  caused  the  vapor  cloud 
to  miss  the  instrument  rafts.  In  low  wind,  the  vapor  rose  rapidly. 

The  water  temperature  in  the  spill  zone  (recorded  by  thermocouples  at  0.5 
inch  to  1  inch  depth)  showed  substantial  heating  of  the  surface  layer, 
with  up  to  50°F  increase  noted.  The  temperature  rise  was  very  gradual, 
and  reached  peak  temperature  in  six  to  eight  minutes.  The  average  diameter 
of  the  surface  boiling  zone  was  about  20  to  25  feet. 

Vapor  cloud  behavior  in  the  wind  is  adequately  explained  by  a  plume 
theory  of  vapor  dispersion,  and  it  has  been  shown  that  the  water  tempera¬ 
tures  observed  could  result  form  an  adiabatic  mixing  of  water  and  liquid 
assEonia. 

1.3. 3. 2  Underwater  Release 

The  results  in  this  series  of  underwater  tests  are  not  much  different 
from  those  obtained  in  the  intermediate  size  tests.  A  10®  to  15®F  rise  in 
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water  temperature  was  noticed.  The  temperature  rise  was  in  phase  with 
the  duration  of  the  spill,  with  the  peak  temperature  being  reached  just 
at  the  end  of  the  spil3.  No  substantial  vapor  release  was  observed,  but 
violent  agitation  of  the  water  was  observed  directly  above  the  release 
pips  whan  the  release  depth  was  36  inches- 

1.4  CONCLUSIONS 

The  following  major  conclusions  can  be  drawn  from  the  results  of  the 
experimental  test  program  and  the  data  analysis: 

1.  The  partition  ratio  is  highest  in  the  case  of  underwater 
release  and  is  as  high  as  0.95.  Very  little  vapor  liberation 
occurs  in  the  case  of  underwater  release  at  depths  from 

10  inches  (for  small  quantities)  to  five  feet  (for  five  gallon 
releases) . 

2.  For  large  surface  spills ,  the  partition  ratio  depends  on  the 
spill  dynamics  and  varies  between  0.5  and  Q.6  (closer  to 
the  latter)  for  instantaneous  release  cn  the  surface.  For 
8 low,  continuous  release  on  the  surface,  partition  ratio 
values  tend  to  be  as  high  as  0.66.  However,  for  small 
spills,  under  controlled  laboratory  conditions  the  partition 
ratio  values  are  between  0.7  and  0  J 5,  with  a  mean  of  0.735. 

3.  In  the  case  of  large  surface  spills,  the  reaction  is  extremely 
rapid  and  liberates  a  dense,  white  fog  of  vapor,  probably 
containing  a  large  fraction  of  aerosols,  in  low  wind  conditions. 

4.  The  vapor  puff  formed  is  very  buoyant  and  rises  into  the  air 
as  it  travels  downwind.  The  rate  of  rise  depends  on  the  wind 
velocity.  Under  low  wind  conditions  the  cloud  forms  a 
characteristic  mushroom  cloud  before  dispersing.  The  path 

of  the  cloud  can  be  estimated  with  reasonable  accuracy  by 
existing  plume  theories.  Because  of  the  rapid  rise  in 
low  wind,  the  toxic  hazard  at  ground  level  is  smaller  for  low 
wind  than  for  high  wind. 

5.  The  boiling  zone  is  reasonably  small  —  about  25  fest  in  dia¬ 
meter  for  a  50  gallon  spill.  The  diameter  of  the  boiling  zone 
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increased  as  the  0.8th  power  of  the  quantity  of  spill  for 
the  tests  conducted.  However,  based  on  the  analysis  of  spread 
of  other  cryogens  on  water,  we  expect  that  a  0.375  power 
relation  with  quantity  is  more  appropriate  for  really  massive 
"instantaneous"  spills. 

6.  The  ammonium  hydroxide  formed  at  the  boiling  zone  stays  close 
to  the  top  of  the  water  surface  and  spreads  radially  at  about 
0.2  ft/sec.  The  depth  of  this  layer  is  of  the  order  of  a 
few  inches  (at  best  one  foot) .  Because  of  this  layered  spreading 
of  NH^OH,  the  hazard  to  aquatic  life  just  beneath  a  spill  site 
is  small. 

In  short,  we  conclude  that  a  reasonable  estimate  of  the  partitioning 
for  a  massive  spill  on  the  water  surface  would  be  0.6  into  water  and  0.4 
into  vapor.  The  downwind  vapor  hazard  at  ground  level  is  not  as  severe 
as  earlier  believed  because  of  the  buoyant  nature  of  aauonia  vapor. 

The  main  achievements  of  the  program  have  been  to: 

1.  Show  that  the  LNG^/water  surface  reaction  is  rapid  and 
confined  to  a  reasonably  small  area; 

2.  Show  that  the  ammonia  vapor  liberated  is  buoyant  and  rises 
rapidly; 

3.  Develop  a  thermodynamic  model  that  accurately  predicts  the 
upperbound  of  the  partition  ratio  for  surface  spills.  We 
feel  that,  though  the  reaction  is  thermally  limited, 

the  dynamics  of  the  spill  do  affect  the  partitioning; 

4.  Develop  a  method  for  indirectly  calculating  the  partition 
ratio  from  vapor  sample  measurements. 

5.  Develop  a  model  that  indicates  that  under  certain  accident  situations 
the  present  design  of  refrigerated  ammonia  tanks  for  barging  is 
inherently  safe  even  if  the  barge  sinks;  and 

6.  Show  that  the  continuous  underwater  release  of  LNHg  at  depths 
greater  than  ten  pipe  diameters  results  in  all  the  liquid 
dissolving  in  water  with  very  little  vapor  production. 
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1.5  UNRESOLVED  QUESTIONS 


Though  ch  has  been  achieved  in  the  present  program,  questions 
have  remain  unanswered,  either  because  of  the  slice  of  the  teats  conducted 
or  due  to  the  difficulties  in  measurement.  The  important  unknowns 
ares 

1.  The  amount  of  aerosol  in  the  vapor  liberated  in  a  surface 
apill  and  its  relationship  to  the  mode  of  spill  (instantaneous, 
continuous,  etc.)  and  the  rise  of  the  vapor  cloud; 

2.  The  rapid  rise  of  the  vapor  cloud  for  a  water  spill  as  against 
the  ground-hugging  clouds ,  possibly  due  to  greater  aerosol 
formation  in  land  spills  reported  in  the  literature; 

3.  The  possibility  of  underwater  explosions  in  the  case  of 
instantaneous  underwater  release  of  lsrge  quantities  of  LNHy 
(The  present  program  experiments  were  limited  to  e  continuous 
release  of  a  maximum  of  five  gallons,  and  it  is  difficult 

to  extrapolate  to  the  consequences  of  a  large  instantaneous 
release);  and 

4.  Validation  of  scaling  laws  for  the  determination  of  the 
maximum  radius  of  the  boiling  zone  (pocl  radius). 

1.6  RECOMMENDATIONS  FOR  ADDITIONAL  TESTS 

As  a  result  of  this  test  program,  we  are  now  in  a  position  to 
extrapolate  the  findings  of  this  work  to  predict  the  hazard  presented 
by  a  3,000  ton  release  of  ammonia  on  water.  Appendix  C  of  this 
report  contains  the  predictions  of  hazards  for  o  3,000  ton  (about 
one  million  gallone)  release  of  LNH^  on  water«  However,  because  of 
some  of  the  unresolved  questions  discussed  above  and  the  initial 
buoyant  behavior  of  the  ammonia  vapor  cloud,  we  feel  that  further 
experimentation  on  a  much  larger  scale  would  serve  to  increase 
confidence  in  the  analytical  models  and  provide  optimum  benefits  from 
the  overall  teat  prograa.  We  recommend,  therefore,  that  a  final 
series  of  testa  he  conducted  utilizing  3,000  gallon  releases  in  a 
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300-500  ft.  diameter  pond.  A  test  program  with  spill  quantity  at 
least,  two  orders  of  magnitude  greater  than  the  largest  tests  in  the 
present  series  would  provide  an  excellent  data  base  with  which  to 
test  the  scaling  laws  and  improve  confidence  in  theories  of  vapor 
movement . 
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2.  INTRODUCTION 

2.1  BACKGROUND 

During  the  past  decade,  maritime  transport  of  liquid  anhydrous 
ammoaia  (LNB^)  on  the  Inland  and  Coastal  waters  of  the  United  States 
has  increased  rapidly.  Anhydrous  anmonia  is  carried  on  barges  as  a 
high-pressure  cargo  in  heavy-wall  steel  tanks,  and  as  a  low-teaperafcure 
cargo  in  atmospheric-pressure  refrigerated  tanks.  There  has  been  very 
rapid  growth  of  refrigerated  transport,  motivated  largely  by  the 
substantial  savings  in  capital  costs  associated  with  onshore  refrigerated 
storage  compared  to  pressure  storage. 

The  frequent  movement  of  bulk  quantities  of  LKH»  on  United  States 
waters  poses  a  potential  hazard  to  the  public,  to  aquatic  life,  and  to 
other  shipping,  should  there  be  an  accidental  release  of  ammonia  either 
on  or  underwater.  A3. though  the  explosion  and  fire  hazard  of  such  a 
spill  is  limited,  the  health  and  pollution  hazard  is  not  Insignificant. 

There  is  a  definite  need,  therefore,  to  study  the  hazards  associated 
with  large  accidental  spills  of  LNH^  and  to  develop  emergency  procedures 
to  neutralize  the  effects  of  such  spills. 

The  United  States  Coast  Guard  (USCG)  recognized  the  hazards  presented 
by  a  spill  of  LNH^  on  the  inland  waterways  and  felt  the  need  to  evaluate 
the  potentialities  of  such  a  hazard.  USCG  therefore  initiated  a 
research  program  to  study  the  important  physical  aspects  of  LNH^ 
spill  on  end  underwater.  The  main  objectives  of  the  study  were 
to  conduct  scale  experiments  and  develop  theoretical  models  so  that  the 
results  of  small-scale  experiments  could  ba  extrapolated  to  predict  the 
hazard  caused  by  spills  of  up  to  3,000  tons  of  LNH^* 

At  atmospheric  pressure,  liquid  ammonia  (LNH^)  boils  at  about  -28°F; 
at  a  vapor  pressure  of  about  108  psla,  it  boils  at  60°F.  The  ammonia  can 
be  as  a  liquid  either  refrigerated  or  at  ambient  temperature  (under  the 
associated  vapor  pressure) ^ .  Ammonia  is  very  soluble  in  water,  and 
dissolution  process  is  exothermic. 
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Releases  of  ammonia  are  harmful  for  two  primary  reasons:  first, 
the  vapor  is  toxic;  and  second,  a  solution  of  ammonia  in  water  (ammonium 
hydroxide)  is  deadly  to  both  flora  and  fauna. 

In  general,  spills  of  ammonia  will  result  in  both  dissolution  and 
vapor  evolution,  and  the  relative  percentage  of  both  forms  depends  largely 
on  the  manner  in  which  LNHj  and  water  combine.  For  example,  a  water 
surface  spill  leads  to  rapid  boiling  of  some  ammonia,  and  simultaneous 
mixing  of  the  remainder  with  the  water.  Thus,  in  some  ways,  the  process 
is  similar  to  spills  of  cryogens  on  water  —  with  the  added  complication 
of  an  exothermic  reaction  occurring  at  the  LNH^/water  interface. 

2.2  HAZARDS  OF  AMMONIA 

2.2.1  Fire  Hazard 

Ammonia  is  not  highly  combustible.  Its  flamability  limits  ir.  dry 
air  are  15  to  28%  (by  volume).  It  is  not  known  how  the  presence  of  water 
vapor  affects  these  limits,  but  most  likely  it  will  tend  to  quench  a 
flame  and  reduce  the  fire  hazard.  Furthermore,  the  ignition  temperature 
of  aamonia  is  relatively  high  (1500°P  in  a  quartz  bomb  and  1204 °F  in  an 
iron  bomb). 

If  liquid  ammonia  remains  on  the  water  surface  for  any  reasonable 

length  of  time  before  it  evaporates,  and  is  accidentally  ignited,  it 

(2) 

may  burn.  Tests  by  Husa  and  Buckely  in  which  LNH^  was  poured  into 
a  3-ft  by  3-ft  by  2-inch  deep  pan  and  ignited  showed  that  brief  flashes 
occurred  when  the  ignition  source  was  brought  near  the  liquid  surface, 
but  that  the  flames  could  not  be  sustained. 

There  are  no  similar  experiments  with  ammonia  on  water.  Even 

(3) 

if  we  assume  that  LNH^  bums  on  water,  existing  correlations  show 
that  ammonia  flame  heights  should  be  lower  than  those  calculated  for 
organic  liquid  fuel  fires  of  the  same  diameter.  Furthermore,  ammonia 
fires  should  not  be  luminous  because  no  carbon  is  present.  Thus, 
with  less  height  and  nonluminousity,  acmonia  fires  should  not  radiate 
as  much  as  organic  fuel  fires. 
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2.2.2  &cplos£ca  Hazard 


The  accumulation  of  flamnable  vapor  ia  an  enclosure  (such  as 
a  house  on  the  shore  or  In  a  vessel)  and  its  subsequent  Ignition  may 
le^d  to  an  explosion.  Ammonia  explosions  are  known  to  have  occurred 
in  industry.  The  force  of  such  an  explosion  depends  on  the  degree  of 

(4) 

confinement.  Buckley  and  Husa  found  that  the  magnitude  of  the 
pressure  ratios  observed  in  ammonia  explosions  is  roughly  three-quarters 
that  ...jund  for  common  hydrocarbon  fuels  under  similar  conditions.  Attempts 
by  the  same  authors  to  detonate  amoonia/air  mixtures  at  atmospheric 
pressure  have  failed. 


2*2.3  Toxicity  Hazard 

The  major  hazard  from  an  ammonia  vapor  cloud  is  its  toxicity,  because 
ammonia  is  an  extremely  irritating  gas  and  is  toxic  in  relatively  small 
concentrations.  Exposure  to  0.25-0. 65 J  NH,  in  air  for  30  minutes  is 
sufficient  to  cause  death  or  serious  injury,  and  most  deaths  from  ammonia 
spill  accidents  on  the  ground  have  been  attributed  to  toxicity.  The 
American  Conference  of  Governmental  Industrial  Hy genie ts  recommends 
a  threshold  limit  value  (TLV)  of  50  ppm  in  air  (18  mgra/m  of  air)  for 
eight-hour  exposures.  Fortunately,  the  human  olfactory  system  is  capable 
of  detecting  the  presence  of  ammonia  at  20  ppm.  long  before  the  toxic 
limit  is  reached 


In  the  case  of  excessive  exposure,  reaction  of  anhydrous  ammonia 
occurs  in  the  upper  respiratory  tract  and  spasms,  inflammation,  or 
edema  of  the  larynx  results.  LNH^  can  cause  severe  caustic  skin  burns. 
Current  practices  of  first  aid  include  removal  to  uncontaminated  regions, 
breathing  pure  oxygen  for  short  periods,  and  flushing  skin  bums  and  clothes 
with  water  ^ .  In  case  of  after  effects  (usually  laryngeal  spasms) 
medical  examination  Is  recommended. 


Hazards  Hue  to  Dispersion  in  Water 


The  short-term  effects  of  large  quantities  of  aeaonia  dissolving 
in  water  are  the  rise  in  temperature  and  pH  level  of  the  water.  Both 
of  these  can  cause  severe  harm  to  aquatic  life  in  the  vicinity  of  the 
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the  spill.  The  water  temperature  .increase  could  also  cause  pressure 
buildup  in  sunken,  slowly  leaking  low-pressure  refrigerated  tanks. 
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2.2.5  Explosion  Hazard  Due  to  Large  Underwater  Release 

Saall-scale  experiments on  tjje  release  of  LNH^  a*»d  other 
cryogens  under  water  have  determined  that  there  is  a  possibility  of 
pressure  explosions  occurring  underwater.  The  physical  aspects 
of  such  explosions  are  not  fully  understood  at  this  time,  but  a 
possible  explanation  lies  in  the  hypothesis  of  extensive  superheating 
of  the  cryogen. 

2.3  PREVIOUS  INVESTIGATIONS 

Few  academic  studies  have  been  carried  out  in  which  one  liquid 
is  boiled  by  contacting  it  with  another  hot  liquid.  Those  that  have  been 
reported  involve  steady  state  experiments  with  organic  liquids  over 
hot  mercury.  They  are  summarized  briefly  in  a  paper  by  Nakanishi  and 
Reldv  1 .  In  this  same  paper,  the  available  experience  involving 
liquid  natural  gas  (LNG)  spills  is  also  discussed  and  pertinent 
references  given.  The  LNG  spill  studies  were  carried  out  primarily 
by  industrial  concerns and  by  the  Bureau  of  Mines The 
major  objective  was  to  ascertain  the  rapidity  of  boil-o£f  should 
LMG  (b.p.  *  -258°F)  be  spilled  in  an  accident.  Such  data  could  be 
used  aa  input  to  a  vapor  dispersion  program  to  predict  downwind 
concentrations.  However,  another  significant  objective  was  to  clarify 
an  unusual  phenomenon  first  observed  by  the  Bureau  of  Mines;  i.e., 
in  a  few  cases  where  LNG  was  spilled  on  water,  a  pressure  explosion 
resulted.  The  subsequent  studies  by  Shell and  Esso as  well 
as  by  Nakanishi  and  Reid^\  indicated  that  such  explosions  might  be 
expected  if  the  cryogen  could  superheat  and  initiate  vaporization  by 
homogeneous  oucleation. 

(13) 

A  major  investigation  using  ammonia  is  reported  by  Resplandy  . 

The  tests,  conducted  in  France,  studied  releases  on  land  for  three 
situations;  namely,  liquid  released  from  a  pressurized  tank,  the 
dispersion  of  vapor  released  from  a  pressurized  tank,  and  the  spill 
of  liquid  ammonia  into  earthen  dikes.  The  results  from  these  tests 


indicated  that  for  pressurized  release  (up  to  300  kg/*in  froa  a 
6-atawspbere  pressure  tank)  a  large  fraction  of  the  vapor  was  in  the 
fere  of  aerosols  of  1C  to  30-sieron  droplets.  It  was  also  reported  that 
up  to  402  of  the  liquid  released  flashed,  resulting  in  the  formation 
of  stable  aerosols.  It  is  not  very  clear  froa  the  paper  what 
fraction  of  the  vapor  foraed  was  in  the  aerosol  phase. 

At  low  wind  conditions,  the  cloud  rose  initially  to  a  height 
of  about  20  aeters  in  a  0.5  a/sec  wind,  then  fell  to  the  ground  and 
dispersed  at  ground  level.  The  visible  vapor  cloud  traveled  just 
about  100  aeters  before  vanishing.  Further,  the  cloud  boundary  and  the 
boundary  of  the  odorous  region  were  very  dose.  In  th*  case  of  liquid 
sasonia  spill  into  enrthern  dykes,  aerosol  fo nation  is  said  to  have 
been  noticed  in  the  first  few  seconds  after  spill.  However,  within 
a  few  alnutes,  the  vapor  generation  ceased,  leaving  a  stagnant  pool 
of  LHHy  The  paper  also  describes  the  teste  to  contain  the  vapor  by 
water  spray  resulting  in  no  major  improvements  in  the  cloud  travel.  Also 
described  was  a  test  to  ignite  the  cloud,  which  failed. 

Small  scale  tests ^  involving  the  pouring  of  cryogenic  liquids 
such  as  nitrogen,  ethane,  LEG,  etc.  on  a  water  surface  Indicate  the 
occurrence  of  agitated  boiling  at  the  water  surface.  However,  when 
asall  quantities  of  methane  were  poured  on  ethylene  glycol  (antifreeze), 
there  was  a  hiss /crack  noise  and  evaporation  was  quickly  over.  Laboratory 
testa  conducted  by  Reid  et  al. ^  show  that  the  LRH^  water  interaction 
restcbled  sore  the  methane/glycol  case  than  the  sethane/wster  case. 

The  vapor  produced  did  not  contain  any  water  fog. 

These  results  suggest  that  the  vapor  cloud  f rocs  a  spill  of 

(13) 

LNH,  on  water  say  differ  froa  the  results  reported  by  Resplandy  , 

"(14) 

Ball  ,  et  al.»  for  the  experiments  in  which  sasonia  land  spills 
were  studied.  For  a  spill  of  the  same  quantity  of  aaaonia,  the  vaporiza¬ 
tion  rate  on  water  is  ouch  greater  than  for  a  land  spill.  This  could 
give  a  shorter,  high-concentration ’’puff"  which  .sight  travel  further 
downwind  before  dispersing  safely  even  though  the  total  fraction  of 

szsxmia  vaporized  sight  be  less  because  of  solution  of  ssaonia  in  th- 

(15) 

water  phase.  Further,  limited  observations  of  land  spills 


indicate  the  presence  of  a  thick  white  fog  for  large  spills,  which 

(7) 

was  missing  in  the  small-scale  experiment  .  The  fog  is  caused  by 
cooling  the  surrounding  air  below  its  dewpoint.  This  phenomenon 
depends  largely  on  the  relative  humidity  of  the  air,  and  once  the  water 
vapor  in  air  precipitates  out,  it  dissolves  the  gaseous  ammonia  in 
the  vapor  cloud.  This  may  result  in  a  fog  which  csd  ,e  "heavy" 
enough  to  stay  close  to  the  ground. 

A  few  small  scale  underwater  release  test  results  have  been 
reported^*^ .  The  results  show  that  most  of  the  ammonia  dissolved 
in  water,  under  suitable  conditions  (4  to  5  inches  release  depths 
for  4  ml  LNH^) .  However,  it  is  reported  that  in  one  test  with  15°C 
water,  a  violent  explosion  destroyed  the  glass  vessel  containing 
the  water.  Based  on  Apfel's^^  superheat  explosion  theory,  it  was 
found  that  LNHg  may  be  superheated  to  a  temperature  of  342°K, 
corresponding  to  a  vapor  pressure  of  460  psia.  This  pressure  is  high 
enough  to  cause  an  explosion. 

Very  little  theoretical  analysis  is  available  in  the  literature 
concerning  the  ammonia-water  interaction.  For  the  dispersion  of  vapor 
ir.  air  and  NH^OH  in  water,  the  usual  air  dispersion  and  water  dispersion 
models  are  used. 

2,4  PRESENT  PROGRAM 

The  present  program  was  undertaken  with  a  view  to  answering  certain 
basic  questions  such  as:  how  much  spilled  LNH^  dissolves  in  water? 
what  is  the  relationship  of  the  dissolution  fraction  (called  the  partition 
ratio)  to  the  dynamics  and  the  quantity  of  spill,  if  any?  and  what 
are  the  features  of  the  vapor  dispersion  in  the  atmosphere  and  ammonium 
hydroxide  dispersion  in  water? 

To  achieve  the  objectives  of  the  program  indicated  in  Section L2, 
and  to  obtain  answers  to  the  above  questions,  a  three-stage  experimental 
program  was  undertaken.  The  primary  parameter  of  interest  in  all  the 
experiments  was  the  petition  ratio  (that  is,  the  fraction  of  the  LNH^ 
spilled  that  goes  into  solution  with  water) .  The  dispersion  of  vapor 
in  the  atmosphere  was  also  studied. 
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First,  a  laboratory  program  was  completed  in  which  up  to  2.5  liters 
(0.7  gallons)  of  LNG^  were  spilled  on  the  water  in  a  small  aquarium  tank. 

Both  surface  spills  and  underwater  releases  were  conducted.  The  details 
of  the  experiments  and  the  results  are  given  in  Chapter  3.  Second,  spills 
were  conducted  in  the  open  atmosphere.  A  small  swimming  pool  was  used 
to  contain  the  water.  The  vapor  dispersion  was  measured  by  sampling  the 
ammonia  vapor  in  air.  In  this  series,  up  to  18  liters  (5  gallons)  of  LNH^ 
were  spilled.  Chapter  4  contains  the  details  of  the  second  stage  of  testing. 
Chapter  5  describes  the  final  part  of  the  experimental  program.  The 
tests  were  conducted  in  a  pond  and  quantities  of  up  to  180  liters  (50  gallons) 
of  were  spilled  on  the  water  surface.  The  partition  ratio  was  es¬ 
timated  using  the  measured  vapor  concentration  in  the  air. 

Several  theoretical  analyses  have  been  made  to  describe  or  predict 
the  various  phenomena  that  occur  following  a  spill  of  LKHg.  Chapter 
6  deals  with  particular  aspects  of  the  spill  problem.  Section  6.2 
gives  a  vapor  flash  modal,  useful  for  calculating  the  amount  of  vapor 
produced  by  flashing  when  liquid  aaanonia  .leaks  from  a  pressurized 
container.  Section  6.3  describes  a  therodynamic  model  for  the  mixing 
of  LNH^  and  water.  It  can  be  seen  that  the  two  liquids  can  react  in 
three  different  ways,  leading  to  three  different  answers.  In  reality, 
all  three  reactions  may  take  place  simultaneously.  In  fact, 
some  of  the  experimental  results  tend  to  confirm  this.  Section  6.4 
presents  a  detailed  heat-transfer  model,  to  predict  the  pressure  rise 
inside  an  ammonia  tank  that  is  submerged  in  water.  Section  6.5  and 
6.6  treat  the  problem  of  dispersion  in  air  and  water  of  ammonia  vapor  and 
ammonium  hydroxide.  Section  6,7  discusses  two  theories  of  the  rise 
of  buoyant  vapor  in  the  atmosphere.  It  is  found  that  one  of  these 
theories  describes  quite  well  the  rise  of  ammonia  vapor  in  the  atmosphere. 
Section  6.8  gives  an  order-of-magnitude  type  analysis  for  underwater  release. 

Chapter  7  is  devoted  to  the  steps  to  be  taken  to  neutralize  an 
ammonia  spill  area.  Various  methods  are  discussed  from  containing  the 
area  with  booms,  to  air  sparging. 
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3.  LABORATORY  EXPERIMENTS 

3.1  PURPOSE  AND  SCOPE 

The  main  purpose  of  the  laboratory  testing  was  to  study  the  liquid 
ammonia-water  interaction,  under  controlled  conditions,  and  to  develop 
an  analytical  model  to  explain  such  interaction.  The  key  experimental 
result  was  the  partition  ratio,  i.e.,  the  fraction  of  the  total 
weight  of  spilled  ammonia  that  goes  into  solution  with  water* 

However,  in  addition  to  studying  the  dependence  of  the  partition 
ratio  on  the  test  variables,  we  also  tried  to  develop  an  under¬ 
standing  of  the  physical  mechanisms  involved. 

The  basic  experimental  steps  were  to: 

1.  Spill  a  predetermined  quantity  of  LNH^  on  water; 

2.  Allow  the  reaction  to  complete,  and  the  vapor  to  disappear,  then 
thoroughly  mix  the  water;  and 

3.  Titrate  the  water  samples  to  obtain  the  amount  of  ammonia 
dissolved  in  water  (and  hence  the  partition  ratio). 

3.2  APPARATUS  AND  PROCEDURES 
3.2.1  Test  Facility 

Figure  3-1  is  a  schematic  of  the  laboratory  test  apparatus. 

The  system  consisted  mainly  of  an  LNH^  transfer  system,  a  spill 
system,  and  a  test  tank.  Ammonia  was  transferred  as  a  two-phase 
fluid  from  a  storage  cylinder  to  a  separator.  The  vapor  formed  was 
vented  out  through  a  hood  and  the  saturated  liquid  was  collected  in 
a  spill  dewar.  The  ammonia  dumping  system  consisted  of  a  pivoted 
spill  dewar,  a  tipping  lanyard,  and  a  funnel  under  the  dewar  that 
had  interchangeable  discharge  tubes  for  controlling  delivery  of 
LNH^  to  water  tank.  The  testing  section  was  a  six-foot-long 
glass  tank,  18  inches  wide  and  20  inches  high,  half  filled  with  water. 

Both  the  separator  and  the  test  tank  were  enclosed  in  a  vented 
laboratory  hood  having  glass  doors  to  permit  observation  of  the  spill 
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FIGURE  3-1  LNH,  LABORATORY  TEST  FACILITY 


tests.  A  circulating  puipp  (not  shown  in  Figure  3-1)  was  provided  for 
rapid  mixing  and  draining  of  water  in  the  tank. 

3.2.2  Instrumentation  and  Measurement 

A  listing  of  the  techniques  used  and  the  parameters  measured 
during  the  tests  is  in  Table  3-1.  Phenolphthalein  was  added  to  the 
water  to  obtain  a  qualitative  indication  of  ammonia  concentration 
and  its  extent  in  the  test  volume.  In  several  of  the  tests,  concentra¬ 
tion  gradients  in  the  water  were  measured  with  pH  probes.  Following 
each  test,  the  partition  ratio  was  obtained  by  titrating  samples  of 
the  final  mixed  solution  from  the  tank. 

Temperature  profiles  for  both  the  liquid  and  vapor  phases  were 
recorded,  through  four  high-speed  thermocouple  channels  (about  20 
milliseconds  response)  and  three  slower-speed  channels,  on  a  multi¬ 
point  recorder.  (The  time  between  printing  of  each  point  was  about 
3.6  seconds.)  One  high-speed  channel  was  attached  to  the  spill  discharge 
nozzle  to  provide  a  spill  initiation  reference.  A  second  high-speed 
channel  was  attached  to  thermocouple  Rake  A,  and  the  remaining  two 
high-speed  channels  were  attached  to  Rake  C.  For  some  tests.  Rake  A 
and  Rake  C  were  combined  to  form  Rake  D,  a  three- thermocouple  rake 
with  three  high-speed  channels.  The  three-thermocouple  read-out 
on  the  multi-point  recorder  was  attached  to  a  common  Rake  B. 

Photography  has  provided  a  permanent  visual  record.  Many  tests 
were  recorded  with  real-time  movies,  while  for  others  the  high-speed 
Fastex  camera  was  used.  A  series  of  35  mm  stills  was  taken  of 
each  test,  particularly  during  the  latter  stages  of  the  test  program. 

3.2.3  Test  Procedures 

Preceding  each  test,  the  water  tank  was  flushed  several  times  to 
remove  any  residual  traces  of  ammonia  and  then  a  phenolphthalein  solution 
(typically  15  ml  of  a  solution  of  2.5  g  of  phenolphthalein  to  109  ml 
of  ethanol)  was  added  to  the  clean  water. 

After  the  thermocouples  were  checked,  the  exhaust  hood  was  secured. 
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Table  3-1 


Laboratory  Teat  Techniques  and  Parameters 

CONCENTRATION 

Gradients  -  Phenolphthalein  Indicator,  pH  Probe 
Final  Mix  -  Titration 

TEMPERATURE  -  LIQUID  (VAPOR) 

Nozzle  -  1  Hi-speed 
Rake  A  -  1  Hi-speed 
Rake  B  -  3  Multi-point  • 

Rake  C  -  2  Hi-speed 
Rake  D  -  Rake  C  +  Rake  A 

PHOTOGRAPHY 

Real-Time  Color  Movies 
Slow-Motion  Color  Movies 
Ultra-Slow-Motion  (Fastex)  Color  Movies 
35  mm  Color  Stills 


and  the  area  cleared  of  non-test  personnel,  ammonia  was  transferred 
from  the  storage  bottle  outside  the  test  building  into  the  spill 
dewar.  This  transfer  was  accomplished  by  opening  a  valve  in  the  neck 
of  the  LHHg  bottle  and,  controlling  the  flow  with  the  needle  valve 
in  the  transfer  line,  filling  the  graduated  spill  dewar  to  a  predetermined 
level.  With  the  dewar  filled,  a  countdown  was  begun  and  the  thermo¬ 
couple  recorders  and  cameras  were  started  at  predetermined  intervals 
before  the  actual  spill. 

Each  test  was  initiated  by  tipping  the  spill  dewar  with  a  lanyard, 
so  that  the  LNH^  dumped  into  the  spill  funnel  flowed  through  a  dis¬ 
charge  tube  onto  the  surface  of  the  water.  The  rates  and  angles  of 
LNH^  discharge  were  controlled  by  the  size  and  angle  of  the  discharge 
tube  below  the  spill  funnel.  Temperature  recording  and  photographic 
surveillance  were  usually  maintained  for  about  five  minutes  following 
completion  of  the  spill.  (Depending  on  spill  tube  size  and  spill 
quantity,  a  spill  was  normally  completed  within  15  to  45  seconds.) 

At  five  minutes  into  each  test,  mixing  of  the  tank  was  initiated. 

A  mixing  period  of  two  minutes  was  allowed  -  although  visual  observa¬ 
tion  of  the  dye  and  temperature  measurement  indicated  that  mixing 
was,  in  fact,  essentially  complete  after  about  40  seconds  -  then  analysis 
by  titration. 

The  total  quantity  of  ammonia  dissolved  in  the  water  was  determined 
from  the  measured  ccncentration  of  ammonia  in  the  final  mixed  solution 
and  the  known  volume  of  the  final  mixed  solution.  This,  in  turn,  was 
divided  by  the  initial  spill  quantity  to  obtain  the  partition  ratio. 

3.3  SURFACE  TESTS 

Eight  groups  of  tests  were  conducted  to  determine  the  dependency 
of  the  partition  ratio  on  various  physical  parameters.  Figure  3-2 
shows  the  relationship  of  the  five  test  categories  in  the  laboratory 
test  program.  Several  test  parameters  were  varied  for  each  category 
to  study  possible  influence  on  the  partition  ratio. 
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FIGURE  3-2  SCHEMATIC  DIAGRAM  OF  THE  LABORATORY  TESTS 


3.3.1  Group  I  -  Basic  Effects 

The  Group  1  testa  were  conducted  to  investigate  the  basic  effects 
of  spill  quantity,  rate,  and  method  of  delivery.  The  tests  are  listed 
in  Table  3-2.  In  this  group,  spill  quantity  ranged  from  800  to  2400  ml; 
spill  rates  were  varied  from  a  nominal  50  ml/sec  (from  a  3/8-inch 
diameter  discharge  tube)  to  a  nominal  100  ml/sec  (from  a  1/2-inch 
diameter  discharbe  tube) ;  and  the  spill  angle  was  either  parallel 
with,  or  vertical  to,  the  liquid  surface.  We  felt  that  Group  I 
would  be  the  key  to  establishing  good  base  data  and  understanding 
the  liquid-ammonia/water-surface  interaction.  Thus,  a  relatively 
large  number  of  tests  (22)  were  conducted  in  this  group. 

Table  3-2  lists  calculated  partition  ratios  that  are  plotted  in 
Figure  3-3  as  a  function  of  spill  quantity  and  plotted  in  Figure  3-4 
as  a  function  of  spill  rate.  The  results  show  that  the  partition 
ratio  is  essentially  Independent  of  spill  si2e,  spill  rate,  or  angle 
of  delivery.  With  two  exceptions,  all  experimentally  determined  parti¬ 
tion  ratios  fell  between  0.65  and  0.82,  or  within  s  0.11  of  the 
0.74  mean  value.  The  two  tests  that  did  not  fall  in  this  range  were 
rerun,  and  then  did  fall  within  the  expected  values. 

The  fact  that  considerable  difference  in  the  ammonia  delivery 
method  seemed  to  have  little  effect  on  partition  ratio  suggested 
that  the  partition  ratio  might  be  limited  more  by  thermodynamic  considera¬ 
tions  than  by  surface  reaction  rates. 

Test  observation  and  the  test  movies  indicated  that  ammonia 
evaporation  occurred  in  a  rather  small  boiling  zone,  the  diameter 
of  which  could  be  roughly  measured  by  observing  surface  waves  and 
the  location  of  the  vapor  cloud  leaving  the  surface.  Outside  of  this 
small  zone,  which  typically  seemed  to  be  about  eight  inches  in  diameter, 
little  vapor  generation  was  noted.  Surface  temperature  measurements 
and  the  phenolphthalein  dye  both  indicated  a  warm  layer  of  fluid, 
evidently  aqueous  ammonium  hydroxide,  propagating  along  the  surface 
away  from  the  boiling  zone.  Typical  values  for  the  thickness  of  the 
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FIGURE  3-4  PARTITION  RATIO  VERSUS  SPILL  RATE 


layer  ranged  from  two  to  four  inches.  The  dye  propagated  much  faster 
horizontally  than  vertically.  Figure  3-5  shows  two  photographs 
of  a  spill.  The  dye  dispersion  and  the  boiling  zone  can  both  be 
clearly  seen. 

Temperatures  were  measured  for  all  tests  in  the  liquid  near  the 
surface,  and  also  in  some  instances  in  the  vapor  space.  Figure  3-6 
shows  a  typical  liquid  temperature  vs.  time  history  at  1/8,  3/8, 
and  7/8  inch  below  the  liquid  surface.  The  p  otted  test  had  1550 
ml  of  ammonia  spilled  through  the  3/8-inch  diameter  flow  nozzle 
parallel  to  the  liquid  surface.  The  measurement  point  was  23.5  inches 
downstream  from  the  spill  point  (downstream  implies  the  injection 
direction  of  the  spill).  As  shown  by  the  plot,  a  rather  abrupt  rise 
in  surface  temperature  occurred  approximately  three  seconds  after 
3pill  initiation.  High-frequency  temperature  variations  occurred 
during  the  spill  (about  25  seconds),  particularly  near  the  surface.* 

This  phenomenon  appears  to  be  associated  with  wave  motion  induced  by 
the  spill.  After  completion  of  the  spill,  the  temperature  variations 
were  much  smoother,  which  suggested  a  much  smoother  NH^OH  flow.  Typically, 
a  maximum  5°  to  10* F  surface  temperature  rise  was  noted  following 
a  spill. 

It  was  possible  to  plot  the  propagation  of  the  leading  edge 
of  the  warm  ammonium  hydroxide  layer  leaving  the  spill  zone  from 
photographic  and  thermocouple  records.  Figure  3-7  shews  that  the 
front  had  a  mean  propagation  velocity  of  about  0.22  foot  per  second. 

This  plot  seems  to  roughly  correlate  with  other  observations 
indicating  the  radius  of  the  boiling  zone  to  be  about  eight  inches 
(at  zero  seconds), 

*The  temperature  trace  for  the  7 /8-inch  dnpth  was  manually  replotted 
from  another  chart.  In  this  process,  some  of  the  high-frequency 
components  may  have  been  lost;  however,  in  general,  temperatures 
measured  at  greater  depths  exhibited  less  change  and  less  high-frequency 
oscillation. 
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neuae  3-5  thotcgwk  showing  the  boiling  zons 

AND  THE  SPREADING  OF  THE  AMMONIUM 
HYDROXIDE.  [Exp#ri««jt  IS,  11/1V72) 
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FIGURE  3-8  TYPICAL  WATER  TEMPERATURE  RECORD 


o  Dye  Observations 

•  -  Temperature  Data  (At  Depths  indicated) 
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FIGURE  3-7  HORIZONTAL  PROPAGATION  Qr  NH4OH  SOLUTION 


3.3.2  Group  II  -  Temperature  and  Salinity  Effects 


The  Group  II  tests  were  conducted  to  investigate  the  influence 
of  water  temperature  and  salinity.  Table  3-3  lists  the  three  tests 
of  this  group  -  two  to  investigate  the  influence  of  water  temperature, 
and  one  to  investigate  the  influence  of  salinity.  Before  each  high- 
temperature  test,  the  water  was  nominally*  at  94°F;  prior  to  the  low- 
temperature  test,  the  water  was  37°F.  In  previous  tests,  water 
temperature  had  varied  between  45°  and  60°F.  Common  salt  (NaCl)  was 
added  to  the  test  water  to  achieve  a  salinity  comparable  to  sea  water 
(3.5%) . 

Figure  3-8  shows  the  partition  ratios  obtained  in  this  set  of 
tests.  Note  that  the  partition  ratio  was  not  significantly  affected 
by  either  water  temperature  or  salinity.  However,  high  temperature 
water  tended  to  give  a  slightly  lower  partition  ratio.  This  observa¬ 
tion  is  in  keeping  with  theoretical  calculations.  (See  Section  6.2.) 
Although  water  temperature  had  no  major  effect  on  the  partition  ratio, 
the  surface  temperature  rise  (in  the  propagating  NH^OH  solution) 
outside  the  boiling  zone  was  much  less  at  higher  water  temperatures. 
This  result,  shown  in  Figure  3-9,  also  correlates  with  predictions  of 
the  theoretical  analysis. 

3.3.3  Group  III  ~  Containment  Effects 

The  Group  III  tests  were  conducted  to  investigate  the  influence 
of  containing  the  spilled  ammonia  (and/or  propagating  ammonium 
hydroxide  layer)  on  both  the  liquid-ammonia /water-surface  reaction 
and  the  partition  ratio.  Containment  was  achieved  by:  1)  restricting 
the  spread  of  NH^OH  to  a  fraction  of  the  total  tapJt  surface  area  with 
a  plywood  baffle  across  the  water  tank;  and  2)  conducting  spill  tests 


*Ncminal  water  temperatures  for  each  test  were  measured  with  a 
mercury  in  glass  thermometer  prior  to  closing  the  hood  and  transferring 
the  ammonia  into  the  spill  dewar.  Typically,  about  five  minutes 
might  elapse  between  measuring  nominal  water  temperature  and  the  actual 
spill.  During  the  high-temperature  test,  some  surface  temperature  drop 
occurred  between  measurement  of  the  nominal  temperature  and  the  actual 
snlll  time  (so  the  surface  temperature  at  the  spill  initiation  was  closer 
to  90°F  than  to  94°F) . 
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FIGURE  3-8  PARTITION  RATIO  VERSUS  SPILl.  SIZE  FOR  DIFFERENT  WATER 
TEMPERATURES  AND  SALINITY 


in  a  six-inch  diameter  beaker.  Table  3-4  lists  the  data  for  the  six 
tests  conducted  in  this  group. 

Figure  3-10  shows  the  partition  ratio  data,  which  indicates 
that  putting  a  baffle  at  either  station  27  (37%  of  the  tank  left 
available  to  the  spill)  or  at  station  42  (58%  of  the  tank  left  available 
to  the  spill)  had  no  appreciable  effect  on  the  partition  ratio.  A 
beaker  test,  with  a  total  area  even  smaller  than  that  previously 
observed  as  the  boiling  zone  in  the  tank,  also  did  not  significantly 
influence  the  partition  ratio.  However,  the  reproducibility  of  the 
beaker  spills  was  not  as  good  as  it  had  been  in  the  larger  tank. 

Figure  3-11  shows  the  temperature  data  from  one  of  the  beaker 
tests.  Note  that  during  the  test  the  temperature  measured  in  the 
vicinity  of  the  surface  was  colder  than  the  initial  water  temperature; 
while  below  the  surface  the  temperature  was  somewhat  warmer  than  the 
initial  temperature;  and  at  about  the  midway  level  in  the  beaker  the 
temperature  is  uninfluenced  by  the  superimposed  layer  of  ammonia  and/or 
ammonium  hydroxide.  This  type  of  vertical  temperature  profile  probably 
occurs  because  of  stable  stratification  of  the  liquids,  with  the  layer 
of  coldest  LNH3  at  the  top,  followed  by  a  slightly  warmer  (but  never¬ 
theless  heavier)  NH^OH  solution,  overlaying  the  water.  The  vertical 
diffusion  of  the  NH^OK  solution  is  very  slow,  as  can  be  seen  by  the 
photographs  in  Figure  3-12. 

3.3.4  Group  IV  -  Interface  Reactions 

The  Group  IV  tests  were  conducted  to  study  the  influence  of 
water  and  air  motion  on  the  ammonia/water-surface  reactions  and  the 
resulting  partition  ratio.  Another  variable  studied  in  these  tests 
was  the  effect  of  isolating  the  liquid  surface  downstream  of  the  spill 
to  determine  if  any  appreciable  portion  of  the  dissolved  ammonia 
resulted  from  vapor/liquid  reactions  remote  from  the  actual  spill 
zone.  Table  3-5  lists  the  tests  and  data  of  this  group. 

Figure  3-13  shows  the  test  facility  modifications  required  for 
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FIGURE  3-10  EFFECT  OF  SPREAD  CONTAINMENT  ON  THE  PARTITION  RATIO 


Figure  3  12  Seconds  After  Spill 


Figure  3-1235  375  Seconds  After  Spill 


FIGURE  3-12  PHOTOGRAPHS  ILLl'STRATtNG  THE  SLOW  DIFFUSION  RATE 
OF  NH40H  SOLUTION  IN  THE  VERTICAL  DIRECTION 
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the  Group  IV  spill  tests.  A  transition  section  was  placed  in  the  vent 
hoed,  so  that  the  air  flow  induced  by  the  exhuast  fan  would  pass  directly 
over  the  water  tank  end  tend  to  purge  the  surface  of  a^onia  vapor. 

(In  prior  tests,  air  flew  had  been  around  the  tank,  hut  not  across, 
the  water  surface.)  A  horizontal  baffle  was  plaead  approximately 
one  inch  below  the  liquid  surface,  and  the  inlet  and  suction 
lines  from  the  circulation  pump  were  modified  to  induce  water  motion 
at  the  surface.  With  these  modifications,  it  was  possible  to  achieve 
an  air  velocity  of  about  15  feet  per  second  and  a  water  velocity  of 
about  0.5  foot  per  second.  The  water  surface  outside  tfca  spill  sons 
was  isolated  from  the  test  area  by  placing  plywood  sheets  on  the  water, 
except  for  about  eight  inches  on  either  side  of  the  spill  zone. 

Figure  3-14  shows  the  partition  ratios  obtained  in  these  tests. 

Mote  that  almost  all  of  the  experiments  had  a  partition  ratio  of 
0.67.  This  is  slightly  less  than  the  mean  of  0,74  obtained  in  the 
unmodified  experiments.  However,  neither  air  notion,  nor  water  motion, 
nor  surface  isolation  affected  the  partition  ratio  to  any  substantial 
degree  when  considered  separately. 

3.4  UNDERWATER  TESTS 

Fine  underwater  tests  were  conducted  to  determine  the  influence 
of  submergence  and  spill  mode  on  the  partition  ratio.  Table  3-6 
lists  the  data  for  these  tests. 

A  gravity  head  device  was  used  for  two  tests,  and  a  pressurised 
transfer  device  was  used  for  the  other  tests.  .Figure  3-15  shows 
schematics  of  both  pieces  of  apparatus.  The  gravity  head  device  was 
an  extension  of  the  spill  tube  used  in  the  previous  surface  tests. 
However,  it  was  limited  to  shallow  depths,  because  boiling  in  the  sub¬ 
merged  portion  of  the  tube  resulted  is  slow  discharge  response  (i.a,, 
rapid  subsurface  sswonia  discharge  did  not  occur  for  about  30  to  45 
seconds  after  liquid  aosoals  bad  been  dipped  into  the  spill  funnel). 
Further,  to  ensure  a  consistent  outflow  of  liquid  swoaia,  it  was 
necessary  to  bend  the  nozzle  upward  so  that  bubbles  forming  near  the 
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(b)  Spill  time  refined  after  viewing  test  movies. 

(c)  Very  preliminary  test  -  incomplete  data. 

(<*)  .03  orifice  in  discharge  line. 

<e)  .05  orifice  in  discharge  line. 


end  of  the  nozzle  would  enter  the  water  rather  than  back  up  Into  the 
discharge  tube.  The  pressurized  discharge  apparatus  was  not  depth 
limited,  nor  was  it  necessary  to  bend  the  discharge  tube  upward  to 
ensure  positive  outflow. 

Note  that  the  partition  ratios  obtained  with  gravity  discharge 
(at  a  depth  of  5.5  inches  below  the  water  surface)  are  not  very  different 
from  the  partition  ratios  obtained  with  the  surface  spills.  This 
agreement  may  be  due  to  the  shallow  release  depth,  because  the  ammonia 
vapor  may  not  have  had  time  enough  to  dissolve.  Thus,  the  results 
may  be  similar  to  those  with  surface  boiling.  However,  with  deep 
discharge,  most  of  the  ammonia  vapor  was  dissolved  by  the  water. 

Table  3-6  lists  the  increased  partition  ratios  obtained  for  pressurized 
releases  under  water.  Note  that  the  partition  ratio  is  high  for 
pressurized  release  even  at ' relatively  shallow  depths  (except  for 
one  test).  This  may  be  due  to  better  mixing  caused  by  the  high  velocity 
jet  that  issues  from  the  end  of  the  pipe  during  such  releases. 

3.5  FINDINGS 

3.5.1  Surface  Spills 

The  partition  ratio  seems  to  be  essentially  a  constant,  not 
influenced  by  changes  in  a  number  of  physical  variables.  Specifically, 
we  observed  that; 

9  In  over  90%  of  the  tests,  the  partition  ratio  was  between 
0.65  and  0.82.  In  almost  all  Instances,  where  data  was 
outside  this  range,  test  reruns  resulted  in  ratios  within  the 
expected  range. 

©  Observation  of  the  test  tank  indicates  that  evaporation  seemed 
to  be  confined  to  an  eight-inch  zone  surrounding  the  spill. 
Phenolphthalein  dye  and  temperature  measurements  showed  that 
a  warm  layer  of  ammonium  hydroxide  propagated  from  the  spill 
zone  along  the  water  surface.  The  layer  was  about  an  inch 
thick,  and  its  surface  temperature  was  about  5®  to  10°F  above 
the  initial  water  temperature. 


e  Within  the  range  of  the  experiment f  the  water  temperature  did 
not  seem  to  have  a  strong  effect  on  the  partition  ratio, 
although  for  higher  water  temperatures  a  smaller  temperature 
rise  was  measured  in  the  propagating  layer  of  ammonium  hydroxide  * 

•  Neither  the  partition  ratio  nor  the  physical  model  just  described 
seemed  to  be  strongly  affected  by  variables  which  might  be 
expected  to  influence  interface  heat  and  mass  transfer,  such 

as:  spill  quantity,  spill  rate,  spill  delivery  angle, 
air  and  water  motion,  or  salinity. 

0  Restricting  the  total  area  of  the  spill  site  did  influence 

propagation  of  the  ammonium  hydroxide  layer,  but  if  the  evapora¬ 
tion  zone  remained  unobstructed,  there  did  not  seem  to  be  any 
marked  affect  on  either  the  evaporation  process  or  the  parti¬ 
tion  ratio. 

•  When  the  area  available  to  the  spill  was  reduced  to  less  than 
the  nominal  eight-inch  boiling  zone  (such  as  in  the  beaker 
tests),  the  upper  layers  of  the  ammonium  hydroxide  solution 
underwent  a  temperature  reduction  rather  than  the  temperature 
rise  previously  noted.  However,  although  the  partition  ratio 
varied  more  in  these  tests,  the  average  value  was  not  greatly 
affected. 

3.5.2  Underwater  Release 

We  observed  that: 

•  When  the  release  was  at  low  velocity  and  close  to  the  water 
surface,  there  was  no  substantial  change  in  the  partition 
ratio  compared  with  the  mean  value  for  the  surface  spills. 

•  When  the  release  was  at  high  velocity,  and  at  depths  of  four 
and  nine  inches,  the  partition  ratio  values  increased  and 

were  much  closer  to  1.0.  This  may  be  a  consequence  of  increased 
turbulence  and  hence  better  mixing  caused  by  the  high  velocity 
jet  of  liquid  ammonia  exiting  from  the  tube.  Also,  the 


increased  residence  tine  of  the  vapor  in  the  water  (due  to 
great  depth)  results  in  redissolution  of  vapor  in  water, 
giving  a  higher  partition  ratio, 

3.6  DISCUSSION 

It  appears  that  the  partition  ratio  for  small-scale  tests  is  not 
very  dependent  on  the  physical  or  dynamic  variables  in  the  experiment 
and  remains  nearly  constant  -  within  a  reasonable  experimental  scatter  - 
at  an  average  value  of  0.735.  The  mode  or  total  quantity  of  surface 
release  do  not  seem  to  influence  the  partition  ratio.  This  holds 
true  even  when  the  spread  of  NK^OH  is  restricted  (as  during  the  beaker 
and  baffle  tests).  On  the  other  hand,  underwater  release  below  a 
certain  depth  seems  to  increase  the  dissolution  of  NH^  vapor  and 
thereby  increase  the  partition  ratio. 

In  the  Group  IV  tests  conducted  with  slightly  changed  apparatus, 
the  partition  ratio  was  consistently  at  0.67,  It  is  not  clear  whether 
this  was  due  to  airflow  over  the  water  surface,  or  water  motion,  or 
both.  However,  the  6%  decrease  in  partition  ratio  from  that  of  the 
surface  tests  might  very  well  represent  some  experimental  problems 
rather  than  any  physical  phenomenon. 

The  analytical  model  presented  in  Section  6.2  predicto  a  partition 
ratio  of  0.735  for  a  70®F  water.  Most  of  the  experiments  were  carried 
out  with  50°F  to  70®F  water,  and  the  agreement  between  experiment 
and  theory  is  remarkable.  This  indicates  that  the  reaction  between 
small  quantities  of  ammonia  and  water  is  purely  a  thermodynamic 
phenomenon  and  is  independent  of  other  kinetic  and  dynamic  variables. 

The  maximum  temperature  of  the  NH^OH  solution  measured  is  a 
function  of  the  initial  temperature  of  the  water.  The  temperature 
rise  is  about  15°F  for  water  at  35®F  and  almost  zero  for  water  at 
90°F.  Note,  however,  that  the  temperatures  were  measured  at 
about  three  to  four  reaction-zone  diameters  from  the  center  of  the 
spill.  Temperature  measurements  made  in  the  reaction  zone  did  not 
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give  very  meaningful  results  due  to  the  extremely  noisy  data  stream 
caused  by  violent  boiling  and  agitation.  The  thermodynamic  model 
predicts  that  the  temperature  rise  of  NH^QR  (14®F  for  32°F  water, 
and  no  rise  for  70°F  water)  is  a  function  of  the  water  temperature. 

The  measured  and  predicted  values  agree  fairly  closely. 

Visual  observation  of  th?.  NH^OH  layer  spread  indicates  that  it 

spreads  horizontally  with  a  mean  velocity  of  about  0.22  foot  per 

second,  as  shown  in  Figure  3-7,  The  vertical  propagation  of  NH^OH 

is  minimal.  Based  on  visual  observation  of  the  vertical  diffusion 

of  a  dye  marker  during  the  first  few  seconds  after  a  spill,  wa  calculate 

2 

that  diffusion  coefficient  for  NH.OH  is  between  0.1  and  0.2  cm  /sec 

2  4 
(0.4  to  0.8  tz/hr) . 


4.  INTERMEDIATE  SIZE  EXPmHEKTS 


4.1  PURPOSE  AND  SCOPE 

The  purpose  of  the  icteyaediate-scale  tests  was  to  measure  both 
the  partition  ratio  and  vapor  concentration  in  air,  simultaneously, 
on  as  large  a  scale  as  possible.  By  doing  so,  we  Intended  to  see  if 
the  partition  ratio  could  be  predicted  by  measuring  the  integrated 
vapor  concentration  in  air. 

The  small  spill  quantities  and  unnatural  environment  of  the 
laboratory  tests  (a  constricting  of  the  hood  and  vertical  air  move¬ 
ment  in  the  vent)  precluded  meaningful  measurement  of  vapor  concentra¬ 
tions  in  the  gas  phase.  Conversely,  during  any  large-scale  spill  it 
is  not  possible  to  directly  measure  the  partition  ratio,  because 
of  the  difficulty  in  obtaining  a  uniformly  mixad  water  sample.  There¬ 
fore,  indirect  methods  were  needed  to  determine  the  partition  ratio. 

These  experiments  were  based  on  the  premise  that  it  would  be  possible 
to  estimate  total  vapor  flow  by  obtaining  vapor  samples  over  a  specified 
time  interval  with  collecting  devices  at  a  few  locations  downwind  of  the 
spill.  Comparing  the  estimated  total  vapor  mass  liberated  during  a 
spill  with  the  quantity  of  LNH^  spilled  on  water  would  allow  calculation 
of  the  partition  ratio.  This  value  of  the  partition  ratio  can  then 
be  compared  with  the  value  obtained  by  titrating  samples  of  the  thoroughly 
mixed  water  on  which  the  LNH^  was  spilled. 

4.2  APPARATUS  AND  PROCEDURE 
4.2.1  Test  Facility 

Figure  4-1  shows  the  essential  features  of  the  test  facility:  a 
swiimlng  pool,  a  vapor-sampling  instrument  rake,  and  a  spill  bucket. 

The  aluminum-sided,  plastic-lined  swimming  pool  -  20  feet  in  diameter 
and  2  feet  deep  -  was  buried  in  the  ground  so  that  its  top  was  just 
above  ground  level  and  the  pool  water  surface  was  at  about  ground  level. 
The  LNH^  transfer  system*  tipping  mechanism,  discharga  nozzle,  water 
circulating  pirop,  and  drain  loop  were  vary  similar  to  tho.se  used  in  the 
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laboratory  test  facility.  The  poc.1  was  generally  filled  with  water 
using  a  fire  hose  connected  to  the  city  water  supply.  However, 
the  water  from  an  adjacent  pond  was  used  for  some  of  the  early 
tests.  The  test  facility  was  located  on  the  testing  grounds  of 
Arthur  D.  Little,  Inc. 

4.2.2  Instrumentation 

Vapor  concentration  measurements  were  made  with  impinger  (sniffer) 
devices  mounted  on  a  cross-shapad  rake,  as  shown  in  Figure  4-2, 
generally  placed  about  15  feet  downwind  of  the  pool  center.  The 
rake  also  contained  thermocouples  to  measure  ammonia  vapor  temperature . 
Pool  water  temperature s  were  also  recorded. 

The  impinger,  shown  in  Figure  4-3,  draws  a  metered  flow  of  sir 
through  a  boric  acid  solution,  which  dissolves  aspirated  ammonia  vapor. 
The  mass  of  dissolved  NH^  car.  be  determined  later  by  laboratory  titra¬ 
tion  techniques.  A  steady  air  flow  of  2.3  liters  per  minute  through 
the  device  was  maintained  -  at  least  for  the  duration  of  the  ammonia 
vapor  flow  over  the  rsks  -  by  connecting  the  imp inge-rs  outlets 
to  vacuum  tanks. 

In  addition  to  the  above  general  arrangement  of  apparatus  and 
Instruments,  the  following  modifications  ware  ms-la  to  perform  the 
surface  and  underwater  releases  of  liquid  ammonia. 

k . 2, 2.1  Surface  Spills 

Both  one-gallon  and  five-gallon  surface  spills  were  performed. 

The  spill  bucket  and  tipping  arrangement  used  in  the  ona-gallcn  tests 
is  shown  in  Figure  4-1.  This  arrangement  was  identical  to  that 
used  in  the  laboratory  tests,  (See  Figure  3-1.)  A  special  bucket, 
shown  in  Figure  4-4,  was  fabricated  for  the  five-gallon  surface  spill 
tests.-  The  bucket  had  a  faucet  to  facilitate  slow  release  of  LNH^ 
(continuous  spills). 

As  shown  in  Figure  4-5a,  three  liquid  temperatures  and  two  vapor 
temperatures  were  measured  approximately  two  feet  downstream  of  the 

5v. 


e 


To  Critics! 

Flow  Orifice 
and  Vacuum  Tank 


1 


Vaoor 


Boric  Acid  Solution 


FIGURE  4-3  SKETCH  OF  IMPINGER  DEVICE  FOR  VAPOR  SAMPLING 


spill  point.  Pool  temperatures  were  measured  at  1/8  inch,  1/2  inch, 
and  3/4  inch  below  the  water  surface.  Vapor  temperatures  were  measured 
at  about  one  inch  and  two  inches  above  the  water  surface  during  the 
one-gallon  spills,  and  at  three  inches  and  12  inches  above  the  surface 
during  the  five-gallon  spills.  In  addition,  five  vapor  temperatures 
were  measured  at  the  downstream  rake. 

4. 2. 2. 2  Underwater  Release 

Figure  4-5b  shows  the  thermocouplearrangement  used  during  the 
underwater  release  erperimants .  Figure  4-6  shows  the  special  release 
system  built  utilizing  a  one-gallon  pressure  vessel  (stainless-steel 
fire  extinguisher),  a  3/8-inch  stainless-steel  ball  release  valve,  and 
0.. 5-inc.h  ID  stainless-steel  submerged  release  tube.  A  similar  system 
employing  a  smaller  release  valve  and  a  0.364-incn  ID  submerged 
release  tube  was  also  fa-ricated. 

The  pressure  vessel  was  filled  through  a  T-fitting  attached  to  the 
release  valve,  and  once  the  tank  was  filled  to  the  desired  liquid 
level,  the  supply  was  valved  off  at  the  I.  The  release  valve  was 
also  closed  to  allow  a  pressure  rise  in  the  tank.  At  about  15  psig, 
ammonia  was  released  through  the  submerged  tube  into  the  pool. 

The  depth  at  which  the  liquid  was  released  was  varied  by  raising 
or  lowering  the  submerged  pipe. 

Nine  underwater  thermocouples  were  placed  1.75-inch  apart  at 
the  spill  site,  and  were  monitored  during  the  underwater  reiea.  »s. 

One  thermocouple  above  the  water  surface  was  also  used.  (See  Figure  4-5b. 

For  ail  tests,  wind  velocity  and  wind  direction  were  measured, 
adjacent  to  the  pool,  with  an  anemometer  and  a  direction  vane.  How¬ 
ever,  only  graduated-dial  readouts  were  available,  rather  than  graphical 
strip  recordings.  Colored  movies  and  35  ms  stills  were  made  of  all 
tests.  Some  tests  were  recorded  on  infrared  film.  Underwater  color 
movies  were  made  of  certain  subsurface  tests. 


FIGURE  4— 6b  ARRANGEMENT  AND  POSITION  OF  THE  THERMOCOUPLE  BANK  FOR  WATER 
TEMPERATURE  MEASUREMENT  DURING  UNDERWATER  RELEASE  OF  LNH, 


1  GALLON  UNDERWATER  DISCHARGE 


Table  4-1  lists  the  techniques  and  numbers  of  parameters  measured 
during  these  tests. 

4.2.3  Procedures 

4. 2. 3.1  Surface  Spills 

After  the  pool  was  filled  with  water,  the  thermocouple  assembly 
was  placed  in  position,  and  the  lmpinger  rake  was  positioned  downwind 
along  the  mean  wind  direction.  The  vacuum  tanka  were  evacuated, 
and  liquid  ammonia  was  transferred  to  the  spill  bucket.  When  the  proper 
volume  of  liquid  had  been  transferred,  and  the  wind  was  in  the  proper 
direction,  the  test  countdown  was  begun. 

About  60  seconds  before  the  spill,  the  temperature  recorders 
were  started,  and  at  five  seconds  before  the  spill,  the  impinger 
suction  was  started  with  a  solenoid-operated  valve.  The  LNH^  in 
the  spill  bucket  was  emptied  onto  the  water  surface,  all  at  once  in 
the  case  of  instantaneous  spills,  and  slowly  through  the  faucet 
in  the  case  of  continuous  spills.  After  the  vapor  cloud  had  disappeared 
and  about  10  minutes  had  elapsed,  the  water  in  the  pool  was  mixed 
thoroughly  by  operating  the  water  recirculating  pump  for  about  15 
minutes.  Five  water  samples  were  taken,  each  from  a  different  point 
in  the  pool,  for  laboratory  analysis  of  dissolved  ammonia  content. 

The  ammonia  vapor  dissolved  in  the  impinger  boric  acid  solutions  was 
also  analyzed. 

4. 2.3. 2  Underwater  Release 

The  procedure  was  similar  to  that  used  for  surface  spills,  but 
LNH^  was  transferred  to  a  pressure  tank.  The  tank  and  contents  ware 
weighed  continuously  and  when  a  pre-set  weight  of  LNH^  (15  to  18  lb3) 
had  been  transferred,  the  transfer  was  terminated.  Because  the  pressure 
tank  was  uninsulated,  the  pressure  would  immediately  start  to  rise. 

When  the  pressure  inside  the  tank  reached  about  15  psig,  the  test  count¬ 
down  was  begun.  When  the  wind  was  oriented  correctly,  the  vessel's 
valve  was  opened  and  LHH^  was  released  underwater.  The  impingers  were 
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switched  on  at  the  same  time,  and  the  rest  of  the  procedure  was 
similar  to  that  described  for  surface  tests. 

4.3  SURFACE  TEST  FINDINGS 

4.3.1  One-Gallon  Surface  Spills 

Only  instantaneous  spills  were  tested  in  this  category.  The 
data  and  results  are  shown  in  Table  4-2.  The  average  partition  ratio 
is  0.5  (0.552  if  only  the  results  with  city  water  are  used)  compared 
to  the  0.735  obtained  with  laboratory  spills  of  comparable  spill 
quantity.  (Figure  4-8  shows  these  partition  ratios  plotted  against 
spill  quantity.  A  discussion  of  the  above  results  is  deferred 
to  Section  4.6.) 

The  temperature  data  obtained  in  the  tests  are  given  in  Table  4-3. 
Note  that  within  about  a  second  after  the  spill  the  temperature  recorded 
by  the  immersed  thermocouples  dropped  to  nearly  the  LNH^  temperature 
and  then  began  to  rise.  The  temperature  increase  (of  5°  to  10°F 
above  water  temperature)  continued  for  three  to  five  minutes.  For 
one  test,  22CF  above  water  temperature  was  reached.  The  vapor  phase 
thermocouples  followed  a  similar  pattern.  However,  in  many  cases 
the  maximum  temperatures  indicated  by  the  vapor  phase  thermocouples 
were  20°  to  40°F  above  the  air  temperature.  Such  values  are  considerably 
higher  than  the  increases  shown  by  the  water  thermocouples.  The 
thermocouples  on  the  instrument  rake  (15  feet  downwind)  recorded  ambient 
temperature  except  when  the  ammonia  cloud  passed  over  the  rake.  This 
occurred  in  many  tests,  and  the  temperature  in  the  cloud  was  indicated 
to  be  between  15°  and  20°F  lower  than  air  temperature.  Such  low  tempera¬ 
tures  lasted  for  only  about  one  to  two  seconds,  indicating  rapid 
passage  of  the  cloud  past  the  rake.  (Unfortunately,  the  recording 
instrument  to  which  rake  couples  T^,  T2,  and  T^  were  connected  did 
not  record  the  event  because  of  its  slow  response  time  and  noncontinuous 
printing.)  For  all  cases  when  the  cloud  passed  through  the  rake  the 
temperature  indicated  by  a  thermocouple  at  a  higher  elevation 
lover  than  the  temperature  indicated  by  the  one  beneath  it.  (This 
phenomenon  is  shown  better  for  t'ue  five-gallon  spills.) 
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Table  4-2 


One-Qalloa  Surface  Spill  Tessa 


Spill 

Spill 

Quantity 

Spill 

Partition 

Water 

Date 

No. 

<liters)+ 

Funnel 

Ratio 

Source 

6/1 

1 

3.25 

w 

... 

6/8 

1 

3.50 

- 

0. 30 

Pond 

6/11 

1 

2.50 

- 

0.45 

Pond 

6/14 

1 

3.50 

- 

0.51 

City  Water 

6/14 

2 

3.5L 

1/2" 

0,49 

City  Water 

6/14 

3 

3.50 

- 

0.65 

City  Water  ! 

6/15 

1 

3.50 

- 

0.56 

City  Water 

6/15 

2 

3,50 

- 

0.55 

City  Water 

mean  *  d . 5 


aaaa  for 

City  water  »  0,552 


Note:  1  gallon  =3.6  litres 
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used  for  Tj^  Tj  and  T-j  could  not 
register  thio  brief  event.  The 
water  tcop.  first  dips  momentarily 
to  59*?  end.  then  rise*. 


Figure  4-7  shows  plots  of  the  total  mass  of  the  NH^  vapor  collected 
in  each  impinger.  Figure  4-7a  shows  the  distributions  of  vapor  collected 
along  the  horizontal  beam  three  feet  off  the  ground,  and  Figure  4- 7b 
shews  vertical  distributions.  The  figure  also  shows  the  best  possible 
Gaussian  profiles  that  pass  through  the  experimental  points  in  each 
experiment . 

4.3.2  Five-Gallon  Surface  Spills 

A  set  of  eight  experiments  were  conducted  with  spills  of  about 
five  gallons;  six  were  instantaneous  spills,  and  two  were  continuous 
spills  at  65  ml/second.  Table  4-4  gives  the  details  of  the  tests  in 
this  category,  such  as  spill  quantities,  partition  ratio  measured, 
etc.  Note  that  the  average  partition  ratio  for  an  instantaneous  spill 
is  less  than  that  for  a  continuous  spill.  Figure  4-8  shows  these 
partition  ratios,  together  with  those  obtained  for  the  one-gallon 
spills . 

Table  4-5  shows  the  temperature  data  obtained  for  the  five-gallon 
tests.  Note  that  the  vapor  phase  thermocouples  were  three  inches  and 
12  inches  above  the  water  level.  The  temperature  data  are  very  similar 
to  those  obtained  for  one-gallon  spills.  The  large  increase  in  vapor- 
phase  temperature,  compared  to  a  relatively  small  change  in  water- 
temperature,  can  be  clearly  observed. 

Figure  4-9  is  a  plot  of  both  vapor-concentration  impinger  data 
and  best-fit  Gaussian  curves.  Note  that  for  the  vertical  distribu¬ 
tion  the  maximum  concentration  occurred,  in  two  experiments,  at  an 
elevated  level  rather  than  at  ground  level.  This  corresponds  well 
with  the  temperatures  recorded  by  the  vertical  array  of  thermocouples. 

The  lowest  temperature  was  recorded  by  the  10-foot-high  thermocouple 
during  passage  of  a  cloud.  Tbi.  seems  tc  indicate  that  NH^  vapor 
rises  rapidly  in  the  atmosphere  during  dispersion. 

4.4  UNDERWATER  RELEASES 

Seven  underwater  tests  were  conducted,  each  releasing  about  one 
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One-Gallon  Spill 


Experiment 

Date  Symbol 


-5-4-3  -2  -1  0  1  2  3  4  5 


Feet 


FIGURE  4— 7a  HORIZONTAL  ARRAY  ASPIRATOR  DATA  {D,  E,  F,  G! 
13  Ft.  Above  Grcnnd) 


One-Gallon  Spill 


C 


«- 

D 


r  B 


777777 


F 


G 


Experiment 

Date 

Symbol 

6/7 

°i 

a'n 

°  ! 

6/14 

O 

6/14 

A 

6/14 

tv 

6/15 

O 

6/15 

O 

Darkened  symbols  are 
estimated  data. 

C 


8 


Feet 

ight  Above  Ground 


Curve; 
due  t; 
data. 


Table  4-4 


Small  Scale  Tests  Above  Water 

Spill  Avg.  Spill  Wind 

Quantity  Partition  Rate  in  Speed 


Test 

(litres) 

Ratio 

ml/sec 

in  moh 

I2£e 

7/1S-1 

17.55 

0.592 

Instantaneous 

5-6 

Above 

7/19-2 

21.45 

0.707 

Instantaneous 

5-6 

Above 

7/20-1 

19.50 

0.637 

61.5 

0-5 

3/8"  pipe 

7/20-2 

21.45 

G.686 

61.5 

0-5 

3/8"  pipe 

7/23-1 

17.06 

0.543 

Instantaneous 

10 

Above 

7/23-2 

19.50 

0.526 

Instantaneous 

10 

Above 

8/3-1 

21.00 

0.597 

Instantaneous 

0-5 

Above 

8/3-2 

21.00 

0.430 

Instantaneous 

0-1 

Above 

mean  *■ 

0.59 

Mean  for  instantaneous 

spill  only  ** 

0.566 

Mean  for  continuous 

spill  only  “ 

0.662 

67 
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FIGURE  4-8  PLOT  OF  PARTITION  RATIO  VERSUS  OUANTITY 
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gallon  of  LNH^  at  depths  between  14  and  18  Inches  below  the  water 
surface.  Table  4-6  gives  the  details  of  the  testa.  During  the 
process  of  release,  only  minimal  agitation  of  the  water  surface  was 
noticed,  except  when  the  release  was  at  seven  inches.  In  that  case 
the  water  surface  was  punctured  with  mini  fountains  (four  to  five 
inches  high),  which  were  almost  directly  over  the  release  point. 

Some  "crackling"  sound  was  recorded  and  the  ground  very  close  to  the 
edge  of  the  pool  rattled,  as  if  there  were  a  small  earthquake. 

The  recordings  of  the  crackling  noises,  made  with  a  hydrophone  in 
the  pool,  have  not  provided  any  useful  information. 

Note  that  the  measured  partition  ratios  are  all  very  high, 
indicating  almost  total  dissolution  of  the  LNH^  in  water.  In 
fact,  in  the  first  two  tests,  he  impingers  did  not  register  any 
ammonia  vapor  at  all.  Therefore,  in  subsequent  tests,  no  impinger 
data  were  taken.  The  low  values  of  the  partition  ratio  obtained  for 
the  first  two  tests  were  probably  due  to  incorrect  sampling  of 
the  water  after  the  test.  (The  pool  water  was  not  mixed  well  before 
samples  were  obtained  for  those  tests.) 

The  water  temperature  was  measured  using  a  bank  of  thermocouples 
along  a  vertical  line  about  five  inches  from  the  exit  of  the  underwater 
discharge  pipe.  (See  Figure  4-4b.)  Figure  4-10  shows  typical 
records  of  the  temperatures  for  14-inch  release  depth  and  7-inch 
release  depth.  These  indicate  that,  close  to  the  point  of  release,  large 
temperature  variations  occurred.  Although  very  low  temperatures  were 
recorded  close  to  the  release  point,  the  high  temperatures  occurred 
only  above  the  release  point.  Also,  the  maximum  increase  in  temperature 
at  any  point  seemed  to  be  not  raor  than  about  10°F,  while  up  to  30°P 
drops  have  been  seen.  The  latter  temperatures  may  be  due  to  cold  LNB^ 
being  released  in  the  vicinity  of  the  thermocouple.  No  particular 
pattern  can  be  ascribed  to  the  temperature  variations  with  time. 

However,  there  is  a  definite  decrease  in  the  effect  as  one  goes  farther 
firm  the  release  point.  The  thermocouple  in  the  air  recorded  a  more 
disturbed  temperature  profile  that  varied  as  much  as  15°F. 

The  spread  radius  was  estimated  from  movies. 
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4.5  ANALYSIS 


From  the  impinger  data,  we  can  calculate  the  total  vapor  mass 


as  follows: 


m(y,z)  “  moles  of  vapor  collected  at  the  location  y,z  on  the 
instrument  rake  at  distance  x 


where 


y  »=  direction  normal  to  the  wind  direction  and  parallel 
to  the  ground 

z  =  vertical  direction 


A  »  Area  under  the  measured  y-molar  distribution  =1  ®(y,z)  dy 

y  3  J  -  « 


Area  under  the  measured  z-molar  distribution  *= 


m(y,z)  dz 


Assuming  that  the  distributions  in  other  regions  of  the  y-z  plane 
are  similar  to  the  measured  distributions,  it  can  be  easily  shown  that 
the  total  mass  of  vapor  liberated  is 


A  A 

M  =  li  —  -P-  V 

v  q  m(yo,Zo) 


<&-!> 


U  15  molecular  weight  of  NH^ 
q  ■>  constant  suction  rate 
U  a  wind  velocity 

In  Eq.  4-1,  m(yQ,zo)  represents  the  value  of  moles  of 
vapor  cc llected  by  the  impinger  at  the  common  point  between  the 
vertical  and  horizontal  arrays,  Once  the  mass  of  vapor  is  calculated, 
£hen  the  partition  ratio  can  be  calculated  with 


F  «  1  - 


(4-2) 
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Table  4-7  shows  the  calculated  values  of  the  partition  ratio, 
made  using  the  impinger  data  of  selected  tests  by  the  above  method. 

These  are  compared  with  the  partition  ratio  values  obtained  from 
water  sampling.  Note  that  for  the  tests  where  the  cloud  hit  the 
instrument  rake  squarely,  the  mass  of  vapor  liberated  in  the  spill 
could  be  estimated  fairly  accurately.  A  comparison  of  the  values 
of  the  partition  ratios  both  calculated  and  measured  indicates  that 
in  most  cases  they  agree  reasonably  well. 

Note  that  the  mass  of  vapor  calculated  by  Eq.  4-1  is  quite  sensitive 
to  wind  velocity.  During  the  pool  experiments,  wind  velocity  was  not 
measured  accurately.  The  numbers  for  wind  velocity  given  in  Table  4-7 
are  best  estimates,  based  on  readings  from  a  meter.  It  is,  however, 
encouraging  to  note  that  even  if  the  wind  velocity  was  actually 
slightly  different  from  the  above-quoted  values,  the  mass  of  vapor 
estimated  would  be  correct  within  a  factor  of  2  for  the  worst  case. 

Table  4-8  indicates  the  same  vapor  concentration  data  in  terms 
of  the  standard  deviations  of  the  best  Gaussian  profiles  that  fit  the 
distribution  data.  The  mean  values  of  the  horizontal  standard  devi¬ 
ations  (o  )  are  1.7  feet  and  3.7  feet  for  the  one-gallon  and  five-gallon 
spills,  respectively.  Similarly,  the  values  in  the  vertical  direction 

(a  )  are  2.2  and  2.9  feet  respectively.  These  values  are  compared  with 
z  __ 

the  values  of  o's  estimated  from  Pasquill-Gifford  curves  for  a  distance 

of  15  feet  from  the  source.  (See  Figure  6-17.)  Table  4-9  lists 

both  of  these  values.  The  agreement  is  reasonably  good.  This  indicates 

that  probably  the  regular  vapor  dispersion  models  could  be  used  to 

describe  the  NH^  vapor  dispersion.  (This  subject  is  discussed  better 

in  Section  5.5  in  the  light  of  large-scale  experimental  data,) 

No  systematic  and  meaningful  analysis  of  the  temperature  date  could 
be  made.  A  general  discussion  of  the  implications  of  the  results  of 
temperature  measurements  is  given  in  the  next  section, 

4.6  DISCUSSION 

There  is  a  definite  pattern  to  the  partition  function  results. 


The  partition  function  is  high,  close  to  0.9,  when  the  LHH^  is 
released  underwater  and  is  low,  about  0.57,  when  a  large  quantity 
is  released  instantaneously  on  the  water  surface.  When  the  release 
is  continuous  and  on  the  water  surface,  the  partition  function 
ranges  from  0.66  to  0.7.  These  results  conform  to  what  one  might 
expect,  based  on  the  physics  of  the  release  problem. 

In  the  case  of  underwater  release,  the  LNK^  and  the  vapor  generated 
have  enough  tine  to  completely  dissolve  in  water.  The  vapor  that 
is  generated  at  the  point  of  release  dissolved  during  its  travel  upward 
through  the  water.  Hence,  the  partition  function  is  high. 

In  the  case  of  a  large  instantaneous  release  on  the  surface, 
two  phenomena  may  be  occurring  simultaneously.  One  is  the  boiling 
of  saturated  LNH^  to  produce  saturated  ammonia  vapor,  while  the 
LNH^  and  water  reaction  also  contributes  vapor.  The  other  phenomenon 
is  the  aerosol  formation,  as  fine  drops  of  ammonia  are  thrown  up 
into  the  air  by  the  brief  but  violent  agitation  that  occurs  during 
the  surface  boiling  process.  Because  of  this,  the  total  quantity  of 
liquid  that  goes  into  solution  with  water  is  low  and  Che  partition 
ratio  is  small.  We  observed  that  the  mean  partition  ratio  of  the 
surface  spills  was  quite  below  that  obtained  in  the  laboratory  experi¬ 
ments.  We  speculate  that  in  the  laboratory  experiments  the  confining 
walls  of  Che  test  tank  may  have  allowed  liquid  ammonia  droplets  to 
condense  from  the  vapor  and  return  to  the  liquid,  resulting  in  an 
increased  partition  ratio. 

The  very  low  partition  ratios  obtained  when  pend  water  was  used 
in  the  pool  may  be  due  to  the  initial  high  pH  cf  the  water  (resulting 
in  love’-  amounts  of  LKH^  in  solution) .  However,  because  the  pond 
water  was  not  analyzed  before  the  tests,  it  is  hard  to  give  any  definitive 
reasons  for  the  observed  low  values.  In  the  continuous  release  case, 
the  partition  ratio  value  is  between  the  underwater  and  surface  spill 
values,  for  obvious  reasons.  In  fact,  the  partition  ratio  for  the 
continuous  release  is  closer  to  the  predicted  value  (see  Section  6,2) 
of  0.73  based  on  the  thermodynamic  mixing  model. 


Immediately  after  a  surface  spill,  the  thermocouples  in  the  water 
showed  a  marked  temperature  decrease,  with  the  one  nearest  the  surface 
indicating  the  lowest  temperature.  Within  15  seconds,  the  tempera t-uic. 
was  everywhere  above  the  ambient  and  slowly  rising.  A  maximum 
temperature  was  reached  in  two  to  three  minutes,  with  the  top  thermocouple 
again  showing  the  maximum  value.  In  many  experiments,  the  maximum 
was  about  120°F,  a  value  very  close  to  the  maximum  (122°F)  predicted 
by  the  second  adiabatic  mixing  model  (see  Section  6.2).  The  observable 
reaction  was  complete  within  about  5  to  10  seconds,  whereas  the  rise 
in  water  temperature  continued  for  about  two  to  three  minutes  and 
stayed  at  the  high  temperature  for  over  five  minutes  before  starting 
to  decrease.  This  slow  increase  in  water  temperature  is  probably  the 
result  of  a  slow  reaction  of .concentrated  NH^OH  solution  at  the  water 
surface  with  the  layers  of  pure  water  below.  Once  the  boiling  had 
stopped,  the  mixing  was  brought  about  by  molecular  diffusion, 
and  this  is  a  long-term  process.  The  persistence  of  the  high  water 
temperature  at  the  surface  is  due  to  the  slowness  of  natural  cooling. 

It  is  interesting  to  see  that  not  much  lateral  mixing  takes  place, 
even  though  hot  water  (or  the  NH4OK  solution)  is  less  dense  than  the 
surrounding  water. 

The  temperature  rise  recorded  by  the  thermocouples  in  the  vapor 
phase,  seems  to  be  more  or  less  in  phase  with  water  thermocouple  results. 
Initially  there  was  a  sudden  decrease  in  temperature  -  to  almost  -27°F 
-  and  then  a  steady  rise.  The  low  temperature  was  indicated  only  for 
about  1  to  1.5  seconds.  This  is  probably  the  result  of  the  passage 
of  NH^  vapor  past  the  thermocouple.  However,  the  maximum  temperature 
registered  by  these  thermocouples  is  far  above  the  value  registered  by 
the  water  thermocouples,  though  the  maximum  readings  occur  at  about  the 
same  time  and  persist  as  long.  Also,  in  many  cases  the  thermocouple 
farther  away  from  the  water  surface  (two  inches  above  or  12  inches 
above)  indicated  higher  temperatures  than  the  one  nearer  to  the  water. 

The  reasons  for  this  behavior  are  nob  known  precisely,  but  there  are 
various  possibilities.  One  possibility  is  that  a  small  drop  of 
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water,  or  concentrated  NH^OH  solution  thrown  up  into  the  air  by  the 
boiling  process  lodged  itself  on  the  thermocouples  and  the  reaction 
between  the  various  radicals  was  completed  on  the  thermocouple  bead. 

Though  the  temperature  rise  can  be  explained  by  this  ’'model,"  the 
duration  of  temperature  rise  cannot  be  explained  because  LNH^  and  water 
react  rapidly  any  time  they  are  mixed  together.  A  more  likely  possibility 
is  that  of  a  reaction  between  a  droplet  of  water  on  the  thermocouple 
bead  and  a  continuous  stream  of  vapor  that  is  being  liberated  from 
the  water  surface  during  the  slow  reaction  process  in  the  water  mass. 

In  this  case,  one  cannot  only  explain  the  high  temperatures  recorded 
(adiabatic  reaction  of  NH^  vapor  dissolving  in  water  results  in  very 
high  mixture  temperatures)  but  also  the  duration  to  attain  these 
maximum  values.  It  is,  however,  noteworthy  that  in  all  experiments, 
except  one,  the  maximum  temperature  rise  above  the  ambient  value  recorded 
by  the  vapor  phase  thermocouples  was  never  more  than  52C'F.  This 
value  is  the  maximum  rise  possible  by  the  adiabatic  mixing  of  LNH^  and 
water  (see  Section  6.2,  Model  2).  However,  for  the  adiabatic  dissolu¬ 
tion  of  gaseous  ammonia  (GNH^)  in  water,  theoretically  the  mixture 
temperature  can  almost  reach  the  temperature  of  boiling  water  (2J.2°F). 
Though  the  exact  nature  of  the  reaction  that  brings  about  such  a  high 
temperature  reading  of  the  vapor  phase  thermocouple  is  not  known,  we 
experimentally  verified  that  is  is  no£  due  to  any  ionization  phenomenon 
caused  by  charged  vapor  flow  over  the  dry  thermocouple  bead. 

The  temperatures  indicated  by  the  thermocouples  or.  the  instrument 
rake  are  useful  to  the  extent  that  they  indicate  the  passage  of  the 
cloud  past  the  rake.  The  entire  cloud  passed  the  rake  within  a  matter 
of  two  seconds.  (Because  one  of  the  recorders  used  for  measuring 
these  temperatures  was  of  the  Integrating  type,  with  duration  of 
integration  considerably  longer  than  two  seconds,  no  temperature  change 
was  registered  by  the  printout  of  this  recorder.)  For  many  tests, 
the  temperature  indicated  by  the  highest  thermocouple  was  less  than 
the  temperature  read  by  the  couple  beneath  it.  In  the  same  tests 
the  vapor-concenti-ation  impinger  data  also  indicated  a  maximum  at 
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the  top  impinger.  There  results  Indicate  that  the  center  of  the  cloud 
was  far  above  the  ground  surface  and  was  still  fairly  cold  (indicating 
the  presence  of  partially  mixed  GNH^).  It  also  indicates  that  the 
cloud  rises  rapidly  during  the  process  of  dispersion. 

The  relative  agreement  between  the  vapor  mass  flow  estimated 
from  the  impinger  data  and  that  obtained  from  water  sampling  indicate 
that  it  is  possible  to  estimate  the  partition  ratio  from  measure¬ 
ments  of  vapor  concentration  in  air.  However,  note  that  the  calcula¬ 
tion  of  the  partition  ratio  is  very  sensitive  to  the  accuracy  of  the 
wind  velocity  data  and  that  the  instrument  rake  be  hit  squarely  by 
the  cloud.  Therefore,  while  in  principle  it  is  possible  to  obtain  the 
mass  of  vapor  flow  by  time-integrated  concentration  measurements , 
cooperation  by  nature  is  very  essential  to  assure  success  in  thxs 
method.  Results  shown  in  Section  4.5  also  indicate  that  the  standard 
deviations  (or  dispersion  parameters)  for  the  measured  concentrations 
could  he  reasonably  predicted  from  existing  methods  in  the  literature, 

A  complete  discussion  oa  the  dispersion  results  Is  given  in  Section  5. 

The  temperature  results  of  the  underwater  tests  indicate  that 
the  maximum  temperature  change  occurs  near  the  exit  section.  The 
maximum  increase  in  water  tempera'.ure  is  not  more  than  about  I0°F, 
whereas  no  temperatures  were  less  than  30°F  below  water  temperature. 
Since  no  details  of  the  SH^OH  concentration  could  be  measured, 
the  temperature  data  by  itself  does  not  give  much  quantitative  informa¬ 
tion.  However,  the  partition  ratio  was  not  very  affected  by  the  release 
at  7  inches  or  14  inches  depth.  In  both  cases  the  release 
is  at  a  depth  greater  than  10  times  the  diameter  of  the  discharge 
pipe.  The  result  is  in  qualitative  agreement  with  the  result  of  the 
simple  order-of-raagnitude  theoretical  analysis  presented  in  Section  6.$. 
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5 ,  LARGE-SCALE  EXPERIMENTS 


5.1  PURPOSE  AND  SCOPE 

The  purpose  of  the  large-scale  field  test  program  was  to  obtain 
downwind  ammonia-vapor  concentration  data  from  as  large  a  test  as 
practical  in  an  unconfined,  natural  environment.  It  was  hoped  that 
such  data  (in  addition  to  that  collected  from  the  .intermediate  size 
experiments)  would  allow  extrapolation  to,  and  prediction  of,  the 
downwind  vapor  concentrations  to  be  expected  from  releases  of  LNH^  amounts 
ranging  up  to  and  including  3,000  tons. 

These  experiments  were  based  on  the  premise  that  a  certain  portion 
of  the  NH^  vapor  cloud  might  dissolve  in  water  as  the  cloud  traveled 
considerable  distances  downwind  over  open  stretches  of  water.  In 
addition,  it  was  assumed  that  the  vapor-sampling  devices  (which  provided 
adequate  data  for  indirectly  calculating  the  partition  ratio  for  the 
intermediate-scale  experiments)  could  also  be  relied  on  to  provide 
adequate  data  for  large-scale  tests. 

5.2  APPARATUS  AND  PROCEDURES 
5.2.1  Test  Facility 

The  large-scale  experimental  program  was  conducted  on  a  pond, 
located  at  Space  Research  Corporation  headquarters,  along  the  United 
States-Canada  border.  Seven  8-foot  by  8-foot  instrumentation  rafts 
and  an  8-foot  by  8-foot  spill  platform,  with  25-foot  by  4-foot  walkway* 
were  constructed  on  the  pond  and  anchored  in  the  configuration  shown 
in  Figure  5-1.  A  partial  side  view  of  spill  platform  and  iaft  A  are 
shown  schematically  in  Figure  5-2.  The  instrumentation  v-£'ts  were 
outfittad  with  cross-shaped  instrument  rakes  that  supported  eight 
impinger  devices,  as  shown  in  Figure  5-3.  As  for  the  intermediate-scale 
tests,  the  imuingers  on  any  one  rake  were  connected  to  a  set  of  vacuum 
tanks  with  vacuum  tubing  containing  in  line  limiting  orifices.  The 
use  of  solenoid-operated  valves  and  underwater  power  connections  to 
each  raft  allowed  the  entire  system  to  be  remotely  activated  from 
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FIGURE  5-2  PARTIAL  SIDE  VIEW  OF  THE  SPILL  PLATFORM  AND  RAFTS 


shore  just  prior  to  an  LNF.^  release, 

5.2.2  Instrumentation 

Figure  5-3  shows  the  details  of  the  vapor  sampling  instrumentation 
set-up  for  the  seven  rafts  (designated  rafts  A~G  as  in  Figure  5-1). 

In  addition  to  this  arrangement  of  vapor  sampling  apparatus,  the 
following  additions  and  modifications  of  instrumentation  were  made 
to  perform  the  surface  and  underwater  releases  of  LNH^* 

5. 2. 2.1  Surface 

In  all,  ten  surface  releases  of  LNH^  were  performed.  Three 
five-gallon  spills  were  made  to  correlate  the  data  with  that  obtained 
from  the  intermediate-scale  tests,  and  seven  were  50  gallon 
spills,  made  to  provide  additional  data  for  this  phase  of  the  project. 

The  five  gallon  spill  tests  used  the  specially  fabricated  bucket 
from  the  intermediate-scale  tests.  The  details  cf  this  bucket  and 
its  associated  mechanisms  were  shown  in  Figure  4-4. 

The  50  gallon  spill  tests  wera  conducted  with  a  similar  apparatus, 
using  a  55  gallon-capacity  insulated  steel  drum.  A  support  rope  held 
the  drum  during  filling,  and  a  rubberized  cable  was  added  to  lessen 

the  shock  to  the  drum  during  the  sudden  tipping  3nd  dumping  of  the 

contents.  The  floating  foam  block  previously  used  to  indicate  the 
volume  of  fluid  in  the  spill  vessel  was  replaced  by  a  yardstick,  so  that 

the  depth  cf  liquid  could  be  noted  just  prior  to  a  spill. 

In  addition  to  the  vapor-concentration  measurements,  five  liquid- 
temperature  measurements  and  two  vapor-temperature  measurements  were 
made  during  some  of  the  experiments,  liquid  temperatures  were  measured 
0.5  to  1.0  inch  below  the  water  surface.  Vapor  temperatures  were 
measured  six  feet  above  the  water  surface  approximately  60  feet  and 
117  feet  downwind  of  the  edge  of  the  spill  platform.  Figure  5-4  is 
a  schematic  of  the  underwater  thermocouple  position  during  the  surface 
release  tests. 
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FIGURE  5-4  PLACEMENT  OF  WATER  TEMPERATURE  THERMOCOUPLES  FOR  LARGE-SCALE  TESTS 


5 .  ? .  2 . 2  Underwater 


Four  underwater  releases  of  five  gallons  of  ammonia  were  made 
with  a  specially  fabricated  2.5-cubic-foot  capacity  welded-steel 
pressure  vessel  and  appropriate  plumbing.  Figure  5-5  shows  the 
layout  of  the  underwater  release  system,  including  a  1.5- inch 
ID  PVC  pipe  connected  to  a  1.5-to-O. 5-inch  reducer,  to  which  a  length 
of  0.5-ir;:h  ID  pipe  is  attached  at  a  depth  of  five  feet.  Similar 
systems  with  discharge  depths  of  three  or  five  feet  were  also  fabricated 
with  pipes  of  constant  1.5  inch  diameter. 

The  spill  vessel  was  filled  by  placing  it,  minus  external  plumbing, 
on  a  scale  and  adding  LNHg  until  the  desired  weight  was  reached.  The 
vessel  was  then  positioned  and  connected  to  the  piping  system  and 
allowed  to  warm  until  the  internal  pressure  had  increased  to  about 
20  psig.  Then,  the  ball  valve  was  remotely  opened  and  LNH^  was 
released. 

5.2. 2.3  Other  Instruments 

Wind  velocity  and  direction  were  measured  during  all  tests  with 
an  accurate,  sensitive  anemometer  and  direction-vane  sensor  system, 
calibrated  and  connected  to  a  high-speed  strip-chart  recorder.  Ail 
but  two  spills  were  recorded  with  16  mm  color  movies  (during  one 
test  the  camera  power  system  malfunctioned,  and  during  ancther  the 
camera  jammed).  T.n  addition,  35  mm  color  stills  were  taken.  Table 
5-1  lists  the  test  measurements. 

5.2.3  Procedures 

The  test  procedures  for  both  surface  and  underwater  releases 
were  similar  to  those  used  in  the  intermediate  size  experiments. 

Only  minor  changes  -  dictated  by  the  location  and  size  of  the 
experiments  -  were  made.  These  primarily  were  because  the  pond  wafer 
could  not  be  thoroughly  mixed  and  sampled  as  the  swimming  pool  water 
had  been. 
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Vi-Inch  Fill  Lina  and  Valve* 


Table  5-1 


^^Rfr-Scale  Test  Measurements 


Surface 

Tests 


Concentrations 

Downstream  Vapor  Samples  8-56 

Temperatures 

Liquid  5 

Downstream  Vapor  2 

Sound  1 


Photography 

Color  Movies 
Infrared  Movies 
35  mm  Stills 

Wind 


Underwater 

Tests* 


8 


5 

1 


Velocity 

Direction 


*dot  all  measurements  were  made  for  every  test. 


5.3  TEST  RESULTS 


5.3.1  Surface  Spills 

In  all,  three  five  gallon  and  seven  50  gallon  instantaneous  spill 
tests  were  conducted  in  this  series.  Table  5-2  summarizes  the  test 
data,  ambient  conditions,  and  vapor  cloud  behavior  during  the  tests. 

Wind  velocity  data  were  obtained  by  timing  the  vapor  cloud  movement 
in  the  movies,  as  also  from  the  chart  records.  Table  5-3  shows  the 
wind  characteristics  obtained  from  strip-chart  records  of  wind  velocity 
and  direction  (measured  by  the  anemometer  and  the  direction  vane) . 

The  table  indicates  the  variable  character  of  the  wind. 

Figure  5-6  shows  typical  water-temperature  data  obtained  during 
a  50  gallon  test.  The  water  temperature.  wa3  considerably  higher 
than  was  observed  during  the  intermediate-scale  tests.  However,  during 
the  five  gallon  spills  of  the  large-scale  tests,  the  maximum  tempera¬ 
ture  increase  over  the  initial  water  temperature  recorded  by  the  came 
set  of  thermocouples  was  about  5°  to  6°F.  Note,  however,  that  the 
line  of  the  thermocouples  for  the  latter  case  was  about  five  feet 
from  the  center  of  the  five-gallon  spill  bucket.  Therefore,  it  is  likely 
that  the  thermocouples  were  not  within  the  boiling  zone. 

The  radii  of  the  boiling  zones  were  measured  primarily  by  scaling 
from  the  movies.  Figure  5-7  shows  a  plot  of  boiling  zone  radius 
as  a  function  of  the  spill  quantity.  Because  of  the  difficulty 
of  exactly  locating  the  edges  of  the  boiling  region  (due  to  vapor-cloud 
overlap)  a  range  of  values  has  been  given.  The  plot  shows  that  the 
estimate  of  a  five  gallon  spill  spread  from  the  large-scale  tests 
does  not  seem  to  agree  with  that  from  the  intermediate-scale  size  tests. 

The  impinger  data  from  several  typical  tests  (5,6,8,  and  13) 
are  shown  plotted  in  Figure  5-8  through  5-10.  These  represent  the 
data  from  50  gallon  (nominal)  spills  obtained  under  both  very  high  and 
very  low  wind  velocity  conditions.  Figures  5-8a  and  5-8b  show  the 
horizontal  and  vertical  distributions  of  the  vapor  mass  collected  at 
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FIGURE  B— 0  TEMPERATURE  VARIATION  OF  WATER  IN  THE  SPILL  ZONE 

FOR  A  50  GALLON  SPILL 


Moles  of  Vapor  Collected 


Cross  Raft  Distance  -  Fe<rt 


FIGURE  S-Ss  HORIZONTAL  IWPINSER  DATA  ON  THE  FIRST  ROW  RAFT  A 


F9GURE  5-8*  HORIZONTAL  IMPINOER  DATA  ON  THE  SECOND  ROW  OF  RAFTS 
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raft  A.  (The  figures  also  jhow  the  best-fit  Gaussian  curves, 
wherever  possible.)  Similarly,  the  other  figures  show  the  horizontal 
distributions  for  the  second  and  third  row  of  rafts.  Figures  5-9a 
and  5-9b  show  the  vertical  distribution  data  just  of  those  rafts  in 
the  second  and  third  row  through  which  the  vapor  cloud  passed.  The 
figures  indicate  that  for  the  tests  In  which  the  wind  velocity  was 
high  (tests  8  and  13)  the  horizontal  distributions  can  be  reasonably 
represented  by  Guassian  profiles.  For  low  wind  velocity  (tests  5  and 
6)  the  distribution  seems  to  be  quite  different  from  Gaussian,  and 
tends  to  be  more  uniform  for  the  rafts  at  the  same  distance  from  the 
source.  Also,  in  general,  the  concentrations  of  NH^  vapor  collected 
were  much  lower  for  the  tests  with  low  wind  that  those  with  high 
wind.  The  vertical  distributions  show  that  during  high  wind  veloci¬ 
ties  the  peak  concentration  occurred  at  a  lower  altitude  than  during 
low  wind.  However,  non-unifcra  distributions  also  occurred,  especially 
at  low  wind  speeds.  (See  test  5  in  Figure  5-8b.) 

It  is  clear  from  the  vertical  concentration  distribution  data, 
obtained  during  the  present  series  of  experiments,  that  the  ammonia 
vapor  plume  rises  above  the  water  surface  as  it  disperses  downwind. 

Thin  has  also  been  confirmed  by  observation  of  the  visible  vapor 
cloud  movement  in  the  test  movies.  Tne  rate  of  rise  of  the  plume 
also  seems  to  depend  on  the  wind  velocity  and  the  size  of  the  spill. 
Figure  5-lla  shows  the  positions  of  peak  concentration  values  obtained 
from  the  vertical  distribution  of  Raft  A  impinger  data,  plotted  against 
wind  velocity.  (The  values  from  other  rafts  in  the  second  and  third 
row  have  not  been  plotted  because  of  the  uncertainties  in  their  "peak” 
values.)  Tne  theoretical  predictions  based  on  plume  theory  are  also 
shown.  (The  details  of  this  theory  are  given  in  Section  6. 7.  Figure 
5-lla  suggests  that  to  explain  the  observed  variation  of  the  peak 
concentration  position  with  height,  allowance  has  to  be  cade  for  the 
possible  presence  of  aerosols  in  the  vapor  cloud.  However,  there  is 
considerable  data  scatter.  The  cloud  motions  were  observed  from  the 
movies  and  the  motion  of  esch  puff  center  over  the  different  rafts 
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Heic  ht  of  Plume  Max.  Concentration  Above  Water  (Feet) 


was  noted.  Figure  5-llb  shows  the  heights  of  the  cloud  centers 
over  the  three  rows  of  rafts.  The  theoretical  predictions  based  on 
plume  theory  are  shown  on  the  same  plot.  (A  discussion  on  the  results 
is  given  in  Section  5.5.) 

Some  general  qualitative  observations  made  from  the  movies 
can  be  summarized  as  follows: 

•  In  quite  a  few  tests,  during  the  spread  of  LNH^  on  the  water 
surface  an  uneven  spread  could  be  observed;  that  is,  the 
spreading  pool  edge  was  not  a  uniform  convex  front,  but  had 
long  elongated  "fingers."  This  may  have  been  caused  by  che 
dynamics  of  the  spilling  process. 

•  The  cloud  formation  was  rapid,  and  in  all  cases  the  cloud 
hit  the  first  raft.  Characteristic  mushroom  cloud  forma¬ 
tion  was  seen  under  low  wind  conditions.  The  clouds  dispersed 
by  a  rapid  rise  at  the  top  of  the  mushroom,  followed  by  a 
stretching  of  the  stem  of  the  mushroom.  The  stretched 

cloud  seemed  to  pass  over  the  rafts. 

•  The  vapor  cloud  was  visible  for  a  longer  duration  (further 
distance  traveled  downwind)  when  the  atmospheric  humidity  was 
high. 

•  The  vapor  cloud  did  not  exhibit  any  special  affinity  for  water 
by  dispersing  close  to  its  surface.  In  fact,  the  tendency 
for  the  cloud  to  rise  in  the  atmosphere  was  very  evident. 

©  After  the  vapor  had  completely  disappeared,  a  distinctly 

visible  pool  spread  on  the  water  surface  at  the  spill  position. 
This  pool  was  distinguishable  from  the  surrounding  water  by 
the  lack  of  small  surface  waves  and  slightly  different  optical 
reflective  properties.  It  is  very  likely  that  this  represented 
the  spreading  NH^OH  solution. 


5.3.2  Underwater  Release 


Four  underwater  teats  were  conducted.  Approximately  five  gallons 
of  LNG^  were,  released  in  each  of  the  tests  about  four  feet  below  the 
water  surface.  The  temperature  in  the  water  was  measured  at  five 
locations.  Figure  5-5  shows  a  schematic  of  the  thermocouple  loca¬ 
tions  relative  to  the  mouth  of  the  discharge  pipe,  as  well  as  the 
general  arrangement  of  the  underwater  release  apparatus.)  Both  0.5-inch 
and  1.5-inch  discharge  pipes  were  used,  each  in  two  experiments. 

Table  5-4  lists  the  experimental  data,  together  with  some  observa¬ 
tional  remarks,  for  this  test  series.  Note  that  water  temperature 
could  increase  by  as  much  as  50*^  (though  in  one  case  it  was  about 
70°F),  and  that  for  the  smaller-diameter  pipe  discharge,  the  maximum 
temperature  observed  was  much  lower  than  when  the  discharge  was 
from  a  large  pipe  diameter.  This  may  indicate  that  the  rate  of  discharge 
ia  important  for  the  type  of  reaction  that  can  occur  between  water  and 
LNH^.  No  vapor  was  sensed  by  the  sniffers  in  the  first  raft,  nor 
did  observations  of  the  area  directly  above  the  release  point  indicate 
any  vapor  clouds.  Only  violent  bubbling  and  agitation  of  the  water 
in  these  regions  were  noticed.  The  thermocouples  close  to  the  release 
pipe  seemed  to  first  indicate  temperatures  as  low  as  0°  to  -3°F 
(probably  liquid  ammonia)  and  then  begin  to  indicate  increased  tem¬ 
peratures. 

5.4  ANALYSIS 

5.4.1  Partition  Ratio 

Partition  ratios  were  calculated  from  the  impinger  data  of  each 
raft  over  which  the  vapor  cloud  passed,  using  the  method  described 
in  Section  4.5,  Table  5-5  lists  the  calculated  partition  ratios  and 
dispersion  variance  parameter  values  for  the  impinger  data.  The 
latter  were  obtained  by  fitting  the  best  Gaussian  curves,  wherever 
possible,  to  the  impinger  data  for  the  horizontal  and  vertical  arrays. 

It  can  be  seen  from  the  table  that  the  partition  ratios  are  far  from 
being  unique.  In  fact,  in  some  cases  the  vapor  masses  calculated  from 
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the  impinger  data  are  very  small  compared  to  the  mass  of  liquid 
spilled.  In  general,  high  partition  ratios  were  obtained  for  those 
experiments  in  which  the  wind  velocity  was  very  low.  This  indicates 
that,  for  such  cases,  the  major  portion  of  the  cloud  missed  passing 
by  the  impingers  and  the  readings  obtained  by  the  impingers  were 
due  to  the  residual  vapor  that  passed  by  then.  (See  the  remark  about 
the  passage  of  the  stem  of  a  mushroom  cloud  in  Section  5.3.1.) 

Table  5-5  also  indicates  that  the  partition  ratio  values  calculated 
from  the  data  of  the  first  raft  seem  to  be  within  reasonable  limits 
(between  about  0.5  and  0.7  for  a  majority  of  the  experiments), 
but  no  definite  pattern  can  be  ascribed  to  the  data  of  the  second- 
and  third-row  rafts.  It  is  known  from  observational  data  that  in 
almost  all  the  tests  the  cloud  hit  the  first  raft*  almost  squarely, 
yet  did  not  necessarily  hit  the  rafts  in  the  second  and  third  rows. 

This  information,  coupled  with  the  large  differences  in  the  partition 
ratio  obtained  from  second-  and  third-row  raft  data,  is  a  further 
indication  of  the  rapid  rise  of  the  vapor  cloud  in  the  atmosphere. 

5.4.2  Impinger  Concentration  Distribution 

Most  of  the  horizontal  distributions  could  be  reasonably  fit  by 
Gaussian  profiles.  The  variances  obtained  from  these  profiles  are 
given  in  Table  5-5,  along  with  the  dispersion  parameter  (a)  values 
obtained  from  the  Gifford-Pasqulll  curves  (see  Figure  6-17)  for  the 
appropriate  atmospheric  type  and  the  downwind  raft  distance.  A 
comparison  of  the  measured  values  and  those  obtained  from  the 
Gifford-Pasquill  data  are  shown  in  Figure  5-12.  Note  that  the  measured 
ay'8  tend  to  be  slightly  higher  for  the  near  raft  and  smaller  for  the 
second  row  raft  than  the  Gifford  a  ’s.  For  larger  distances,  this 

y 

difference  seems  to  be  reduced.  Also  note  that  the  latter  values 
were  obtained  by  extrapolating  the  Gifford  curves,  which  have  o  values 
listed  only  for  distances  greater  than  100  meters.  Hence  allowance 
must  be  made  for  possible  errors  in  extrapolation.  Considering 
this  uncertainty,  the  agreement  seema  to  bs  quite  good. 
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In  the  case  of  vertical  distribution  of  the  measured  impinger 
data,  the  Gaussian  profile  seems  to  be  more  an  exception  than  the 
rule.  The  most  likely  explanation  for  this  behavior  is  tha  rapid 
rise  of  tha  cloud  at  low  wind  velocities.  In  fact,  note  that  in 
Table  5-5,  the  agreement  between  measured  and  Gifford  a,,  values  is 
good  only  for  tests  in  which  the  wind  velocity  was  high  (about 
10  mph). 

5.4.3  Ammonia  Vapor  Rise 

The  path  of  the  center  of  the  vapor  cloud,  liberated  by  the  LNH^ 
spill  cn  water,  can  be  described  by  the  theories  of  vapor  cloud 
dispersion  in  the  atmosphere  discussed  in  detail  in  Section  6.7. 
Essentially,  the  cloud  can  be  treated  either  as  a  maintained 
plume  or  as  a  puff  or  vapor.  Both  theories  are  applied  to  describe 
the  ammonia  vapor  cloud  dispersion,  and  their  results  are  then  compared 
with  the  experimental  observation. 

For  comparing  the  predictions  and  experimental  values  of  vapor 
cloud  rise,  the  following  experimental  values  are  used: 

M  **  mass  of  liquid  ammonia  spilled  =>  284.5  pounds 
p  «  partition  ratio 

R  *»  maximum  radius  of  spread  *»  12  feet 
max  r 

t  *  duration  of  evaporation  -  10  seconds 

T  »  ambient  temperature  (60°F) 

x  *  distances  to  the  rafts  =  21,  55,  and  110  feet 

g  “  acceleration  due  to  gravity  <=  32.2  feet/second 

The  physical  property  values  for  both  liquid  ammonia  and  the 
saturated  vapor  are  taken  from  Appendix  A.  The  following  atmospheric 
conditions  are  assumed  tor  the  sake  of  calculations: 


U  *  wind  velocity  »  8  mph  «  11.73  feet/second 
~  s-  potential  temperature  gradient  =  0.01  °K/mster 

3 

p  «*  density  of  air  •  0.07488  lba/ft 

p  ■  density  cf  saturated  ammonia  vapor  at  atmospheric  pressure 
v  ■=  .05556  Iba/ft-* 
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Bence ,  we  have,  assuming  no  aerosol  formation. 


d  “  1  -  ~  “  density  defect  *  0.258 
pa 

V  «  total  volume  of  vapor  generated  »  ]8  x  0.4  x  284.5  =  2050  ftJ 

5«4*3.i  Predictions  from  Plume  Theory  (See  Section  6. 7. 2.1.) 

Because  of  the  expanding  character  of  the  pool  source  of  ammonia 
vapor,  we  assume 

b  -  R  |  f~2 
i  max  /v] 

F  ■  buoyancy  flux  = 

fit 

£  «  buoyancy  length  =  — ■„  *  0.336  feet 
P  11. 73^ 

using  equation  6-56  we  have,  with  B«*l, 

2  -  0.795  x2'3. 

The  heights  of  the  vapor  cloud  predicted  by  this  theory  for  x  values 
of  21,  55,  and  110  feet  are  6,  11.5,  and  18.3  feet. 

1/2 

=  Brunt  Vaisalla  frequency  *  0.0185  rad  jans/second 

S  »  stratification  parameter  -  -ry.lla73  *  1891. 

• jju  x  • UXbj 

Hence,  from  equation  6-61 

za  o  maximum  plume  rise  height  »  0.336  x  1891  2/3  -  51.4  feet. 
Approximate  downwind  distance  within  which  this  height  is  reached  is 

x*  =  443  feet 

5. 4. 3. 2  Predictions  from  Puff  Theory  (Ses  Section  6. 7.2. 2) 

3i  -  buoyancy  volume  released  =  Vgd  =  17030  ft^/s^ 

g 

*c  -  buoyancy  length  -  ~  ■  11.12  ft 

Assuming  b  0.3  and  using  equation  6*62  ,  we  have 
1/2 
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z  *=  6.86  x 


The  heights  of  the  vapor  cloud  predicted  by  the  above  equation  by 
this  theory  for  x  values  of  21,  55,  and  110  feet  are  respectively 
31.5,  51,  and  72  feet. 

S  ■  stratification  parameter  -  *  57 

11.12  x  .OlSo 

Assuming  again  e  ■  0.093  and  a<  *  0.3,  we  have  from  equation  6.65 
z  «  273  feet 

CO 

This  height  is  reached  first  at  a  distance  x*  given  by  equation  6.66 

x*  «*  1790  feet. 

We  compared  the  values  predicted  by  the  above  two  theories  with 
the  experimental  values  as  shewn  in  Table  5-6.  The  values  in  the 
table  indicate  that  the  continuous  plume  theory  seems  to  correlate 
with  experimental  data  extremely  well.  The  experimental  data  plotted 
in  Figure  5-12b  (visual  observation  data)  indicate  that  height 
decreases  inversely  with  wind  velocity.  This  is  precisely  the  law 
for  the  continuous  plume  theory,  rather  than  the  puff  theory  in  which 
the  height  goes  down  as  the  inverse  square  root  of  the  velocity. 

Results  computed  for  other  wind  velocities  at  the  three  raft  positions, 
based  on  the  plume  theory,  are  shown  in  Figure  5-I2fc. 

Similar  calculations  of  cloud  rise  made  with  plume  theory, 
assuming  two  different  percentages  of  aerosol  in  vapor  (0Z  and  20Z) 
are  shown  in  Figure  5-12a  and  are  compared  with  those  obtained  from 
impinger  data.  The  comparison  indicates  that  maximum  concentration 
data  could  probably  be  explained  if  we  allow  for  aerosols  in  the  vapor 
cloud. 

5,4.4  Comparison  of  the  Vapor  Concentration  Data  with  Predicted  Values 

In  this  section,  we  show  a  method  of  calculating  the  vapor  concen¬ 
tration  downwind.  The  analysis  uses  the  results  of  a  vapor  dispersion 
model  shown  in  Section  6.5  and  the  plume  rise  theory  indicated  earlier. 
The  results  obtained  are  then  compared  with  the  experimentally  obtained 
NH^  vapor  concentration.  For  the  calculations,  the  following  data 
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m  ! 


are  assumed: 


Quantity  of  spill 

Partition  Ratio 

Distance  of  raft  downwind 
at  which  the  vapor 
concentration  is  to  be 
calculated 


4.0  gallons 
0.6 


55  feet 

?.5  mph  =>  11  feet/second 

q  *  2.25  llters/’minute 
1.324  x  10“3  feet^/second 

0% 

mv  **  40  x  5.69  x  0.4  *  91.04  pounds 


Wind  velocity 

Impinger  suction 
rate 

Fraction  of  aerosol 
in  vapor 

Hence  mass  of  vapor  liberated 

From  Figure  5-llb,  the  height 

(zq)  of  the  center  line  of  the 

vapor  plume  at  55  feet  from  the 

spill  center  =>  15  feet 

Gifford  values  of  dispersion  parameters  from  Figure  6-17  for 

atmosphere  B  are:  a y  ®  10  feet,  o z  «  7.5  feet. 

Assuming  that  the  wind  is  steady,  and  that  the  maximum  concentra¬ 
tion  at  any  downwind  location  occurs  in  the  same  spatial  position  during 
the  passage  of  the  cloud,  we  can  derive  the  following  relationship 
between  the  moles  of  vapor  collected  by  an  impinger  at  height  z^, 
directly  downwind. 

From  Section  6.5,  the  concentration  of  vapor  at  any  location  (y,z) 
at  a  distance  x  downwind  can  be  adequately  represented  by 


(y-yj  <2-0 


m  (t-x/U) 

c(x,y,x,t)  «  - 7-t - -t-t—-  e 

v  J  ’  it  a  (x)  o  (x)  U 

v  z 


2  a., 


(z+zo) 


2a  1 
e  z 


+  e 


2o 


and  m. 


lv*ulo  c(x»y0»zt»t>  dt 


(5-1) 


(5-2) 


where 


mv  =  moles  of  vapor  collected  by  an  impinger 

m^  «  rate  cf  evaporation  of  vapor  (a  function  of  time) 


124 


c  »  vapor  concentration  (in  density  units) 


yo.<20  "  position  of  the  maximum  concentration 
t  *  time 

q  «  impinger  suction  rate 

U  »  molecular  weight  »  17  gras/mole  *  0.375  pounds/mole. 

Assuming  that  the  peak  concentration  occurs  at  the  position  of  the  center 
line  of  the  cloud,  we  can  show  that 


M 


m  q 
vH 


v  t To  a  Up  e 

y  2 


(Vzo  >' 

~2a~T~ 

z 


Substituting  the  values  with  »  21  feet 


we  get  M  « 


91.04  x  1.324  x  10~3 
V  tt  x.10  x  7.5  x  11  x  0.0375 


(21-15) 

2  x  7. 5" 


"  90  x  10  **  moles  (5-3) 


This  value  is  compared  with  the  experimentally  observed  value  (see 

-5 

Figure  5-Sa)  of  60  x  10  moles.  The  cheoretical  value  is  about  50% 
more  than  the  experimental  value.  Wind  direction  is  never  steady 
in  the  atmosphere.  This  might  be  one  reason  for  lower  mass  collected 
by  the  impinger.  Also,  because  the  raft  is  floating,  the  impinger 
direction  may  not  have  been  parallel  to  the  wind  direction,  leading 
to  lower  mass  of  NK^  vapor  collected.  T.t  is  also  possible  that  the 
impinger  itself  was  sucking  air  at  a  rate  less  than  2.25  liters/minute 
(because  of  the  long  tubing  in  the  vacuum  line) ,  Other  possibilities 
are  discussed  in  Section  5.5. 

5.4.5  Spill  Radius 

The  maximum  spill  radius  measured  as  a  function  of  the  spill  quantity 
is  shown  plotted  in  Figure  5-7.  Note  that  a  straight  line,  of  elope 
0.8  on  the  logarithmic  plot,  correlates  most  of  the  data.  However, 
if  the  same  power  law  with  volume  is  assumed  to  hold,  then  the  radius 
of  a  3,000  ton  spill  would  be  about  36,000  feet  -  a  highly  Improbable 
number.  Hence  care  should  be  exercised  when  the  small-scale  test 
results  are  extrapolated  to  large-scale  spills.  A  theoretical  model 
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has  bean  worked  out  and  reported  elsewhere  to  predict  the  maximum 

radius  of  spill  for  the  spread  on  water  of  instantaneously  released 
liquid  natural  gas.  The  results  of  the  analysis  are  as  follows: 


R 


max 


(5-4) 


t 


e 


0.674 


(5-5) 


where 

V  *»  volume  of  spill 

y  «  linear  liquid  regression  rate  (assumed  to  be  a  constant) 

/  plnh3 

6  *  affective  gravity  *  g  (  i  - - — • 

\  p water 

t  °  time  to  complete  the  evaporation 

The  3/8  ths  power  law  variation  of  radius  with  spill  size  is  also 
plotted  in  Figure  5-7.  Note  that  the  line  seems  to  correlate  the 
large-scale  field  observations.  Based  on  equation  (5-4)  and  the 
experimental  results  of  the  50-gallon  tests,  we  find  the  "liquid 
regression"  rate  (y)  to  be  2.8  inches/minute.  Using  this  value,  the 
time  to  complete  evaporation,  calculated  from  equation  (5-5)  is  9.5 
seconds.  This  time  is  compared  to  the  experimental  va3ues  of  between, 
five  and  11  seconds.  The  agreement  seems  to  be  very  good.  Using  ..he 
above  formulae  for  prediction,  a  3,000-ton  instaxitaneous  spill  of 
LNKg  on  water  would  spread  to  a  maximum  radius  of  about  475 
feet  and  evaporation  would  be  complete  in  about  120  seconds. 

5.5  DISCUSSION 

The  partition  ratios  have  been  calculated  from  the  impinger  data 
using  the  "mass  balance"  method.  The  validity  of  this  method  for 
obtaining  partition  ratios  has  already  been  tested  in  the  case  of 
intermediate-scale  tests.  The  method  gives  reasonably  accurate  values, 
provided  the  impinger  data  are  good  and  wind  velocity  is  known  with 
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good  precision.  (See  Table  4-7.)  In  the  large-scale  testa,  the 
partition  ratios  calculated  vary  from  a  low  of  0.497  to  a  high  of  0.98. 
The  latter  high  values  have  been  obtained  in  tests  where  either  the 
impinger  data  are  peculiar  or  visual  observations  indicated  the  pas3ege 
of  the  vapor  cloud  above  the  rafts.  The  logical  conclusion  is,  there¬ 
fore,  that  in  such  cases  a  major  mass  of  vapor  did  not  pass  through 
the  impinges,  resulting  in  high  values  of  the  partition  ratios.  How¬ 
ever,  in  the  tests  where  the  wind  was  good  and  the  bulk  of  vapor  cloud 
passed  thr.ough  the  rafts,  the  partition  ratios  are  reasonable,  giving 
credence  to  the  method  of  calculating  the  vapor  mass  from  just  the 
time-integrated,  spot-sampled  vapor  concentration  data.  Recall  that 
the  major  assumption  in  this  method  of  calculation  vapor  mass  is  that 
the  concentration  distribution  profiles  everywhere  in  apace,  at  a  given 
location,  must  be  similar.  This  assumption  should  be  reasonably  valid 
for  steady  wind  conditions.  Even  in  cases  where  the  wind  meanders 
very  much,  the  above  assumption  will  be  valid,  provided  the  time 
of  passage  of  a  cloud  at  any  given  location  is  small  compared  to  the 
period  of  meandering  of  the  wind. 

The  Gaussian  nature  in  the  distribution  of  the  vapor  concentrations 
(impinger  data)  ia  clearly  indicated  in  the  sample  data  shown  in 
Figures  5-8  through  5-10.  The  horizontal  distributions  conform  to 
the  Gaussian  better  than  do  the  vertical  distributions.  In  the  low- 
wind  tests,  very  peculiar  vertical  distributions  have  been  noticed.  As 
explained  earlier,  the  main  reason  for  this  is  the  passage  of  the  cloud 
wall  above  the  rafts.  Another  reason  may  be  the  movement  of  the 
rafts  (which  are  floating)  during  the  cloud  passage  resulting  in  the 
impinger  mouths  being  not  in  line  with  the  wind  direction.  While 
the  possibility  always  existed  that  the  wind  direction  was  not  always 
parallel  to  the  Impinger  Intake  tube,  the  movement  of  the  rafts 
within  the  duration  of  the  passage  of  vapor  cloud  is  so  small  as  not 
to  materially  affect  the  impinger  data. 

The  plume  rise  data  has  been  shown  in  two  graphs.  One  of  than 
has  been  obtained  from  the  peak  concentration  locations  in  the  z  direc¬ 
tion  distribution  of  impinger  data.  In  view  of  the  comments  in  the 
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previous  paragraph,  it  is  very  likely  that  these  peak  location  data 
are  not  very  accurate.  In  fact,  the  2  direction  distributions  seam 
to  be  reasonable  only  for  the  first  raft.  Figure  5-lla  is  based  on 
this  data.  This  plume  rise  data  has  to  be  compared  with  that  obtained 
from  observations  of  the  visible  cloud  behavior  (plotted  in  Figure  5-llb) . 
The  two  figures  give  the  impression  that  the  maximum  vapor  concentra¬ 
tion  occurs  below  the  center  of  the  visible  cloud.  However,  before 
accepting  this  as  always  being  the  case,  note  that  the  data  are 
really  insufficient  to  draw  any  conclusions.  In  fact,  the  impinger 
data  are  time- integrated  values  of  concentration  and  therefore  may 
not  reflect  the  true  position  of  the  peak  concentration  during  the 
cloud  passage.  Also,  there  are  other  uncertainties  in  the  z  direction 
Impinger  data,  as  discussed  earlier. 

The  plume  theory  predictions  agree  with  the  observation  data 
(Figure  5-llb)  very  well.  In  addition,  the  experimental  data  and  the 
plume  theory  have  the  same  power-law  relationship  with  wind  velocity. 

This  agreement  with  "plume"  theory,  instead  of  with  "puff"  theory  is 
intriguing,  because  in  the  experiments  the  vapor  was  observed  to 
disperse  in  the  form  of  a  spherical  cloud  moving  downwind  rather 
than  as  a  plume.  However,  It  is  likely  that  in  spite  of  the  fact 
that  the  vapor  seems  to  be  liberated  in  the  form  of  a  puff,  the  time 
duration  of  10  seconds  (for  50  gallons)  is  long  enough  for  the 
vapor  to  move  as  a  plume.  It  Is  quite  uncertain  how  the  yapor  liberated 
by  a  really  massive  spill  (occurring  in  an  extremely  short  period)  would 
behave.  However,  application  of  puff  theory  to  the  prediction  of  the 
movement  of  cloud  originating  from  an  atomic  blast,  has  proved  to 
be  very  successful  . 

Because  of  the  location  of  the  source  at  ground  level,  the  surface 
boundary  layer  wind  structure  is  important  in  the  dispersion  of  the  vapor. 
The  wind  velocity  in  this  surface  boundary  layer  varies  from  zero  cn 
the  ground  (small  but  finite  value  on  the  water  surface  because  of 
water  movement)  to  an  essentially  constant  value  at  the  top  of  the 
boundary  layer  (of  height  50  to  100  meters) .  A  plume  rise  model  incor¬ 
porating  the  variable  character  of  the  wind  has  been  worked  out,  but 
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because  of  lack  of  information  on  the  vertical  velocity  distribution 
of  the  wind,  no  comparison  cpuld  be  made  with  experimental  vapor  cloud 
rise  observations.  However,  the  simple  pltsne  theory  baa  predicted  the 
experimental  observations  quite  well.  In  obtaining  the  "predictions," 
the  mean  pool  radiua  was  assumed  to  be  0.707  of  the  maximum  pool  radius 
end  the  partition  ratio  was  aasumed  to  be  60%. 

To  explain  the  results  shown  in  Figure  5-lla  by  the  plume  theory, 
one  has  to  assume  a  large  fraction  of  aerosol,  up  to  about  20% 
by  weight,  in  the  vapor..  Even  though  it  is  likely  that  the  vapor 
liberated  has  considerable  amounts  of  aerosol,  20%  by  weight  seems  to 
be  rather  a  large  fraction.  (Unfortunately  there  are  no  direct 
measurements  of  the  amount  of  aerosol  in  the  vapor  cloud*)  It  is 
therefore  concluded  that  visual  observation  data  give  the  best  indica¬ 
tion  of  the  location  of  the  peak  vapor  concentration. 

The  dispersion  parameter  values  (o's)  obtained  from  the  experimental 
data  and  those  obtained  from  Gifford  curves  are  plotted  in  Figure  5-12. 

Note  that  Gifford  curves  were  developed  from  a  large  number  of 
tracer  experiments  with  continuous  sources.  Also,  Gifford  curves 
are  given  only  for  distances  greater  than  100  meters.  For  comparison 
with  the  experimental  data,  these  curves  have  been  extrapolated  to 
the  raft  distances,  which  are  much  closer  than  100  meters.  Also, 
the  atmosphere  type  for  the  days  of  the  experiments  were  inferred 
from  the  sunshine  data,  the  type  of  vind,  kinds  of  clouds,  etc.  Since- 
these  are  all  subjective  information,  it  is  possible  to  make  some 
errors  in  judgment  and  misconstrue  the  class  of  the  atmosphere. 
Considering  all  these  problems,  the  agreement  between  experimental 
and  Gifford-dispereion  parameter  values  are  quite  excellent.  There¬ 
fore,  for  predicting  the  NHg  vapor  concentrations,  Gifford  values  can  be 
used. 

A  possibility  exists  that  all  the  atmospheric  classes  indicated 
in  Table  5-5  are  inecrreou.  When  s  cold  vapor  such  as  saturated 
HKg  vapor  disperses  ever  a  water  surface,  it  may  create  its  own  local 
inversion  condition  in  she  atmosphere.  Since  dispersion  of  the  cloud 
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of  vapor  depends  to  a  considerable  extent  on  the  local  atmospheric 
condition,  and  since  inversion  represents  a  stable  situation,  it  is 
feared  that  the  vapor  cloud  may  not  disperse  as  quickly  as  expected. 
However,  experimental  observations  do  not  seem  to  justify  these  fears. 

The  dispersion  example  worked  out  in  Section  5. 5. A  shows  agreement 
between  the  predicted  concentration  and  the  e.v' rimental  value. 

Note  that  a  simple  continuous  point  source  forsi'.  la  was  used  in  the 
concentration  prediction.  Thus,  the  answers  are  expected  to  be  valid 
for  distances  greater  than  about  two  times  the  maximum  pool  diameter 
from  the  pool  center.  The  example  shown  in  Section  5.4.4  indicates 
an  error  of  about  50%  in  prediction.  Though  this  error  seems  to  be 
high,  it  is  not  unreasonable  considering  the  fact  that  we  are  comparing 
two  numbers  which  are  the  time  integrated  values  of  the  concentration. 
Experimentally,  there  are  many  variables,  the  important  ones  of  which 
are  the  wind  velocity  and  direction.  The  assumptions  made  in  the 
dispersion  model  as  to  the  steadiness  of  the  wind  may  not  be  that 
accurate.  However,  that  the  predicted  and  experimental  values  agree 
that  closely  is  in  itself  a  proof  of  the  adequacy  of  the  dispersion 
model  to  predict  reasonably  accurate  concentrations. 


6.  THEORETICAL  ANALYSES 

6.1  INTRODUCTION 

In  the  previous  chapters,  several  experimental  results  have  been 
presented,  analyzed,  and  discussed.  In  this  chapter,  several  theoretical 
analyses  are  given  vhich  are  relevant  to  the  prediction  of  hazards 
presented  by  a  spill  of  liquid  asnonia  on  water.  These  analytical 
models  provide  a  foundation  for  suitably  extrapolating  small-scale 
experimental  data  to  bigger  scale  spills  that  might  occur  due  to  an 
accident . 

Section  6.2  is  a  discussion  of  the  vapor  flash  model,  which 
provides  a  bawls  for  calculating  the  amount  of  vapor  produced  when 
liquid  ammonia  from  a  pressurized  tank  leaks  out  tc  the  atmosphere. 

Section  6.3  presents  a  thermodynamic  mixing  aodei  for  the  mixing 
cf  ammonia  and  water,  where  it  can  be  seen  that  q;iite  different 
answers  are  obtained  depending  on  whether  ammonia  mixes  with  water 
or  vice  versa.  In  an  actual  spill,  both  types  of  mixing  may  take 
place.  The  results  of  the  model  are  the  partition  ratio  and  the 
temperature  of  the  mixed  solution. 

Section  6.4  contains  a  detailed  heat  transfer  analysis,  made  to 
obtain  tho  pressure-time  relation  for  the  tank  pressure  of  a  liquid 
ammonia  tank  submerged  in  water.  Such  a  situation  may  occur  following 
the  collision  of  an  ammonia  barge  and  its  subsequent  sinking 

Sections  6.5  and  6.6  present  a  compendia  of  models  for  the 
dispersion  of  ammonia  vapor  in  air  and  ammonium  hydroxide  in  water ;  i.e., 
the  various  formulae  used  in  different  situations,  rather  than  detailed 
derivations.  Most  of  the  formulae  have  proved  accurate  for  predicting 
various  other  gas  dispersions  and  water  pollutant  concentrations. 

Section  6.7  presents  the  predictions  of  the  plume  ana  vapor  puff 
rise  theories,  which  are  based  on  the  analysis  of  the  motion  of  a  buoyant 
cloud  o£  gas. 


An  underwater  release  reaction  theory  in  presented  in  Section 
6.8.  This  model  is  based  on  an  anology  between  the  behavior  of  turbulent 
flames  and  LNK^  Jet* 

6.2  FLASH  MODEL 

Ammonia  is  transported  as  a  liquid,  either  under  pressure  at  ambient 
temperature,  or  at  atmospheric  pressure  as  a  cryogenic  liquid.  When  the 
pressurized  liquid  is  released,  a  fraction  of  the  liquid  vaporises 
very  rapidly  and  escapes  into  the  atmosphere.  The  remaining  liquid 
cools  to  below  its  atmospheric  pressure  liquefaction  temperature. 

The  present  analysis  is  to  obtain  a  method  of  calculating  this 

fraction  of  vapor  produced.  In  deriving  the  following  equation,  we 

have  assumed  that  the  process  is  adiabatic.  It  can  be  shown  that 

.  .  * 

the  process  is  also  isenthalpic.  Thus,  we  can  equate  the  Initial 
and  final  enthalpies  of  a  unit  mass  of  LNH^.  Assuming  no  supercooling 
of  the  liquid,  we  can  write  the  enthalpy  equation  as 


»!  -  fh° +  «-» 


where 


(6-1) 


f  =  fraction  of  released  LKK^  that  vaporizes 
h  »  specific  enthalpy  (Btu/lbm) 

A°  «  heat  of  vaporization  at  one  atmosphere  saturated  condition 

superscript  o  «  fiual  (outside  the  tank)  value  (one  atmosphere 
pressure  saturated  condition) 

superscript  i  =  initial  (inside  the  tank)  value 

subscript  £  *  liquid 

subscript  v  =  vapor 


Hence, 
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Illustrative  examples  have  been  worked  out  for  three  cases  of 
tank  conditions,  and  the  results  are  given  in  Table  6-1.  Note  that 
the  fraction  of  vapor  formed  increases  as  tank  pressure  increases. 
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Table  6*1 


Vapor  Flash  Calculations  for  Three  Tank  Conditions 


Ambient 

Temperature 

(°f) 

Tank 

Pressure 

(pai't 

(¥> 

Enthalpy  of 
Tank  Liquid 
C&tu/lba) 

heat  of 
Vaporijsstitnj 

X° 

at  atmospheric 
pressure 
(Bfiu/lbm) 

1 

Vapor  fraction  (f )  j 

formed  at  ‘ 

1  ata.  pressure  . 

0 

30.42 

42.9 

589.3 

1 

1 

0.051 

i 

40 

73.32 

86.8 

589.3 

0.126  | 

100 

211.90 

155.2 

589.3 

0.242  ; 
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A  computer  program  has  also  been  written  to  calculate  the  vapor 
fraction,  given  the  ambient  temperature.  The  program  utilizes  an 
"ammonia  properties  subroutine"  to  evaluate  property  values. 

6.3  THERMODYNAMIC  MODELS 

6.3.1  Model  Types 

When  LNE^  is  spilled  on  water,  it  reacts  with  water  in  a  relatively 
small  area  and  forms  ammonium  hydroxide  (NH^OH)  solution  and  gaseous 
ammonia  vapor  (GNH^).  The  reaction  is  quite  rapid  and,  for  all 
practical  purposes,  heat  loss  from  the  reaction  can  be  neglected. 

Thus,  the  mixing  can  be  considered  to  be  adiabatic.  The  mixing  process 
can  be  modeled,  using  an  equilibrium  thermodynamic  approach. 

However,  depending  on  the  details  of  the  mixing  process,  different 
answers  are  obtained.  This  behavior  is  not  unlike  the  differences  found 
between  adding  drops  of  water  to  concentrated  sulfuric  acid  and  adding 
drop3  of  sulfuric  acid  to  water.  In  this  section,  the  three  different 
thermodynamic  models  that  explain  the  various  mixing  reactions  are  derived 
and  discussed.  The  first  case  is  most  representative  of  surface  spills  of 
LNH-j  on  water.  It  can  explain  most  of  the  experimentally  observed  phenomenon 
with  the  exception  of  certain  temperature  measurements  observed  in  the 
larger  size  spill  tests. 

Case  1.  A  "packet"  of  LNB^  mixes  with  increasing  amounts  of  water, 
thereby  undergoing  dilution,  and  liberating  ammonia  vapor  (GNH^). 

The  ammonium  hydroxide  (NH^OH)  formed  mixes  again  with  water,  under¬ 
going  further  dilution.  The  "packet"  is  assumed  to  be  at  uniform 
concentration  throughout  its  volume,  and  the  vapor  produced  by  dilution 
is  assumed  to  be  removed  out  of  the  reaction  zone. 

A  model  derived  on  these  assumptions  is  presented,  and  an  expression 
for  the  maximum  vapor  liberated  by  a  unit  mass  of  LNH7  undergoing 
infinite  dilution  is  obtained.  The  results  of  the  model  are  compared 
with  experimental  data  and  discussed. 

Case  2.  A  given  mass  of  LNH^  reacts  adiabatically  with  a  known  mass 
of  water.  The  final  condition  obtained  is  an  equilibrium  state,  with 
the  final  NH^OH  solution  in  equilibrium  with  its  vapor  at  saturation. 
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The  derived  model  presents  expressions  for  calculating  both  the 
total  mass  of  vapor  produced,  and  the  temperature  of  the  final  mixture 

(NH^OH)  for  various  final  concentrations.  The  results  of  this  model 
are  compared  with  some  experimental  results. 

Csss  3.  The  reaction  is  visualized  to  be  between  e  given 
mass  of  water  and  increasing  quantities  of  Pfocess  quite  the 

opposite  of  Case  1) .  A  model  based  on  this  premise  has  been  worked 
out  and  its  results  are  given. 

6.3.2  Continuous  Water  Mixing  (Case  1) 

Figure  6- la  shows  a  schematic  representation  of  Case  lt  in  which 
water  is  continuously  mixed  with  NtLOH  solution  (starting  with  LNH^; 
i.e.,  100%  UH^OH  solution)  under  adiabatic  conditions  and  at  one 
atmosphere  pressure.  As  vapor  is  liberated,  it  is  removed  from  the 
liquid.  The  solution  remains  saturated  at  its  local  condition  at 
all  times. 

The  following  equations  may  be  written  for  the  conditions  of 
the  solution  before  and  after  mixing  with  the  quantity  "dM"  of  water. 

Sefp.rrin3  to  the  control  volume  (shown  in  Figure  6-lb)  we  have, 

m  =  m(c)  **  mesa  of  ammonia  (NH^)  in  the  solution  (NH^OH)  at  any 
concentration 

M  m  mass  of  water  in  the  solution 

c  ■  concentration  of  the  NK^CH  solution  »  If  (6-3) 

For  the  addition  of  "dM"  mass  of  water 

Initial  enthalpy  of  the  system  -  (m  +  M)I(c)  +  dM  1^  (6-4) 

Final  enthalpy  of  the  system  after  mixing  »  |  (a  +  dm)  +  (M  4-  dM)  |  I(c  +  dc) 

(6-5) 


Where  I  refers  to  the  specific  enthalpy  of  mixture  at  any  concentration 
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Figure  6— la  Thermodynamic  Mods!  for  Mixing  of  LNH^  and  Water 
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Figure  6— 1b  Control  Volume  for  thu  Mating  of  LNH3  snd  Water 
FIGURE  S-1  SCHEMATIC  DIAGRAMS  ILLUSTRATING  THE  THERMODYNAMIC  MIXING  MODEL  1 


136 


Since  the  mixing  is  isenthalpic,  we  can  equate  equations  (6-4)  and 
(6-5).  Expanding  I  (c  +  dc)  in  Taylor  series,  simplifying,  and  rearrang¬ 
ing  using  equation  (6-3)  gives 


dm 

dM 


m 

c 


j  i1  - c  §  -  V  l 


or  alternacively 


(6—6) 


.  H  I  [I  '  c  H  “  1w]  |  (6-7) 

dc  C  (  [(1”C)  (I-Iw)  -  cX]  ) 

where 

X  =  X(c)  =  Heat  of  vaporization  of  NH^OH  at  any  concentration  c. 

The  terms  within  brackets  on  the  righthand  side  of  equations  (6-6)  and 
(6-7)  are  functions  of  only  the  concentration  (and  hence  are  obtainable 
from  the  NH^-water  system  property  data)  and  the  condition 
of  the  water.  The  two  equations  can  therefore  be  solved  (for  a  given 
initial  value  of  m,  the  mass  of  LNHg)  to  obtain  m  as  a  function  of 
the  mass  of  water  added  (M) ,  or  the  NH^OH  concentration  (c) . 

Without  loss  of  generality,  and  for  the  sake  of  a  simple 
solution,  we  assume  m  =  1  at  c  *  1  for  the  Initial  condition: 
i.e.,  a  unit  mass  of  LNH^  to  begin  with.  Hence  the  total  mass 
of  vapor  produced,  by  the  time  the  LNHg  is  diluted  to  a  concentration 
c,  is 

in  ■*  1  -  m(c)  (6-8) 

The  physical  meaning  of  the  term  I  -  c  ~  (called  "partial  euthalpy 
of  water")  in  the  righthand  side  of  equation  (6-7)  is  shown  in  Figure 
6-2;  the  enthalpy  of  the  saturated  NH^OH  solution  at  atmospheric  pressure 
is  represented  by  point  E.  Point  D  is  the  ii  tersecticn  of  the  ordinate 
axis  and  the  tangent  (DE)  at  E  to  the  atmospheric  pressure  enthalpy 
curve,  while  BD  represents  the  value  of  (I  -  c  ^)  at  the  concentration 
c.  Point  A  on  the  ordinate  axis  represents  the  condition  of  water. 
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Distance  AD  represents  the  magnitude  of  the  term  (I  -  c  -  1^) . 

As  the  3et^  diluted,  the  intersection  point  of  the  tangent  line 
with  the  ordinate  moves  up  the  ordinate  until  the  point  D  coincides 
with  point  A  at  a  "critical  concentration"  c*. 

For  c  >  c*,  >0,  indicating  thereby  that  a  decrease  in 

concentration  (dilution)  results  in  the  liberation  of  vaport. 

At  c  »  c*,  -j-  «  0. 

*  dc 

rim 

For  c  <  c*,  -3—  <  0;  i.e.,  for  a  further  dilution  below  s*,  the 

uC 

mass  of  ammonia  in  the  solution  should  increase. 

The  last  case  is  not  allowable,  because  for  that  to  happen  some 
ammonia  vapor  should  dissolve  in  the  solution.  However,  the  basic 
assumption  in  the  model  is  that  any  vapor  produced  is  removed  immediately 
and  is  not  allowed  to  redissolve.  Because  of  this  stipulation,  dilu¬ 
tion  of  the  solution  below  the  concentration  c*  results  in  unsaturated 
HH^OH  solution,  without,  however,  changing  the  mass  of  ammonia  that  is 
in  Che  solution  at  concentration  c*.  The  solution  temperature  is  no 
longer  given  by  the  saturation  temperature  corresponding  to  the  concen¬ 
tration. 

Hence  the  maximum  amount  of  vapor  produced  by  the  above  mixing 
process  is  given  by 

Total  vapor  produced  *  m^Cc*) 

mass  of  ammonia  ir  solution  after) 

Therefore  partition  ratio  p  *  ammonia  spillid"  (6-9) 

=  1—nVy.Cc*)  =  m(c*)  / 

The  solution  of  equation  (6-7)  is  obtained  numerically  by  curve  fitting 
for  the  properties  that  occur  on  the  righthand  side.  Figure  6-3  shows 


fNote  that  vapor  is  liberated  when  "dm"  is  negative. 
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FIGURE  6-3  PLOT  OF  PARTIAL  ENTHALPY  OF  WATER  AS  A  FUNCTION  OF  CONCENTRATION 
FOR  NH,  -WATER  SYSTEM 


dX 

the  plot  of  the  partial  water  enthalpy  I  -  plotted  against  the 
concentration.  The  critical  concentration  is  determined  by  equating 
the  enthalpy  defect  and  the  initial  water  enthalpy  (for  70°P  water, 
c*  a  0.46).  The  solution  obtained  in  terms  of  total  mass  of  vapor 
liberated  i a  plotted  in  Figure  6-4  for  two  values  of  water  tempera¬ 
ture.  The  partition  ratios  given  in  Table  6-2  were  obtained  for 
two  different  water  temperatures. 

for  any  concentration  c  below  c*,  the  temperature  of  the  NH^OH 
solution  is  not  equal  to  the  saturation  temperature  corresponding  to 
the  concentration  c.  The  following  analysis  indicates  the  method 
used  for  calculating  the  temperature  solution  below  c*. 

As  before,  the  enthalpies  of  the  system  before  and  after  the 
addition  of  a  ndM"  mass  of  water  are  equated.  Hence  we  get 

(m*  +  M)  [1(c)  -  c  (Tsat  (c)  -  T)]  +  dM  I 
p  w 


Initial  Enthalpy  of  System 

“  (m*  +  M  +  dM)  [I(c  +  dc)  -  c^  TSat(c  +  dc)  -  (T  +  dT)  ] 


Final  Enthalpy  of  the  mixture 


(6-10) 


where 


m*  =  mass  of  ammonia  in  the  solution  whan  the  concentration  is  c* 

(This  value  remains  s  constant  during  the  subsequent  dilution.) 

«  specific  heat  of  NH,0H  solution  (assumed  independent  of 
temperature  and  concentration) 

(c)  *  saturation  temperature  at  concentration 

“  temperature  of  the  solution  at  concentration  c 

m* 


P 

sat 


Also  c*  * 


and  dM  =* 


m*  +  M 


a* 

2 


ac 


(6-11) 
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Water  Temperature 


I 

i 


I 


I 


Table  6-2 


Theoretically  Predicted  Partition  Ratios 
and  Maximum  NK^CH  Solution  Temperature 


Jiaxiaua 


Mixing 

Hater 

Temperature 

°T 

Critical 

Concentration 

* 

Partition 

Ratio 

Temperature 
of  NTI.0K 
Solution 
°F 

Concentration 

at  which  Maximum 

Temperature  Occurs 
* 
h 

32 

50,5 

0.775 

43 

~30 

70 

46,0 

0.735 

70 

0.0 
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Expanding  the  terms  on  the  rlghthand  side  of  equation  (6-10) 
in  a  Taylor  series,  then  rearranging  and  simplifying  using  equation 
(6-11) ,  it  can  be  shovn  that  equation  (6-10)  reduces  to 


d  (I  -  c  (TSat  -  T)  -  I  ) 

_  P  W 

(1  -  c  (TSflt  -  T)  -  I  ) 
p  w 


(6-12) 


Integrating  the  above  equation  and  substituting  the  initial  condition 
that, 

T  «  Tsat(c*)  at  c  -  c* 

we  obtain  the  following  solution  for  the  temperature 

T(c)  -  Tsat(c)  +A  [(Iw  -  Kc))  +  f*  (I  (c*)  -  Iw)3  (6-13) 

This  result  is  plotted  in  Figure  6-4  for  two  water  temperatures  - 
32°F  and  7QaF  -  along  with  the  saturation  temperature  as  a  function 
of  the  concentration.  The  result  indicates  that  when  cold  water  is 
used,  the  maximum  temperature  attained  by  the  solution  is  higher  than 
the  water  temperature,  and  that  for  moderate  and  high  water  temperatures 
no  temperature  rise  is  observed, 

6 , 3 . 3  Complete  Adiabatic  Mixing  (Case  2) 

In  this  model  we  again  discuss  the  adiabatic  mixing  of  LNE^  and 
water,  but  the  feature  that  distinguishes  this  model  from  that  of 
Case  1  is  its  complete  equilibrium  between  vapor  (if  produced)  and 
the  final  NE^OH  solution.  The  results  are  independent  of  the  ‘'path” 
in  which  the  final  situation  is  attained. 

The  basic  aspects  of  the  model  may  be  seen  in  the  idealized 
experiment  shown  in  Figure  6-5. 

A  certain  mass  of  LHH^  at  saturation  (atmospheric  pressure)  and 
a  unit  mass  of  water  at  given  temperature  are  kept  inside  an  adiabatic 
cylinder,  separated  by  an  Ideal  membrane  as  shown  in  Figure  6-5a  with 
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Figure  6~ 5a  Ffgura  6— Ee 
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the  cylinder  pressure  maintained  at  atmospheric  value.  At  some  instant, 
the  membrane  is  punctured  and  the  LKH^  and  water  are  allowed  to  mix, 
react,  and  come  to  equilibrium  to  form  an  NH^OH  solution  of  a  desired 
concentration.  It  is  desired  to  find  what  initial  mass  of  LNH7  is 
needed  to  attain  a  specified  concentration  end  to  determine  the  mass 
of  ammonia  in  the  solution  at  the  final  condition. 

Since  the  process  is  adiabatic  and  a ^  constant  pressure,  it  can 
be  shown  that  it  is  also  isenthalpic.  If, 

x  =  mass  of  INK,  to  be  taken  initially  to  obtain  a  final  concentration  c 
s  v 

mvj,  m^  »  mass  of  water  in  the  liquid  phase  end  vapor  phase,  respectively, 
in  the  equilibrium  mixture  at  c 

_  £  V 

and  m  ,  a  «  mans  of  ammonia  in  the  liquid  and  vapor  phases,  respectively, 
then  we  csa  writs  the  following  enthalpy  balance  equations. 

If  no  vapor  is  formed 


a 


(6-14) 


ifc  vale's  case,  the  NS^OH  solution  formed  is  not  saturated;  and  if 
vapor  is  formed 

°a.je3  +  +  v i(c>  +  +  #  ww 

To  evaluate  the  value  of  x  from  equation  (6-15) ,  the  following 
additional  equations  are  used: 


(6-15) 


a  e  {-~ — *)  rj 
a  i  -  c  v 


1  ■=  a  -<•  m 

V  "i 


m  =-  a  +  a 
z  a  a 


(6-16) 


v  - ,  -  v 
ana  e.. 
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where  R(c)  is  a  factor  obtained  fro:*  property  data.  In  general, 

y 

R  «  1,  so  that  for  all  practical  purposes  one  can  assume  a  «  0 

i  '  w 

(hence  m  *  1) 

Substituting  the  set  of  equations  (6-16)  in  (6-15),  and  knowing 

all  the  enthalpy  values  from  property  data,  rhe  value  of  can  be 

evaluated.  Partition  ratio  can  then  be  calculated  by 
t 


(6-17) 


When  no  vapor  is  formed  (which  happens  when  small  amounts 
of  ammonia  mix  with  large  quantities  of  water)  the  final  solution 
obtained  is  not  saturated.  In  such  a  case,  the  temperature  of  the 
mixture  can  be  calculated  from 


T  *>  lBat  (c)  4* 


1_ 

c 

p  L 


f  (\  XLNK~  *  IW) 

_ l _  -  1(c) 

[  (1  +  n) 


(fc-13) 


where  is  given  by  equation  (6-14) .  The  temperature  of  the  solution 
which  is  In  equilibrium  with  the  vapors  at  any  concentration  la  equal 
to  the  saturation  temperature  corresponding  to  that  concentration. 

The  results  calculated  for  two  water  temperatures  are  shown  in 
Figures  6-6  and  6-7.  figure  6-6  snows  the  total  mass  of  water 
necessary  to  make  the  neceegary  concentration  solution,  starting 
with  &  unit  mass  of  LNHj,  and  gives  the  total  NH^  vapor  produced  at 
an y  equilibrium  concentrations  Figure  6-7  shows  the  dependence 
of  the  equilibrium  temperature  and  partition  ratio  on  the  final 
concentration.  It  is  noteworthy  that  a  substantial  temperature 
increase  (above  water  temperature)  can  result  if  the  final  concentra¬ 
tion  is  between  15  aed  25%. 

6.3.4  Continuous  fixing  (Case  3) 

One  can  visualise  a  situation  where  increasing  amounts  of  LNH^ 
mix  with  a  given  quantity  of  water.  An  idealized  experiment,  similar 
to  the  one  mentioned  for  Case  1,  can  be  pe.  formed  where  instead  of 
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FIGURE  8-6  MASS  OF  VAFOR  PRODUCED  AND  MASS  OF  WATER  NEEDED  TO  PRODUCE  A  GIVEN  CONCENTRATION 
OF  AMMONIUM  HYDROXIDE  BY  ADIABATIC  MIXING  WITH  UNIT  MASS  OF  LNH~ 


water  mixing  with  NH^OH,  LNH„  mixes  with  NH^OH  solution,  starting 
with  pure  water  (i.e.,  300%  NH^GH  solution).  The  other  conditions 
(adiabatic,  immediate  vapor  removal,  etc.)  remain  the  same  as  those 
given  for  Case  1. 

An  enthalpy  balance  equation  can  be  written,  and  expressions 
can  be  obtained  for  the  mass  of  vapor  produced,  temperature  of  NH^OH 
solution,  etc.  Detailed  derivations  are  omitted,  but  the  results 
are  shown  in  Figure  6-8. 

The  results  of  this  model  show  that  even  though  the 
total  mass  of  vapor  liberated  increases  as  the  total  mass  of  LKH^ 
added  increases,  the  partition  ratio  first  decreases  and  then  increases. 
This  indicates  that  up  to  a  certain  concentration  (about  452  for 
70°F  water)  the  mass  of  vapor  produced  per  unit  mass  of  LNK^  added 
first  increases  and  then  decreases.  However,  the  main  difference 
between  this  model  and  the  model  for  Case  1  is  that  in  the  latter 
case  the  concentration  of  the  mixture  continuously  decreases,  and 
after  a  certain  concentration  no  more  vapor  is  liberated.  However, 
in  the  model  for  Case  3,  there  is  a  continuous  increase  in  the  mass 
of  vapor  liberated,  even  though  the  partition  ratio  itself  reaches 
a  minimum,  then  rises  back  to  unity.  This  is  because  of  the  increasing 
quantity  of  LNK^  that  is  added. 

6.4  HAZARDS  CAUSED  EY  SINKING  AN  LNH-j  5ARCS 

6.4.1  Introduction 

A  major  quantity  of  ammonia  is  carried  on  U.S.  inland  waterways 
in  refrigerated  barges  in  amounts  exceeding  2500  tons.  b.  typical 
2500-ton  barge  built  by  National  Marine  Service  by  Bethlehem  Steel 
is  shown  in  Figure  6-9.  The  barge  has  two,  40-foot-long,  18-foot-diameter 
tanks,  insulated  by  three-inch-thick  urethene  foam  and  an  18-gage 
g.'lvanieed  steel  cover.  The  tanks  are  equipped  with  six- inch-diameter 
pressure  relief  valves,  preset  to  open  at  10  psig.  The  valves  are 
pilot  operated  and  have  no  flutter  when  operating  at  the  set  pressure. 
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FIGURE  6-9  DETAILS  OF  A  REFRIGERATED  LIQUID  AMMONIA  BARGE 


If  such  a  barge  i s  involved  in  an  accident,  and  then  sinks  with 
the  LNH^  tank,  a  serious  problem  exists.  Depending  on  the  amount 
of  tank  insulation  lost  in  the  accident,  a  tremendous  amount  of  heat 
leaks  into  the  LNK^  ^rom  the  water.  The  sudden  heating  of  the  LNH^ 
causes  an  increase  in  tank  pressure  and  venting  of  liquid  qr  gaseous 
ammonia  (depending  on  the  orientation  of  the  barge). 

The  analysis  presented  here  concerns  the  consequences  of  such  a 
disaster  and  whether  the  venting  capacity  normally  provided  is  adequate 
to  prevent  the  tank  over-pressure  and  possible  rupture.  The  consequences 
of  an  explosion  releasing  all  of  the  remaining  LNH^  at  once  underwater 
are  not  yet  very  clear. 

6.4.2  Problem  Description 

To  evaluate  the  potential  hazard  posed  by  a  submerged  tank  filled 
with  LNH^,  we  considered  the  "worst"  of  the  several  possible  cases. 

As  shown'  in  Figure  6-10,  a  sunken  tank  can  assume  a  variety  of  orienta¬ 
tions.  The  particular  orientation  depends  on  the  position  of  the  barge 
after  sinking  (note  that  a  full  tank  of  is  buoyant,  by  itself). 

The  worst  case  is  the  situation  where  the  tank  has  bean  inverted  during 
the  sinking  so  that  the  relief  valves  are  on  the  bottom  of  the  tank 
and  must  pass  liquid  to  relieve  the  tank  pressure.  In  such  a  case, 
the  dissolution  of  liquid  ammonia  in  water  takes  place  around  the  tank, 
heats  the  water  around  the  tank,  and  increases  heat  transfer  rate 
to  the  tank.  A  temperature  rise  of  60°  to  70°F  has  been  observed 
in  a  few  instances  in  the  test  program.  We  also  considered  a  situa¬ 
tion  where  the  urethene  foam  insulation  on  the  tank  was  nullified  by 
detachment  or  permeability  to  water,  and  the  relief  valves  were 
malfunctioning . 

The  parameter  of  interest  in  this  analysis  is  the  tank  pressure 
as  a  function  of  time  subsequent  to  immersion  in  water.  Such  auxiliary 
information  as  ullage  volume,  thickness  of  ice  on  the  tank  outer  wall, 
temperature  of  the  tank  contents,  etc.,  was  also  generated.  Two 
specific  cases  have  been  evaluated: 
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*  Initially,  contents  are  saturated  at  14.7  psia.  After  sinking, 
the  tank  is  Inverted  and  releases  liquid  from  the  relief  valves  at 

10  psig.  The  outside  water  temperature  is  100eF,  and  all  of  the 
insulation  is  inoperative. 

#  A  second  case  supposes  that  the  relief  valves  do  not  open  at 
all  and  the  tank  eventually  ruptures. 

6.4.3  Model  Formation 

6. 4. 3.1  Assumptions 

We  have  stated  previously  that  the  principal  output  from  the  model 
is  a  tank  pressure  profile  versus  time  subsequent  to  immersion  in  100°F 
water.  In  assembling  a  model  of  the  tank's  physical  behavior,  we 
made  the  following  assumptions: 

e  tank  contents  are  always  saturated  at  uniform  pressure; 

«  static  head  of  liquid  is  ignored; 

•  heat  transfer  to  contents  is  uniform  over  the  exposed  area 
of  the  tank  cuter  wall; 

©  ambient  water  is  at  100°F;  ^  . 

©  flow  through  each  of  the  two  six-inch  relief  valves  is  modelled 
as  turbulent  flow  through  a  sharp-edged  orifice  downstream 
from  a  long  contraction  region; 

«  relief  valves  are  fully  open  at  11%  over  set  pressure;  and 

•  thermal  capacity  cf  the  steel  tank  is  ignored  in  comparison  with 
the  thermal  capacity  of  the  contents. 

6. 4. 3. 2  Energy  Equation 

Three  distinct  physical  situations  might  occur  during  Che  transient 
temperature  rise  of  the  tank  contents: 

1.  Relief  valves  are  closed  and  both  vapor  and  liquid  phases 
present  within  the  tank: 


„ . 


2.  Relief  valves  are  cloeed  with  the  tank  full  of  liquid} 

3,  Relief  valves  are  open  and  mass  flows  out  of  the  tank  with 
both  vapor  and  liquid  phases  present  within  the  tank. 

For  each  of  the  above  cases,  the  energy  aquation  is  written  in 
terms  of  the  heat  input  rate  and  the  increase  in  the  internal  energy 
of  the  contents.  These  are  given  below,  using  the  following  nomenclature. 


A 

c 

F 

E 


cv 


sf 


rav 

p 

p 

• 

Q 

t 

T 


■  surface  area  of  the  tank  (ft  ; 

■  specific  heat  of  water  (Btu/lbo-°F) 

*  energy  contained  within  the  control  volume  (Btu) 

**  specific  enthalpy  of  liquid  (Btu/lbm) 

«  specific  enthalpy  of  vapor  (Btu/lbm) 

t 

**  heat  transfer  coefficient,  water  to  ice(Btu/hr  -  ft"1 
**  heat  of  fusion  of  water  (Btu/lbm) 

“  thermal  conductivity  of  ice  (Btu/hr  -  ft  -  °F) 

»  thickness  of  the  ice  layer  (in) 

■  liquid  mass  within  the  tank  (Ibw) 

«*  vapor  mass  within  the  tank  (Ibm) 

=  contents  pressure  (psia) 

Zap 


-  °F) 


ft 


time  derivative  of  pressure  (psia/oin) 


®  rate  of  heat  transfer  to  tank  (Btu/min) 

■  time  (min) 

=»  contents  temperature  (°F) 

■»  temperature  of  ambient  water  (°F) 

=  specific  internal  energy  of  the  liquid  (Btu/lbm) 

=  specific  internal  energy  of  the  vapor  (Btu/lbm) 

3 

=  specific  volume  of  liquid  (ft  /lbm) 
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3 

specific  volume  of  vapor  (ft  /lb®) 


V 

e 

* 

a 


tank  volume  (ft  ; 


volume  fraction  of  vapor  (dimensionless)  0  <  6  <  1 


_3_ 

time  derivative  of  vapor  volume  fraction  (ain  ) 


3, 


»  density  of  ice  (llm/ffc  ) 

Situation  1  Relief  valves  are  closed  and  both  vapor  and  liquid 
phases  present  within  the  tank.  For  this  case*  we  define  a  control 
volume  ao  shown  in  Figure  6-11,  and  write  the  First  Law  of  thermodynamics 
for  the  closed  system. 


ft  (Ec> 


Expanding  this  we  get 


3t  (mv  Uv  +  mf  uf> 


(6-20) 


at 


f  0v  -  %>  - 

V  vf 


(h*-  -  Pv,) 
f  r 


(6-21) 


In  terms  of  the  time  derivatives  of  vaocr  volume  and  oressure, 
•  •  *  ■ 

0  and  P  respectively,  First  Law  equation  becomes 


k  "  ^  hf  "  Pvf 

+  ? 

Vv 

vf 

L  v 

I 

L 

dv 


dP  ..  n  ,  .  g  f  dhv  °'v  , 

— r—  (h  -  Pv  )  -r  —  I - -  p  — : 

,  2  v  v  v  \dP  d? 

* 


-) 


dv. 


^  ®  dp  s'!  /dh,  dv. 

(hc  -  ?v£)  +  (1_~  6U  i  (fc-22) 


'f 


dP 


d\?  V 
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where  wa  have  an  equation  in  terns  of  6,  P,  Q  and  the  thermodynamic 
properties  of  ammonia  at  saturation  conditions.  This  is  the  most 
useful  formulation  in  terms  of  our  interest  in  the  vapor  volume 
fraction  0  and  in  the  pressure  ?  as  functions  of  time. 

We  need  an  additional  equation  to  cake  the  problem  amenable  to 
solution.  Assuming  we  can  later  formulate  the  heat  transfer  an 
additional  relationship  b.-.eween  0  and  P  is  necessary.  Conservation 
of  mass  within  the  control  volume  requires  that 


(sn  +  m,) 
St  v  f 


(6-23) 


i.e.,  the  total  system  mass  is  invariant  in  time.  We  can  rewrite  this 
equation  as 
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We  know  that  --=?  may  not  be  constrained  to  be  zero,  so 


3P  (av  +  raf>  =  0 


(6-25) 


This  expression  may  be  expanded  to  yield 
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or,  in  another  form, 
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Now,  we  can  relate  0  and  P  as  follows: 
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Given  an  expression  for  Q,  we  can  aolvn  for  6  find  P  as  functions 

#  * 

of  time,  given  the  two  first-order  differential  equations  in  B  and  P. 

Situation  2  Relief  valves  closed  with  the  tank  full  of  liquid. 
In  this  situation,  any  head  added  will  only  serve  to  expand  the  liquid 
further  in  the  absence  of  vapor.  The  assumption  of  saturated  contents 
requires  that  there  be  a  finite  vapor  fraction  existing  to  have  two 
phases  in  equilibrium,  but  we  may  simplify  the  results  of  Situation  1 
for  6  «  1  in  the  interests  of  a  simpler  solution.  The  equations  then 
reduce  to 


(hj-PV+H; 


*  2 
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(6-29) 


Strictly  speaking,  the  use  of  property  values  for  the  saturated 
liquid  state  is  not  entirely  valid,  because  we’ll  be  considering  the 
compressed  liquid  region  when  the  tank  becomes  full  of  liquid. 

But  this  formulation  viil  presume  two-phase  equilibrium  in  its  use 
of  property  values  for  the  saturated  liquid,  while  at  the  same  time 
allow  easier  solving  for  the  case  when  B  «  1. 

Situation  3  Belief  valves  are  open  and  liquid  flows  out  of ‘the 
tank.  Here,  we  must  consider  an  open  system  (mass  crossing  the  system 
boundary)  enalvsis  of  the  tank.  The  First  La n  formulation  is  now 
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The  left  side  of  the  equation,  (E^) ,  becomes  after  sons  manipulation. 
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The  right  side  of  the  equation,  (Q  -  mQut  h^)  nmy  be  expanded 
in  the  form 
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If  the  right  and  left  sides  are  combined  in  expanded  fora,  the  energy 
equation  yields 
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6. 4. 3. 4  Heat  Transfer  from  Water 


When  a  tank  containing  cold  <-28°F)  LHH^  Is  immersed  in  water, 
ice  begins  to  form  on  the  outer  walls  of  the  bars  portions  of  the  tank. 

Figure  6-12  illustrates  the  ice  layer  and  the  temperature  distri¬ 
bution  of  this  situation.  Because  of  the  high  conductivity  of  the 
metal  wall,  we  assume  negligible  temperature  gradient  in  the  wall. 

Heat  is  liberated  due  to  water  freezing  and  the  subcocling  of  both 
water  and  ice  and  the  heat  is  transferred  to  the  LNH^.  Bv  ignoring 
the  heat  capacity  of  the  ice  layer  and  wall,  we  can  write  the  following 
equation  for  the  rate  of  heat  transfer  per  unit  area  of  the  wall'*' 
k.A 

Q  -  “  (32-T)  (6-37) 

£i 

where  T  is  the  temperature  of  the  rank  contents  at  any  instant  of  time. 


An  additional  equation  is  necessary  to  quantify  the  ice  thickness 
as  a  function  of  time.  Considering  an  energy  balance  on  the  ice  layer, 
we  have 

Q  »  fL  .  +  Q, 

^  ‘funicn  from  water 


This  equation  can  be  written  in  terms  of  ice  layer  thickness, 
heat  of  fusion,  etc.,  as 
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h  r  +  C  (T  -  32)]  +  h  A  (T  -  32)  (6-33) 
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The  equation  can  be  rearranged  to  give  an  expression  for  the  rate  of 
growth  of  the  ice  layer,  in  terms  of  the  temperature  of  the  tank  contents, 
as 

ki 
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(6-39) 
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+A11  properties  are  expressed  in  British  thermal  units. 
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FIGURE  S— 12  SCHEMATIC  DIAGRAM  OF  HE  ICE  FORMATION  ON 
THE  OUTER  WALL  OF  THE  SUNKEN  LNHg  TANK 
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6. 4. 3. 5  Governing  Equations 

In  obtaining  a  pressure  profile  in  time  for  the  contents  of  the 
tank,  ve  need  to  solve  simultaneously  for  the  rime  dependence  of  the 
vapor  volume  fraction,  $,  as  well  as  the  ice  layer  thickness, 

In  general,  then,  the  solution  reduces  to  that  of  a  set  of  three 
simultaneous,  first-order  differential  equations  of  the  form: 

P  -  ~r  •  f (P.  6,  1,,  P,  T,  V  ,  V  ,  h  ,...t) 
dt  l  1  v  r  v 


(6-40) 
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=  f^  ^P  >  B»  1^»  ^  J  Vv. 


(6-42) 


where  each  time  derivative  is  a  function  of  values  of  all  three  variables 
(but  not  their  time  derivatives) ,  plus  the  thermodynamic  properties  of 
ammonia  at  saturation.  The  thermodynamic  properties  may  bs  correlated 
as  functions  only  of  the  saturation  pressure,  P,  with  the  result  that 
each  time  derivative  is  a  function  of  only  the  three  dependent  variables 
and  of  time.  This  system  of  equations  ia  readily  amenable  to  solution 
on  a  digital  computer  using  one  of  many  numerical  techniques.  The 
following  section  describes  the  particular  computer  program  chosen. 


6.4.4  Solution  to  the  Gcvei 


Equations 


The  governing  equations  that  describe  the  pressure  behavior  within 
the  tank  as  a  function  of  time  involve  the  thermodynamic  properties 
of  anhydrous  ammonia.  These  physical  properties  were  obtained  from 
reference  (19)  and  are  listed  in  Appendix  A.  Unfortunately,  tabular 
data  is  not  readily  amenable  to  use  in  a  numerical  solution  of  the 
governing  equations  in  the  preceding  section;  rather,  an  analytical 
correlation  is  far  more  useful,  particularly  when  calculating  derivatives 
of  properties  along  the  saturation  curve. 

In  deference  to  computational  ease,  the  temperature,  specific  volumes 
of  the  vapor  and  liquid,  and  specific  enthalpies  of  the  vapor  and  liquid 
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«ere  correlated  by  curve-fitting  the  property  value  in  a  series  of  I 

3 

Chebyshev  polynomials  in  saturation  pressure.  This  technique  j 

is  a  standard  one  and  yields  smooth  approximations  of  the  functional  | 

dependence  of  the  property  with  respect  to  saturation  pressure.  | 

F  -  ty  was  fitted  to  the  5th  degree  in  the  Chebyshev  polynomials  j 

actory  reproduction  of  the  property  value  as  well  as  j 

p  _ . . _sion  of  smooth  derivatives  through  the  range  of  pressure  I 

(7.67  to  128.8  psia)  over  which  the  correlation  is  valid.  Accuracy  I 

4 

of  the  computer  calculated  property  values  was  checked.  Smoothness  of  | 

the  derivative  was  assured,  becacoe  the  calculated  first  derivatives  •  \ 

with  respect  to  pressure  are  monutonic,  with  the  second  derivatives  t 

(curvature)  unchanging  in  sign  over  the  range  of  pressure  for  which  the 
correlation  is  applicable.  SUBROUTINE  APROP  contains  the  coefficients  j 

yielded  by  the  curve  fitting,  and  serves  to  calculate  property  values  j 

and  first  derivatives  for  a  given  value  of  the  absolute  pressure. 

The  system  of  three  simultaneous  first-order  differential  equations 
from  Section  6. 4. 3. 5  which  describe  the  transient  characteristics  of  ^ 
submerged  liquid  ammonia  tank,  were  solved  via  a  computer  program 
invoking  one  of  the  standard  predictor-corrector  techniques  for  numerical 
integration  of  such  a  oystem.  The  program  involves  a  short  main  program 
and  four  subroutines,  in  addition  to  SUBROUTINE  APROF.  The  program 
is  extremely  flexible  in  that  it  has  the  capability  to  yary  the  (time) 

•step  size  during  the  integration  to  keep  the  numerical  errors  within' 
specified  bounds.  The  program  elements  are: 

SUBROUTINE  APf.OF  -  Calculates  property  values  and  derivatives 
given  &  pressure  ; ..  *ut . 

SUBROUTINE  FUNCT  -  Evaluates  time  derivatives  of  vapor-volume 
fraction,  pressure,  and  ice  layer  thickness  at  a  particular  time  step, 
given  the  property  valuer  from  APROP. 

SUBROUTINE  RUKU  -  Using  a  Runge-Kntta  technique,  RUKU  performs 
a  single-point  integration  of  the  three  dependent  variables,  given 
derivatives  fro s  "  NOT . 
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SUBROUTINE  KILNE  -  Using  the  Milne  predictor-corrector  method,  MILNE 
integrates  based  on  a  three-point  technique  using  derivatives  at  three 
points  to  find  function  values  at  a  fourth. 

SUBROUTINE  DE  -  This  subroutine  is  called  from  the  MAIN  program. 

Its  basic  function  is  error  control  during  the  integration.  At  the 
start  of  tne  integration,  PE  calls  FUNCT  (which  calls  APROP)  for  values 
of  the  derivatives  of  vapor  volume  fraction,  pressure,  and  ice  layer 
thickness  at  the  origin.  It  then  calls  RUKU  which  performs  a  single¬ 
point  Runge-Kutta  integration  of  this  derivative  to  obtain  the  pressure, 
etc.  at  the  first,  second,  and  third  time  steps,  DE  then  calls  MILNE 
for  a  calculation  of  the  dependent  variables  at  the  fourth  and  all 
succeeding  time  steps,  and  checks  the  error  in  the  predicted  versus 
corrected  values.  Should  this  error  exceed  that  maximum  acceptable 
error  specified  in  MAIN,  DE  will  halve  the  step  size,  employ  RUKU 
to  get  started  again,  and  proceed.  In  the  case  where  the  error  yielded 
by  MILNE  is  sufficiently  low  (as  specified  in  MAIN) ,  DE  will  double 
the  stop  size  to  speed  up  subsequent  integration.  Again,  since  the 
stop  size  is  changed,  RUKU  will  need  to  be  used  to  start  the  integra¬ 
tion  procedure  again. 

Figure  6-13  is  a  flow  chart  for  the  program  that  solves  the  governing 
equations . 

6.4.5  Results  and  Discussion 

The  results  of  the  numerical  solutions  to  the  various  cases  are 
discussed  in  this  section.  In  all  of  the  cases,  it  was  assumed  that 
the  heat  transfer  from  water  to  the  tsnk  was  over  the  entire  surface 
area  of  the  tank  (no  insulation) . 

Figure  6-14  shows  the  pressure-time  history  for  various  initial 
ullage  volumes.  After  the  pressure  in  the  tank  reaches  10  psig, 
the  relief  valve  opens.  The  important  observation  that  can  be  made 
from  the  Figure  is  the  duration  of  time  (about  6-7  hours)  it  takes  for 
the  tank  pressure  to  reach  10  psig.  This  is  indeed  an  important  result 
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FLOW  CHART  FOR  THE  NUMERICAL  SOLUTION  OF  THE  EQUATIONS 
TO  OBTAIN  SUBMERGED  TANK  PRESSURE  AS  A  FUNCTION  OF 
TEMPERATURE 
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FIGURE  6-14  PRESSURE  RISE  IN  TANK  AFTER  IMMERSION. 
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in  that  it  allows  considerable  amount  of  time  in  which  to  take  action 
following  an  accident  and  barge  sinking.  Also  notice  from  the  Figure 
that  for  larger  ullage  volumes  the  rate  of  pressure  rise  io  higher. 

This  is  explained  as  due  to  the  higher  heat  capacity  of  the  liquid 
compared  to  that  of  the  vapor.  The  larger  the  amount  of  liquid,  the 
greater  its  ability  to  absorb  heat,  and  the  slower  the  rise  in  pressure. 

The  liquid  expands  due  to  heating  and  tends  to  occupy  the  ullage 
volume.  Figure  6-15  shows  the  variation  of  the  ullage  volume  with  time 
for  different  initial  vapor  volumes.  If  the  initial  ullage  volume 
is  a  very  small  fraction  of  the  total  volume,  there  occurs  a  time  during 
the  heating  process  when  the  liquid  expands  to  completely  fill  the 
tank  volume.  This  is  shown  clearly  in  Figure  6-15  for  the  case  of  2% 
volume.  However,  in  practice  the  vapor-phase  volume  is  generally 
never  less  than  10%. 

In  the  event  the  relief  valve  fails  to  open,  the  pressure-time 
history  io  as  indicated  in  Figure  6-16a,  assuming  of  course  that  the 
tank  structure  does  not  fail.  One  can  see  that  after  the  pressure 
inside  the  tank  reaches  a  value  corresponding  to  which  the  saturation 
temperature  is  32°F,  the  ice  layer  outside  the  tank  completely  disappears, 
resulting  in  a  sudden  increase  in  heat  transfer.  This  has  detrimental 
affects  in  that  the  pressure  rise  becomes  extremely  rapid  and  eventually 
the  tank  ruptures.  The  pressure  rise  in  the  post-ic^-melt  period  is. 
shown  in  Figure  6-16b.  Observe  the  time  scales  involved  before  rupture. 

6.4.6  Conclusions 

A  theoretical  analysis  of  the  consequences  of  the  submergence  of 
a  2500- ton  ammonia  barge  has  been  made.  The  main  result  of  the  analysis 
is  the  tank  pressure-time  history.  The  analysis  has  been  carried  out 
for  the  worst  possible  case,  that  of  the  tank  inverting  after  submerging 
with  the  relief  valve  at  the  bottom.  Ala.  assumed  in  the  analysis 
was  the  total  loss  of  the  urethene  insulation  around  the  tank  and  the 
formation  of  an  ice  layer  on  the  outside  cf  the  tank  wall.  The  follow- 
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FIGURE  6-15  VARIATION 


FIGURE  6— 1®b  PRESSURE  RISE  AFTER  ICE  MELTS 


ing  conclusions  can  be  drawn  from  the  results; 

g  The  pressure  rise  inside  the  tank,  even  when  no  venting  is 
taking  place  in  the  initial  stages,  is  very  low. 

©  The  relief  valve  opening  pressure  (Id  psig)  is  reached  only 
after  6  to  8  hours,  depending  on  the  initial  ullage  volume. 

The  larger  the  initial  ullage  volume,  the  quicker  the  time  in 
which  this  set  pressure  is  reached.  Even  sc,  the  time  scale 
is  of  the  order  of  several  hours. 

©  If  the  relief  valves  open  at  the  above  pressure,  the  rate  of  ' 
pressure  rise  drops,  eventually  reaching  zero,  and  even 
negative  values.  That  is,  the  relief  vents  provided  are 
sufficient  to  cope  with  the  problem. 

•  If  the  relief  valves  are  stuck  and  do  not  open  at  all,  there 
results  a  time  (after  a  long  period)  at  which  the  pressure 
increases  very  rapidly  and  reaches  values  higher  than  the 
maximum  allowable  pressure  (45  psig)  for  the  tank,  within  a 
matter  of  minutes.  This  may  lead  to  disasterous  results. 
However,  even  in  this  case,  there  is  sufficient  lead  cime  (of 
the  order  300  to  400  hours)  within  which  remedial  measures 
can  be  taken. 

In  short,  the  design  of  the  tank  is  such  that  it  can  withstand 
several  hours  of  submergence  even  without  any  insulation.  However,  a 
word  of  caution  is  expressed  in  relation  to  the  above  analysis. 

The  whole  analysis  rests  on  the  assumption  that  a  uniform  layer  of 
ice  forma  on  the  outer  wall  of  the  tank.  If  for  some  reason  the 
Ice  does  cot  form  and  there  is  a  large  rate  of  heat  to  the  tank,  condi¬ 
tions  will  become  worse  much  sooner. 

6.5  VAP02.  DISPERSION 

6,5.1  Introduction 

Saturated  ammonia  vapor  is  lighter  than  air  and  therefore  rises 
rapidly  when  released  into  the  ataosphare.  While  rising,  the  vapor  is 


also  carried  downwind  by  the  prevailing  wind.  The  latter  phenomenon 
is  dealt  with  in  this  section.  The  parameter  of  interest  la  the 
downwind  ground  level  concentration  of  ammonia  vapor. 


Two  models  are  given,  one  of  which  is  applicable  to  the  vapor 

dispersion  in  a  small-scale  spill.  In  this,  the  lateral  dispersion 

caused  by  the  air-vapor  mixing  at  the  plume  edges  is  important  and  is 

taken  into  account.  The  second  model  is  applicable  to  large-scale 

spills  where  the  size  of  the  cloud  initially  is  so  large  that  lateral 

dispersion  can  be  neglected  in  comparison  with  vertical  mixing.  In 

both  models,  we  assume  that  the  concentrations  within  the  vapor  cloud 

can  be  represented  by  Gausjian  profilee.  Also  used  are  the  Pasquill- 
08) 

Gifford  curves  v‘  *  for  the  description  of  the  variances  of  these  Gaussian 
curves.  The  effects  of  local  topology,  variations  in  wind  velocity  and 
direction  during  the  cloud  dispersion,  meteorological  condition  varia¬ 
tions,  etc.  are  neglected.  (The  rise  of  ammonia  vapor  in  air  has 
been  modeled  and  discussed  in  Section  5.) 

6.5.2  Description  of  the  Models 

6. 5. 2.1  Continuous  Spills  (Model  1) 

The  concentration  at  any  point  downwind  can  be  represented  by 
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m(t  -  ~)  •  mass  rate  of  generation  of  ammonia  vapor 
U  -  mean  wind  velocity 
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h  *  height  of  the  center  of  the  plume  above  the  ground 

(or  water) 

2  »  coordinate  measured  from  the  ground  surface. 

m  can  be  obtained  if  we  assume  a  constant  boiling  rate  of  L!®^  on 
water  and  a  linear  increase  in  area  of  spread  vifch  time.  (See 
Section  5.)  z.Q  is  determined  from  the  plums  theory  (Section  6.7), 
while  and  a 2  are  obtained  from  the  Pasquill-Gifford  curve3  shown 
in  Figure  6-17a  and  Figure  6~17b. 

In  deriving  the  above  equation,  it  has  been  assumed  that  the  vapor 
is  generated  from  a  point  source.  Hence  the  equation  is  applicable 
at  large  distances  compared  to  the  diameter  of  the  spill. 

The  concentration  in  equation  (6-43)  ia  expressed  in  density 
unitB.  The  relationship  between  vapor  concentration  and  the  concentra¬ 
tion  in  density  units  (c)  is  expressed  in  mole  fraction  (x)  of  the  air 
vapor  mixture. 
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molecular  weight  of  air  =  28.9  gm/uole 


molecular  weight  of  vapor  =  17.0  gm/nole 


The  mass  of  vapor  collected  (M f)  by  an  iapinger  with  a  volume  auction 
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where 

m  =>  total  essks  of  vapor  liberated 
v 

and  y.z  ■*  the  cordinates  of  the  iapiugcr  with  respect  to  clcud  center 
(j  ■  0  is  on  the  ground,  and  j  **  0  is  st  the  cloud  center) 
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FIGURE  «— 17  PAffCiUSILL-GIFFORD  DISPERSION  COEFFICIENT  CURVK3  FOR  VARIOUS 
ATMOSPHERE  TYPES 


s.3.2<2 


:g?.  Spills  ( las  t,aa  can  eons  Spills) 


la  large  instantaneous  spills,  the  vapor  cloud  formed  ie  so  massive 
that  under  ordinary  wind  conditions  (less  than,  say,  10  mph)  the  mixing 
with  ambient  air  ia  slow.  Consequently,  vapors  from  larger  spills  will 
rise  to  greater  heights  than  the  vapors  from  smaller  spills,  during 
the  same  time  duration. 

The  dispersion  of  vapors  from  large  spills  depend  to  a  great  extent 
on  the  atmospheric  condition.  In  neutral  and  unstable  atmospheres,  the 
vapor  cloud  rises  continuously  as  it  disperses  downwind.  In  stable 
atmospheres  (''invers^on’‘,)  the  cloud  center  reaches  a  ceiling  height  and 
from  then  on  disperses  almost  entirely  in  the  horizontal  directions  only. 
This  will  occur  because  of  the  lack  of  atmospheric  mixing  process  above 
this  ceiling  height.  From  the  point  of  view  of  ground  level  concentrations 
o p  vapor,  the  worst  condition  during  stable  atmospheric  conditions  occur 
when  there  is  fumigation '■*'  taking  place  (that  is,  when  good  mixing 
occurs  between  the  ground  level  and  the  ceiling  height,  resulting  in  a 
uniform  vertical  distribution). 

Because  of  the  relatively  short  duration  in  which  even  massive 
(instantaneous)  spills  of  LNH.,  evaporate  and  the  comparatively  long 
time  it  takes  for  the  cloud  to  reach  distances  of  the  order  of  a 
kilometer  or  mere  (in  low  wind  velocities),  the  spill  can  be  assumed  to 
generate  all  of  the  vapor  instantaneously.  To  take  this  finiteness  of 
the  source  size  into  consideration, it  is  assumed  that  the  vapor  is 
generated  from  a  "virtual  point  source"  located  5  pool  diameters  upwind. 
Therefore,  in  the  formulae  given  below  the  dispersion  parameters  (o's) 
have  to  be  estimated  from  the  location  of  this  virtual  point  source. 

The  following  formulae  are  used  for  the  two  types  of  atmospheric 
conditions: 


(l)  neutral  and  unstable  atmospheres  (See  Figure  6-18) 
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FIGURE  S— 18 


SCHEMATIC  DIAGRAM  SHOWIf*G  THE  FEATURES  OF  THE 
VAPOR  DISPERSION  MODEL  FOR  LARGE  SPiLLS 


1,e*>  maximum  ground  level  concentration  at  any  x  is 


2m 

Vouad  -  '^pTrT'r-  * 
y  2 


where 


(6 


h  =  the  height  of  the  center  line  of  the  plume  at  downwind 

distance  x,  calculated  by  the  method  discussed  in  section 

6. 7. 2.1 

mv  »  mass  of  vapor  liberated  in  the  spill 


(ii)  Stable  Atmosphere 

There  are  several  types  of  stable  atmospheres.  An  excellent  review 
of  the  different  types  and  their  consequences  on  atmospheric  dispersion 
of  pollutants  is  given  in  reference  18.  The  different  types  arise 
because  of  the  variety  of  mixing  process  that  takes  place  in  the 
atmosphere  depending  on  the  time  of  day,  solar  insolation,  location  of 
the  area  of  concern  (over  water  Gr  land)  etc.  To  calculate  the  ground 
level  ammonia  vapor  concentrations  in  stable  atmospheres,  it  is  necessary 
to  know  the  exact  type  of  stability.  In  the  discussions  below  two 
extreme  types  of  stable  atmospheres  are  treated  and  methods  for  estimating 
the  ground  level  concentrations  are  given. 


a)  Stable:  Fanning 

In  this  type  of  stable  atmosphere  the  temperatu-:e  in  the 
atmosphere  Increases  continuously  from  the  ground  up.  The  ammonia  vapor 
released  from  a  spill  rises  to  a  certain  height  and  then  spreads 
horizontally,  with  virtually  no  vertical  dispersion.  The  ground  level 
vapor  concentration  in  such  a  case  is  very  small.  This  concentration 
is  estimated  by  using  the  same  equation  6-46br  after  substituting  the 
ceiling  height  (calculated  from  equation  6-61)  for  the  parameter  h. 
The  values  of  a  and  corresoond  to  those  in  stable  atmospheres. 
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b)  Stable;  Fumigation 

In  this  type  of  stable  atmosphere  there  ia  good  Lj.xing  between 
the  ground  and  a  certain  height  (ceiling  height) .  Above  this  height  the 
atmosphere  is  inverted.  Because  of  this  mixing  process  below  the  ceiling 
height  the  pollutants  tend  to  disperse  and  increase  the  ground  level 
concentrations.  In  the  case  of  ammonia  vapor  liberated  from  a  ground 
level  spill,  the  vapor  rises  till  it  reaches  the  ceiling  height.  Later 
on  the  dispersion  is  primarily  in  the  horizontal  direction  with  the 
vertical  concentration  distribution  being  practically  uniform  from  the 
ground  to  the  ceiling  height.  From  the  point  of  view  ground  level 
vapor  concentrations,  this  situation  is  by  far  the  worst  condition. 

Given  below  is  an  equation  to  estimate  thia  ground  level  concentration. 


m 

max  =  v 
c ground  =  2w  g  2  h 

y 

where 

h  =  ceiling  height 

Oy  =  horizontal  dispersion  parameter  for  unstable 

condition  (generally  atmosphere  C  is  used) 


(6~47) 


It  is  to  be  noted  that  the  above  ceiling  height  is  an 
atmospheric  parameter  and  has  to  be  known  before  ground  concentration 
can  be  evaluated.  The  maximum  duration  for  which  the  ground  concentration 
level  persists  above  any  tolerable  concentration  C*,  is  given  by; 
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(6-48) 


Specific  example  are  worked  out  in  Appendix  C  to  illustrate  the 
value  of  ground  level  concentrations  obtained  for  a  3000  ton  ammonia 
spill,  during  different  kinds  of  atmospheric  conditions. 
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6.5.3  Discussion 


The  dispersion  models  presented  are  often  used  for  predicting  the 
concentrations  downwind  from  a  vapor  source.  However,  most  previous 
applications  have  been  to  describe  the  dispersion  of  gases  that  are 
either  neutrally  buoyant  or  are  heavier  than  air.  In  the  sections 
jusi  presented,  we  have  attempted  to  describe  the  dispersion  of  ammonia 
vapor  -  a  lighter-than-air  gas  -  by  using  the  same  formulae,  but  including 
the  vertical  rise  of  the  vapor.  The  procedure  essentially  retains  the 
same  turbulent  mixing  characteristics  of  the  atmosphere  while  allowing 
for  the  buoyancy  driven  rise  of  the  vapor  cloud.  The  method  of  calculation 
of  the  latter  has  been  presented  in  Section  5.4.3. 


Two  models  of  vapor  dispersion  were  presented:  a  model  that  should 
be  used  to  describe  the  dispersion  of  vapor  from  a  continuous  spill, 
and  a  model  for  instantaneous  (large)  spills.  The  model  to  fc«  uspj  in 
a  particular  situation  depends  on  not  only  the  nature  of  the  opi.il  but 
also  the  relative  times  of  concentration  observation  and  evaporation 
of  spill.  The  finiter.es3  of.  the  source  area  is  accounted  for  by 
assuming  the  vapor  to  be  released  from  a  point  source  located  five  pool 
diameters  upwind  from  the  spill  point.  This  origin  shift  does  net 
significantly  affect  the  estimation  of  safe  distances  for  toxic  concen¬ 
tration  levels,  a 6  the-sa  distances  are  expected  to  be  large  compared 
to  Che  maximum  pool  sire.  However,  if  the  concentrations  are  very  close 
(within  one  ov  two  diameters)  to  the  spill  desired,  mors?  appropriate 


area  source  models  have  to  be  used.  One  of  these  methods . is  described 


in  reference  17- 


Two  extreme  examples  of  stable  atmospheres  have  been  discussed. 

The  need  to  distinguish  amongst  the  various  classes  of  stable  atmospheres 
arises  due  to  the  buoyant  nature  of  ammonia  vapor.  It  is  because  of 
this  phenomenon  together  with  the  variability  of  stable  atmospheres 
Chat  a  single  equation  can  not  adequately  describe  the  ground  level 
concentrations.  As  a  consequence  of  the  low  vertical  mixing  during 
fanning  stability  the  ground  concentrations  will  be  extremely  small. 


However,  the  concentration  predicted  by  aquation  6-47  for  the  fumigation 
case  are  high  and  should  bs  construed  as  being  extremely  conservative. 
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In  most  realistic  situations,  the  ground  level  concentrations  are 
expected  to  be  lower  (perhaps  by  as  much  as  an  order  of  magnitude)  than 
that  predicted.  This  should  be  so  because,  in  nature,  the  vertical 
concentration  distribution  is  rarely  uniform  and  is  more  likely  to  be 
peaked  at  the  ceiling  height.  Also,  there  will  be  some,  though  not  much, 
dispersion  into  regions  above  the  ceiling  height.  Together,  these  two 
phenomena  will  contribute  to  lessen  the  ground  level  concentration 
predicted  by  the  fumigation  theory. 

It  is  to  be  noted  here  that  no  test  was  conducted  during  the  present 
program  in  a  stable  atmosphere  and  consequently  the  predictions  by  the 
above  equations  could  not  be  tested  against  experiemntal  data  from  the 
present  series. 

6.6  WATER  DISPERSION 
6.6.1  Introduction 

When  liquid  ammonia  is  spilled  on  water,  a  considerable  portion 
of  the  spilled  quantity  goes  into  solution  with  water,  forming  ammonium 
hydroxide  (NH^OH).  This  NH^OH  formed  is  hot  and  lighter  than  water. 
Therefore,  it  tends  to  spread  out  on  the  water  surface  (a  spread 
velocity  of  about  0.22  ft/sec  was  observed  in  the  laboratory-scale 
spill  testa),  However,  NI^OH  is  readily  miscible  with  water  and 

thus  any  mixing  or  stirring  of  the  two  liquids  result  in  a 
diluted  ammonium  hydroxide  solution.  If  the  ammonia  spill  occurs 

in  a  river  or  an  estuary,  the  fiow-created  stream- turbulence  or  the 

wave  action  is  sufficient  to  rapidly  mix  the  ammonium  hydroxide. 

with  water  and  dilute  it. 

The  phenomenon  of  mixing  is  generally  described  theoretically 
by  the  classical  diffusion  equations  with  one  or  more  diffusion  coef¬ 
ficients.  A  comprehensive  survey  of  literature  and  the  derivations 

of  the  equations  from  first  principles  are  described  in  references  (21  ) 

(2^\  *  24  ) 

and  (  22  )  •  Several  investigators  (Fisher''  '  1 ,  Holley  and  Harleaan''  , 

(  21  ^ 

Thatcher  and  Rarleman  ' )  have  correlated,  from  experimental  data, 
the  dispersion  and  turbulent  diffusion  coefficients  for  dispersion 
in  rivers,  estuaries  and  other  regions.  The  dispersion  equations 
derived  are  strictly  correct  for  7,eutrally  buoyant,  conservative  (or  at 
best  first  order  decay)  substances.  Though  concentrated  ammonium  hydroxide 
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is  lighter  than  water,  the  density  of  water  containing  low  concentrations 
of  it  will  be  very  close  to  that  of  pure  water.  Also,  even  though  some 
vapor  is  generated  during  the  dilution  of  NH^OH,  the  quantity  of  vapor  pro¬ 
duced  is  so  small  as  to  make  the  dilution  process  a  mass  conserved  one. 

The  equations  presented  below  predict  the  concentration  of  the 
IIH^OH  solution.  The  spill  is  assumed  to  be  instantaneous.  The  two 
regions  of  water  considered  are  a  non-tidal  river  and  a  tidal 
river. 

5,6.2  Description  of  the  Models 

6. 6. 2.1  Dispersion  in  a  Non-Tidal  River 

When  dissolves  in  water  and  spreads  radially,  the  NH^OH 

formed  has  a  concentration  gradient  in  the  vertical  direction.  Figure 

6-X9  shows  the  distribution  of  NH.OH  solution  concentration  at  the 

4 

instant  the  LNH^  spread  is  a  maximum.  Dye  distribution  observed  in  the 
.laboratory  experiments  indicates  that  a  high  concentration  of  NH^OH  i.3 
at  the  water  surface,  decreasing  to  a  very  low  value  within  a  few 
inches  depth  (three  to  four  inches) .  The  water  dispersion  models 
all  assume  an  initial  uniform  concentration  (in  fact,  infinite 
concentration  for  instantaneous  spill)  condition.  For 
such  a  hypothetical,  uniform  concentration  region,  we  define 


0„7TR  ^ 

i  max 


(6-49) 


*  mass  of  J,KH^  tnat  dissolves  in  water 

P ^  density  (0.628  gms/cm^) 

R  »  maximum  radius  of  spread  of  LNH- 
max  J 

h  **  thickness  of  a  hypothetical,  100%  concentrated  rad.ius-of- 
spread  equal  to  R 


In  general,  the  value  of  h  is  very  small  (1,3  inch  for  a  3000- 
ton  apill  with  50%  partition)  compared  to  the  depth  of  rivers.  Hence 
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th»  source  can  be  essentially  treated  as  an  instantaneous  area  source. 
Two  models  are  discussed  below. 

The  diopersion  of  NR^OH  can  be  visualized  as  follows .  During  the 
initial  stages,  the  spill  pool  of  NU^OH  is  moving  downstream  at  about 
the  same  velocity  as  the  surface  velocity  of  the  river.  During  this 
period  the  NK^OH  is  getting  diluted  by  dispersion,  essentially  in  the 
vertical  direction  only.  (The  lateral  and  longitudinal  dispersions  can 
be  neglected  if  the  spill  diameter  is  of  the  same  order  as  the.  width  of 
the  river.)  After  some  time,  the  vertical  distribution  of  NH^OH 
concentration  becomes  essentially  uniform.  Beyond  this  stage  the 
dispersion  is  essentially  longitudinal.  Given  below  are  two  models  to 
describe  the  above  two  different  regions  of  dispersion.  In  the  first 
case  —  called  the  near  field  analysis  —  the  finite  area  of  the  source 
is  accounted  for.  In  the  second  case  —  called  the  farfield  analysis  — 
the  concentration  Is  assumed  to  be  uniform  throughout  the  cross  section 
of  the  river  and  dispersion  to  be  predoainently  longitudinal. 

(i)  Near  Field  Analysis: 

In  the  near  field  analysis  we  first  propose  a  model  in  which  the 
area  source  Is  taken  into  account  together  with  the  dispersion  in  all 
three  directions.  This  model  is  amenable  to  solution  only  on  a  computer 
Subsequently,  we  give  a  simplified  version  of  this  model  in  wnich  the  . 
longitudinal  and  transverse  dispersions  of  NK^OH  are  neglected. 


General  Model:  The  concentration  at  any  point  within  the  river 
at  time  t  after  the  spill  is  given  by 
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mass  of  NH,0H  par  unit:  area  of  spill  ”  —7 


2  a  e  *> 
x  y  2 


®  depth  of  the  river 
**  width  of  the  river 

*®  location  of  the  centre  of  spill  from  the 
aid  stream 

n  turbulent  diffusion  coefficients  respectively 
in  x,  y  fi:  2  (i.e.,  longitudinal,  lateral  and 
depthwise)  directions. 


The  turbulent  diffusion  coefficients  e  ,,  e  ,  and  e„  are  obtained 

x  y  ** 

from  the  characteristics  of  the  river.  Harleasn  suggests  the  use  of 
the  coefficients  given  in  Table  6-3.  The  concentration  given  in 
aquation  (&-50>  is  in  density  units. 
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Table  6-3 


Turbulent  Diffusion 

Coefficients  in 

a  River  1 

Diffusion 

Coefficients 

Very  wide  rivers 
(w/ri  >  100) 

For  narrow  rivers 
(w/d  <  100) 

Remarks 

5* 

0.067  u*d 

0.06?  u*  Rh 

The  ez  value  is  the 
mean  of  the  vertical 

Cv 

0.1  et 

0.1  eT 

distribution  given  by 

X 

e  =  u*z(l-z/d) 

cy 

0.1  e. 

0.23  u*  Rh 
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Simplified  Model:  Neglecting  the  longitudinal  and  lateral  dispersion 
and  assuming  no  edge  effects  and  using  equation  6-49,  equation  6-50  can 
be  simplified  to  the  following  equations. 


t  “  (x  -  R  )  /  U 
1  max 


C2  ■  <*  +  W  >  U 


U  «*  mean  stream  velocity 


and  for  t^  t  _<  t^ 


■i  l/ff  e  t 
ST  z 

C  »  <  _n 

water  surface  (  [*  ”  e 

P£  h/d  +  j  T.  sin  (n*  ~  )  — - 

n®!  n 


’.lit  e  t  < 

V  2 


2  2 
-n  -n 


(6-51) 


(6-52) 


ana  c  =*  o  for  t  <  t,  or  t  <  tn 

At  any  point  x  downstream  this  concentration  would  la3t  for  a 

time  of  about  2R  /  U.  It  should  be  noted  that  the  minimum  concentration 
max 

predicted  by  this  above  aquation  is  equal  to  h/d.  However,  in  the  river 

the  concentration  goes  below  this  value  due  to  the  longitudinal  dispersion. 

This  is  discussed  below. 

(ii)  Far  Field  Analysis 

In  order  to  calculate  the  far  field  concentration  value  the  following 
equation  can  be  used.  Far  field  analysis  results  can  be  used  when  in 
equation  6-52  the  term  involving  the  summation  of  the  expontially 
decaying  sine  series  is  very  amalj.  compared  to  h/d. 
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r6-53) 


where 


NH.OH  concentration  (maximum) 
4 


£  =  77  n  U 


longitudinal  dispersion 
coefficient 


t  *  time  from  the  instant  of  spill 
and  other  quantities  are  defined  in  Table  6-3. 

6. 6. 2. 2  Dispersion  in  a  Tidal  River 

When  ammonia  splJ  Is  occur  at  the  mouth  of  a  harbor  or 
in  the  esturine  regions,  one  must  consider  the  tidal  character  of  the 
waterway  in  dispersing  the  a’  monium  hydroxide  forced. 

For  a  constant  cross-sectional  area  region  with  a  sinusoidal 
variation  in  the  tidal  velocity,  superimposed  on  a  stream  velocity, 
the  following  equation  can  be  used  to  predict  the  concentration 
downstream. 


<JM  = 


[(x-Uj.+  l-i-  -^cos  easier 6)] 

-JSzT  - 


(6— o4) 


vher  e 


M  »  total  mass  of  NH^OH  that  is  dispersing 

Uf  -  constant  stream  velocity 

U^.  ■  amplitude  of  sinusoidal  wave  velocity 

6  =  phase  lag  between  spill  time  and  high  water  slack  tine 

a  *  frequency  of  oscillations 

E  a  longitudinal  dispersion  coefficient 

A  *  cross-sectional  area  of  the  channel 
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The  concentration  given  by  equation  (( 
area  mean  value*  Local  concentrations  can 
complicated  formulation  of  the  dispersion  ] 
NK^GK  generated  is  assumed  tc  be  uniformly 
cross-sectional  area  of  the  stream  at  the 
these  assumptions ,  equation  (6-54)  can  be  1 
tion  far  downstream  from  the  spill  positio; 
coefficient,  E,  is  an  empirically  detenain 
on  the  flow  characteristics  of  the  river* 
determined  from  the  following  formula: 

E  **  6u*d 

6.6.3  Diocussion 

The  character  of  navigable  waterways 
place  that  no  one  general  water -dispersion 
all  places.  We  have  presented  three  model 
rivers,  and  the  other  for  tidal  regions, 
asr.umed  that  the  primary  mixing  agency  is 
equations  can  be  used  to  predict  the  tempo 
hydroxide  concentration  in  a  stream  for  in 
Different  equations  must  be  used  for  conti 
such  equations  have  been  published  elsewhe 
values  of  the  turbulent  coefficients,  or  t 


It  is  pointed  out  here  that  during  the  present  series  of  experiments 
ammonium  hydroxide  concentrations  were  not  measured.  Consequently,  the 
results  of  the  models  suggested  cannot  be  compared  to  any  experimental 
date.  However,  these  models  have  been  well  established  in  the  literature 
and  have  been  verified  against  eidfk/.inental  data. 

4r  4  m 

6.7  VAPOR  CLOUD  RISE 
6  7.1  Introduction 


When  liquid  ammonia  spills  on  water,  rapid  boiling  ensues,  result¬ 
ing  in  the  liberation  of  a  cloud  of  NH„  vapor.  The  cloud  is  lighter 


than  air  and  rapidly  rises  into  the  atmosphere,  while  at  the  same 
time  drifting  horizontally  on  the  prevailing  wind.  Vapor  liberation 
can  be  continuous,  as  in  a  continuous  spill,  or  more  or  less  in  a  puff, 
in  the  case,  of  an  instantaneous  spill.  The  two  models  given  below 
describe  the  motion  of  a  maintained  plume  and  that  of  a  puff  of  vapor. 

The  plume  rise  model  was  developed  by  Slawson  and  Csanady  ' 

(27' 

refined  recently  by  Hoult  et  d.  The  theory  was  developed 

primarily  to  describe  the  behavior  of  chimney  plumes.  Tbe  puff  model 

(28) 

is  due  to  Morton  et  al.  ,  and  describes  the  motion  of  a  suddenly 
released  cloud  of  buoyant  gas.  In  both  the  theories,  the  motion  of  a 
uniform  density  vapor  in  a  stably  stratified  atmosphere  is  considered. 

6,7.2  Vapor  Cloud  Motion  Models 

6 . 7 . 2 . 1  Plume  Theory 

A  continuous  stream  of  buoyant  gas  released  from  a  maintained 
source  rises  into  the  atmosphere.  The  gas  mixes  with  atmospheric  air 
and  is  diluted,  resulting  in  a  continuous  decrease  of  buoyancy.  In 
the  initial  stages,  the  plume  rs  carried  vertically  by  the  momentum 
of  exit  gases.  During  the  second  stage,  buoyancy  forces  dominate 
the  plume  rise. 


If  the  atmosphere  were  of  constant  density,  theoretically  the 
plume  would  keep  on  rising.  However,  at  a  certain  height  in  a  stably 
stratiiied  atmosphere,  the  plume,  which  is  being  constantly  diluted 
due  to  turbulent  mixing  with  air,  attains  a  mean  density  that  is  equal 
to  the  outside  air  density.  Above  this  height,  the  plume  rise  is 
retarded  by  a  negative  buoyancy s  but  it  continues  to  rise  because  of 
upward  momentum.  However,  after  the  plume  has  risen  a  small  distance, 
upward  velocity  drops  to  zero  and  the  plume  falls  down,  finally 
reaching  an  equilibrium  height.  During  vertical  motion,  the  plume 
is  being  dispersed  horizontally  by  the  wind. 

The  following  equations  describe  plume  motion.  In  deriving  them, 
Slawson  et  al,^*^  and  Hoult  et  al.  assumed  that  the  rate  of 

turbulent  entrainment  of  the  ambient  air  was  proportional  to  the  verti¬ 
cal  velocity  of  the  plume.  Also,  the  plume  was  assumed  to  be  circular 
in  cross  section  and  have  uniform  density  and  velocity. 


The  following  nomenclature  describes  the  definitions  of  the  various 
symbols  used  in  the  equations. 


b 

K 

Um 


c 


d 

F 

P 

g 

£ 

£ 

c 

£ 

P 


radius  of  the  continuous  plume  at  any  position 

radius  of  plume  at  the  place  of  generation 

initial  buoyancy  volume  released  in  the  puff  *»  VgA 

fraction  of  actual  size  for  evaluating  the  position  of  the 
virtual  source  (0.005  <  c  <  0.01) 

distance  along  the  plume  beyond  which  the  initial  momentum 
has  no  affect  on  the  plume  size 

buoyancy  flux  in  the  continuous  plume  **  b^u^gd  « 

acceleration  dus  to  gravity 

buoyancy  length  (used  with  a  subscript) 


buoyancy  length  in  the  case  of  cloud  rise 

3 

buoyancy  length  in  the  plume  rise  **  Fp/U 
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R  *  radius  of  the  sperical  cloud  at  any  position 
S  »  stratification  parameter  -  ^ 

t  *  duration  over  which  c  volume  V  of  vapor  is  released 

T  *»  absolute  temperature  of  the  ambient  air 

u^  »  initial  vertical  velocity  of  vapor  at  release 

U  «  mean  wind  velocity 

V  *  volume  of  vapor  initially  generated 

x  *»  downwind  distance 

x*  *  downwind  distance  where  the  cloud  becomes  neutrally  buoyant 
z  *  vertical  coordinate 

z  a  maximum  height  reached  by  the  vapor  cloud  (or  plume) 
a  -  empirical  parameter  in  equation  (6-62)  (generally  -  0.093) 

B  «*  empirical  parameter  occur ing  in  equation  (6-56)  (0.5  <  B  <  2.0) 


fractional  density  defect  » 


R) 


—  *  absolute  value  of  the  potential  temperature  gradient  in  the 
atmosphere 

k  »  parameter  in  equation  (6-62) 


Pa  «  density  of  ambient  air  at  the  position  of  vapor  release 

Pv  =*  density  of  vapor  at  release 

a)  *  Brunt  Vaisalla  frequency  *  ^  ~ 

nJT  3z 

The  Plume  trajectory  (in  a  constant  density  region  of  the  atmosphere) 


can  be  written  as 


or  z  «  (1.5/ S  ) 


1/3  2/3  I 


(6-56) 


(6-57) 


where  $  is  an  empirical  parameter  found  to  be  in  the  range  0.5  ^  2. 
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(Fay  suggests  the  use  of  f  »  It  based  on  the  analysis  of  a  large 
volume  of  chimney  plume  data.) 

Physically,  the  buoyancy  length,  indicates  the  minimum  distance 

along  the  plume,  measured  from  the  exit,  at  which  the  entrained  fluid, 

rather  than  the  ejected  fluid,  is  balancing  the  buoyant  impulse  of 

the  exit  gases.  The  distance,  d,  along  the  plume  up  to  which  initial 

(2£ft 

gas  momentum  dominates  the  plume  rise  is  given  by  ' , 


p>H 


(6-58) 


Equation  (6-56)  was  derived  assuming  that  the  radius  of  the  "chimney" 
at  the  exit  of  the  gases  is  small.  Essentially,  the  source  was  assumed 
to  be  a  point.  However,  if  the  vapor  source  Is  an  area  source  (as  in 
the  case  of  LNH^  spill  on  water)  the  equation  can  be  easily  modified 
to  account  for  the  release  from  a  virtual  source.  The  height  above 
the  ground  of  the  center  line  of  the  plume  can  be  represented  by  the 
equation  for  the  virtual  source  by 


l..  .L+/3\i/yi.f  a  ,  £\2/3 

ls>  *p  \  [  \)  (  x) 


(6-59) 


where 


;  -  v  2(vr 

\  38ti 

and  e  “  cb^/8 


(6-60) 


c  «*  effective  radius  fraction,  for  the  finiteness  of  the  source 
(.005  <  c  <  .01). 

'//hen  the  atmosphere  is  stably  stratified,  the  limiting  height  reached 
by  the  plume  is  given  by 
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,2/3 


(6-61) 


where  S  is  the  stratification  parameter  related  to  the  potential 
temperature  gradient  in  the  atmosphere t 


Equation  (6-56)  describes  the  trajectory  of  the  plume  center 
fairly  close  to  the  Eource,  while  equation  (6-61}  describes  the 
maximum  height  reached  by  the  plume. 

6. 7.2,2  Puff  Theory 

When  a  puff  of  buoyant  gas  is  released  into  the  atmosphere,  it 
behaves  very  much  like  a  buoyant  gas  plume.  However,  in  the  case  of  a 
gas  cloud,  mixing  with  ambient  air  takes  place  over  the  entire 
surface  area  of  the  cloud.  This  leads  to  much  more  rapid  mixing  that 
that  in  a  continuous  gas  stream. 

(28) 


Morton  et  al.  '  have  derived  the  following  equations  to  describe 
the  behavior  of  a  buoyant  cloud  in  the  atmosphere.  The  cloud  shape 
is  assumed  to  be  spherical,  and  the  initial  momentum  and  size  of  the 
cloud  at  release  are  ignored.  The  cloud  trajectory,  in  the  constant  den¬ 
sity  region  of  the  atmosphere,  is 


1/4 


X‘  H  /  .  . 

[-djr]  W 


.1/: 


(6-62) 


or  in  dimensional  form 


r_J_  1 

L  2v  (a<)3  J 


1/4 


1/2  1/4  -1/2 


(6-63) 


while  the  radius  of  the  soherical  cloud  is 


3cuc 

2tt 


i'V  *  \112 


\  lJ 


(6-64) 
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where  the  value  of  the  parametric  group  ate  has  been  experimentally 
determined  to  be  in  the  range  of  0.27  <  ax  <  0.34.  Morton  suggests  the 
use  of  a  *  0.G93. 

For  the  rise  of  a  cloud  in  a  stably  stratified  atmosphere,  the 
following  formulae  are  used; 


(a<)3J 


(6-65) 


and  for  the  distance  downwind  at  which  the  vapor  cloud  first  reaches 
neutral  buoyancy  the  expression  is 

v*  tt  /k 


(6-66) 


The  vapor  cloud  reaches  the  equilibrium  height  given  by  equation 
(6-65)  after  several  oscillations.  For  practical  purposes,  the  steady- 
state  height  is  a  distance  about  five  times  the  distance  indicated 
in  equation  (6-66) . 

6.7.3  Discussion 

The  plume  model  is  used  to  describe  a  continuous  plume,  and  for 
this  model  the  height  of  the  center  line  of  the  plume  decreases 
inversely  with  the  wind  velocity  for  a  fixed  observation  position  down¬ 
wind,  The  results  presented  by  this  theory  have  been  found  to  correlate 
all  of  the  observed  data  for  chimney  plumes  from  power  plants.  The 
theory  indicates  that  the  maximum  height  reached  depends  on  both 
the  wind  velocity  and  the  temperature  gradient. 

The  second  puff  model  describes  the  ascent  of  a  suddenly  released 
buoyant  vapor  mass.  The  height  of  rise  of  the  center  of  the  cloud 
at  a  given  downwind  distance  is  proportional  to  the  inverse  square 
root  of  the  wind  velocity.  However,  the  maximum  height  reached  is 
independent  of  the  wind  velocity,  and  depends  only  on  the  initial  buoyancy 
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volume  and  the  potential  temperature  gradient  in  tha  atmosphere*  It 
(18) 

has  boen  reported v  '  that  the  result  of  equation  (6-65)  under  predicts 
by  30%  the  actual  peak  height  reached  by  clouds  generated  by  atomic 
blasts.  It  has  been  argued  by  Morten  et  al.  v  that  such  an  effect 
is  due  to  the  strong  initial  turbulent  motion  in  the  interior  of  a 
cloud  whoee  surfree  is  relatively  smooth  and  because  of  this,  there 
is  very  little  mixing  with  ambient  air.  In  effect,  the  cloud  rises 
as  a  bubble  rather  than  as  a  mixing  mass.  Similar  observations  have 
been  made  in  the  case  of  vapor  liberated  during  quiescent  conditions 
from  the  istantaneous  spill  of  50  gallons  on  LNH^  on  water. 

6.8  UNDERWATER  RELEASE  OF  LIQUID  AMMONIA 

6.8.1  Introduction 

A  liquid  ammonia  barge  that  sinks  in  water  may  eventually  release 
ammonia  underwater.  Depending  on  the  orientation  of  the  tank, 
liquid  ammonia  or  gaaeovs  ammonia  may  be  released.  When  the  liquid 
is  released,  a  plume  of  liquid  ammonia  rises  towards  the  water  surface 
and  during  this  rise  a  considerable  quantity  of  LNH^  reacts  with  water. 
The  primary  concern  of  this  section  is  to  predict  a  "critical  depth" 
deeper  than  the  depth  at  which  all  of  the  liquid  reacts  with  water 
before  reaching  the  water  surface  if  the  liquid  ammonia  is  released. 

Trying  to  model  in  all  its  entirety  the  above  process  of  rising' 
ammonia  plume  which  continuously  reacts  with  water  producing  ''apor 
(which  itself  may  dissolve  in  water  at  a  site  far  removed  from  the 
site  of  generation)  is  much  too  complicated.  The  fraction  ?f  the 
released  liquid  ammonia  reaching  the  water  surface  might  be  a  function 
of  the  release  depth,  release  rate  underwater  currents,  mixing  rate 
between  the  water  ard  the  ammonia  plume,  turbulence  level  in  the  stream 
and  the  plums,  effect  of  vapor  bubbles  on  the  mixing  process,  «.tc. 

The  physical  model  presented  below  is  based  on  the  similarity  between 
a  liquid  ammonia  plume  rising  in  watar  and  a  turbulent  diffusion 
flame.  The  model  is  based  on  simple  dimensional  arguments  and  gives? 
a  method  of  calculating  the  "critical  depth"  to  an  order  of  magnitude 
only. 
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6.8.2  Details  of  the  Model 


6. 8. 2.1  Analogy  with  Turbulent  Diffusion  Flame 


The  physical  processer  ssociated  with  the  ammonia  plume  bear 
striking  resemblance  to  tl  if fusion  flame.  If  one  conceives  of  NH^ 

as  a  fuel,  H^O  as  an  oxidii,  and  NH^OE  as  the  combusion  product, 
the  similarity  is  apparent.  Water  diffuses  into  the  LNH,  plume  forming 
a  mixture  of  NH^OH  and  NH^.  LNH^  diffuses  into  the  water  forming  a  solu¬ 
tion  of  KH^CH  and  ^0.  The  mixing  process  also  involves  s  release 
of  heat  though  the  magnitude  of  heat  release  per  mass  of  NH^  Is  clearly 
much  less  than  that  associated  with  a  combustion  process. 


The  modeling  procedure  becomes  especially  convenient  because  of 

the  analogy  with  flame.  For  the  turbulent  flame,  it  has  long  been 

known  that  the  flame  height  is  a  function  only  of  the  burner  diameter 

(  30) 

for  a  given  fuel.  The  height  is  independent  of  volumetric  flow  rate 


If  the  flame  beight  is  considered  as  the  length  of  the  reaction 
zone,  then  by  analogy  the  nixing  height  of  the  ammonia  plume  would  be 
the  height  at  which  the  ammonia  would  have  achieved  complete  solution 
in  the  water.  If  the  release  depth  were  greater  than  thi3  mixing 
height,  one  would  expect  no  ammonia  to  escape  into  the  atmosphere. 

In  the  following  section  an  order  of  magnitude  estimate  calculation  is 
given  for  calculating  the  mixing  height. 

6. 8. 2. 2  Derivation  of  Mixing  Height  Relation 


The  dimensional  arguments  for  estimating  the  height  of  turbulent 
mixing  height  etart  with  the  mass  conservation  equation: 
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(6-67) 


where 


c  “  some  specie  concentration 
D  «■  turbulent  diffusion  coefficient 

U  ■*  velocity  cf  flow  in  the  axial  direction  of  the  discharge  tube 
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z  =•  axial  co-ordinate 
z  *■  cross  axial  cooridinate 

The  solution  of  the  above  equation  will  result  in  the  definition 
of  a  "boundary  Layer  thickness  of  reaction  zone"  d,  which  is  of  the 
following  order  of  magnitude 

a  ~  0  (  ) 

Assuming  that  the  reaction  zone  dimension  at  the  mixing  height  is 
of  the  same  order  of  magnitude  as  the  pipe  diameter,  we  have 


(6-68) 
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i.e- , 
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(6-70) 


In  the  case  of  turbulent  jet  from  a  pipe,  the  turbulent  diffusivity  varies 
as 

D  o  (U  d)  (6-71) 

Substituting  6-71  in  6-70,  we  have 

h  \*  0  (  d  )  (6-72) 

Of  course,  application  of  the  turbulent  flame  concept  requires  that 
the  flow  from  the  pipe  be  turbulent.  Calculations  of  the  liquid 
ammonia  Reynolds  number  at  the  exit  of  a  two-inch  diameter  pipe  for 
release  from  a  10  psig  tank  indicate  the  flow  is  quite  turbulent. 

Therefore,  we  can  write  a  relation  for  the  plume  mixing  height  "h"  as 

h  =  A  d  +  B  (6-73) 

where  A  and  B  ate  experimentally  determined  constants. 
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6.8.3  Discussion 


The  result  derived  in  equation  (6-72)  indicates  that  the  length 
of  the  LNH^  plume,  before  it  is  completely  converted  to  NE^OH  is 
of  the  same  order  of  magnitude  as  the  diameter  of  the  pipe  from  which 
it  is  released.  Within  the  scope  of  the  assumptions  made,  chie  is  an 
interesting  result.  It  implies  that  if  the  mouth  of  the  discharge  pipe 
is  at  a  depth  greater  than  about  5  to  10  diameters  of  the  release  pipe, 
all  of  the  LNH„  would  be  converted  to  NH^OH  and  no  vapor  would  be 
liberated . 

Also  the  analysis  presented  has  tacitly  assumed  that  the  turbulent 
mixing  in  the  jet  is  characterized  by  the  turbulence  in  the  pipe.  How¬ 
ever,  this  may  not  always  be  true.  During  the  discharge  of  LNE^ 
into  a  stream,  the  stream  turbulence  itself  may  be  very  important. 

Also  to  be  considered  is  the  turbulence  caused  by  the  rapid  boiling 
phenomenon  during  the  heating  of  LNH^  by  water. 

Many  other  phenomena  that  may  occur  in  an  underwater  release 
have  not  been  considered.  These  include  possible  explosions,  limita¬ 
tion  cf  water  quantity,  water  currents,  etc.  Further  refinements 
of  the  model  will  have  to  include  one  or  more  of  these  phenomena. 


•7 .  NEUTRALIZATION 


7.1  INTRODUCTION 

The  general  characteristics  of  a  spill  cf  refrigerated  (-28°F) 
ammonia  on  water  may  be  described  as  follows.  A  spreading  pool  of  liquid, 
boiling  ammonia  forms  cn  the  water.  About  half  to  three-fourths  of  the 
total  amount  of  ammonia  spilled  dissolves  in  the  water  by  the  time  the 
liquid  has  completely  boiled  away,  forming  a  layer  cf  ammonium  hydroxide 
which  Is  elevated  in  temperature  between  55F  and  15°F  above  the  ambient 
water  temperature.  After  the  liquid  aasonia  has  disappeared,  more 
ammonia  will  be  generated  from  the  ammonium  hydroxide  layer  and  vaporize 
into  the  atmosphere.  However  the  quantity  of  vapor  produced  is  small 
compared  to  the  amount  of  vapor  liberated  during  the  boiling.  The 
thickness  of  the  ammonium  hydroxide  layer  is  in  the  neighborhood  of 
a  few  inches,  and  the  rate  of  spread  is  roughly  a  few  inches  per  second 
(radial  increase),  depending  on  the  rate  of  spill.  The  maximum  spread 
radius  is  proportional  to  the  three  eights  power  of  the  quantity  of 
spill.  For  example,  if  3,000  tons  of  refrigerated  ammonia  were  spilled 
instantaneously,  the  rate  of  spread  would  be  2-1/2  inches  per  second, 
the  maximum  thickness  of  the  ansnonium  hydroxide  )ayer  would  be  about 
5  inches,  and  the  maximum  radius  of  the  liquid  pool  would  be  approxi¬ 
mately  475  feet.  The  characteristics  cf  underwater  leaks  are  determined 
by  the  leak  rate  and  depth  under  water. 

The  dispersion  of  the  ammonia  vapor  in  the  atmosphere  depends  on  the 
wind  velocity  and  atmospheric  conditions.  Because  of  the  extremely  buoyant 
nature  of  the  ammonia  vapor,  there  results  a  rapid  rise  of  the  vapor 
cloud,  following  the  spill.  It  is  estimated  that  within  about  1000 feet 
of  downwind  travel,  the  center  of  the  vapor  cloud  from  a  3000  ton  spill 
would  reach  a  height  of  between  1500  to  1800  ft. 

Estimate  of  the  time  of  vaporization  of  a  3000  ton  instantaneous 
spill  indicates  that  the  entire  boiling  process  lasts  only  for  about 
120  seconds.  If  therefore  the  time  c-f  spilling  of  this  3000  tona  is 
larger  than  about  10  minutes,  it  should  be  assumed  that  evaporation  of 
the  spilled  mass  occurs  as  soon  as  the  spill  touches  the  water. 


The  approaches  which  can  be  taken  towards  corrective  action  mcy 
be  enumerated  as  followe: 

A.  Ammonium  hydroxide  solution  in  water 

1.  Diopersal 

2.  Containment 

3.  Air  sparging 

4.  Sorption 

5.  Chemical  neutralization 

B.  Ammonia  vapors 

1.  Dispersal 

2.  Watt,  fog 

7.2  AMMONIUM  HYDROXIDE  SOLUTION 

The  preferred  method  of  treatment  for  the  ammonia  dissolved  in  tha 
water  depends  on  the  specific  situation.  If  the  principal  hazard  is 
due  tc  the  presence  of  ammonium  hydroxide  in  the  water,  and  the  presence 
of  ammonia  \*apor  in  the  air  is  cf  secondary  Importance,  containment  of  the 
ammonium  hydroxide  while  allowing  evaporation  of  ammonia  from  the  water 
to  occur,  or  stripping  of  the  ammonia  from  the  water  by  subsurface  air 
sparging,  may  be  desirable.  If  the  reverse  ie  the  case  and  the  signifi- 
cant  hazard  Is  ammonia  vapor  in  the  air,  it  is  important  to  prevent  evapora¬ 
tion  by  dilution,  neutralization  or  absorption  cf  the  ammonium  hydroxide 
In  the  water. 

7.2.1  Dispersal 

If  the  spill  is  on  a  river  or  any  body  of  water  where  rapid  water 
movement  occurs,  the  NIL  OH  in  the  water  will  be  carried  away  with  the 
current.  If  mixing  occurs,  which  ia  likely  under  such  conditions,  the 
dispersal  will  bo  enhanced  by  occurring  in  depth  ss  well  as  horizontally, 
particularly  in  the  case  where  there  Is  significant  wave  action. 

The  effects  of  the  NH^OH  will  be  felt  along  the  entire  swept-out 

path,  although  the  severity  of  these  effects  may  decrease  with  distance 

downstream  due  to  the  reduction  of  the  NH-CB  concentration  as  a  result 

t* 

of  evaporation  of  NH^  from  the  water  surface  and  the  mixing  with  larger 
volumes  of  uncon tasinated  water. 
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la  the  case  of  a  confined  channel,  however,  a  plug  of  NH^OH  may 
sweep  down  the  river  with  only  a  slight  decrease  in  concentration  with 
downstream  distance. 

in  lakes  or  harbor  area3  where  currents  may  be  small,  dispersal 
may  be  enhanced  by  running  small  vessels  through  the  spill,  utilizing 
the  turbulence  caused  by  the  propellers  and  wake  of  the  vessel  to  pro¬ 
mote  mixing  with  the  water.  Care  should  be  taken,  however,  to  temporarily 
block  cff  v/ater  intakes,  and  it  should  be  remembered  that  the  NH.OH  is 
highly  corrosive  to  copper,  copper  alloys,  and  galvanized  surfaces. 

Smaller  spills  in  still  bodies  of  water  may  be  dispersed  with  fire 
hoses  (or  any  high  pressure,  high  volume  water  hose)  played  onto  the 
NH^OH  layer  from  dock  or  shipboard  facilities. 

Since  prompt  action  is  required  for  dispersal  of  an  LNH^  spill, 
the  use  of  water  hoses  or  the  turbulent  wake  of  a  boat  or  vessel  is 
the  only  practical  means  of  achieving  rapid  dispersal. 

7.2.2  Containment 

In  lakes,  ponds,  or  other  less  active  bodies  of  water,  the  NH^OH 
layer  may  remain  intact  near  the  spill  location.  If  it  is  important 
to  localize  Che  adverse  impact  due  to  the  NK^OE  on  aquatic  life 
floating  boom  such  as  is  used  to  contain  oxl  spills  may  be  deployed 
to  contain  the  Nfi^OH  layer.  The  instantaneous  spill  considered  above  does 
not  represent  a  realistic  situation.  In  most  cases,  a  3,000  ton 
spill  would  be  the  result  of  a  tanker  barge  or*  vessel  accident,  and 
the  ammonia  would  most  likely  be  released  over  a  period  of  time  rang¬ 
ing  from  minutes  to  hours.  Consequently,  the  maximum  spill  diameter 
would  be  substantially  less  that  1,000  feet.  The  EH.  OH  spreading 
could  probably  be  contained  by  2,000  feet  of  boom.  Placement  of  a  boom 
would  allow  a  substantial  fraction  of  the  ammonia  to  evaporate  from  the 
NH^OH  layer,  and  would  allow  the  use  of  a  neutralizing  agent  or  cf 
absorption  or  removal  or  the  NH^OH  or  Nrl^. 

There  are  a  large  number  (more  than  65)  of  different  types  or 
floating  booms  for  containment  of  oil  spills,  some  of  which  would 


be  suitable  for  containing  an  LKH^  spill.  In  the  case  of  an  LNH^  spill, 
it  Is  important  to  use  a  boon  which  is  rapidly  deployable  and  is  con¬ 
structed  of  oaterials  which  are  inert  to  NH.GH  and  NH,.  It  has  been 
observed  that  the  currents  in  sheltered  areas  where  cost  spills  occur, 
such  as  at  piers,  in  ships,  at  bulkheads  and  in  coves,  are  ouch  smaller 

than  the  high  velocity  currents  which  occur  in  aidstreao  at  peak  flow 
(32) 

periods.  It  was  found  that  a  boom  having  a  freeboard  barrier  of  6  inches 
and  a  subtaerged  skirt  ox  12  inches  would  suit  many  spill  situations  and  be 
rapidly  deployable. 

The  following  general  performance  characteristics  are  desirable 
for  an  containment  boon: 

1.  The  boca's  fabric  should  contain  the  flotation  pods,  ballast 
and  vertical  stiffeners,  and  be  capable  of  withstanding 
total  tension  stresses  up  to  6,000  lbs  without  additional 
chain  or  cable  tension  aeabers. 

2.  The  fabric  Eust  be  resistant  to  NH^OE,  abrasion,  solar 
heat  and  ultraviolet  radiation,  salt  water  and  anticipated 
environmental  temperature  extremes . 

3.  The  boon  oust  be  towcble  from  one  end  at  speeds  up  tc  5  knots 
without  damage . 

4.  The  boots  should  weigh  less  than  2  lbs  per  running  foot  and 
cost  less  than  $12  per  foot  with  a  ainisuB  life  expectancy  of 
five  years.  It  should  be  bright  yellow  in  color  to  reduce 
accidental  damage  by  uninformed  vessels  in  harbors.  Grab 
handles  should  be  provided  on  the  top  edge. 

5.  The  ballast  and  float  configuration  should  give  stability 
in  15  knot  winds,  0.5  to  1.0  knot  currents  and  2  ft  waves. 

Although  none  of  the  canufacturers  of  oil  containment  booms  appear 
to  meet  all  of  the  above  specifications,  several  of  thea^^  could  be 
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considered  for  containment  of  LNH^  spills:  Acme  0.  K.  Corral  Boon; 

Bennett  20-inch  Harbor  Oil  Boom;  Jaton  3oom,  Centri-Spray  Corp.;  Reynolds 
Aluminum  Oil  Boom;  Sea  Curtain,  Kepner  Plastics  Fabricators,  Inc. 

In  cases  of  currents  above  the  retention  capability  of  the  floating 
boom,  it  rosy  be  desirable  to  divert  the  LNH^  spill  to  keep  it  from 
vulnerable  shoreline  locations  or  to  guide  it  into  an  area  where  it  can 
be  controlled  and  further  corrective  action  can  be  taken. 

Mooring  of  the  floating  boom  should  not  be  as  critical  for  LNH^ 
spills  as  for  oil  spills,  since  the  duration  of  the  spill  is  limited 
by  evaporation,  and  the  full,  twelve  hour  tidal  variation  will  probably 
not  have  to  be  accommodated  by  the  mooring. 

Bubble  barriers  or  air  curtains  are  not  usually  rapidly  deployable 
and  are  not  likely  to  be  available  in  time  to  contain  LNH^  spills.  The 
effect  of  such  barriers  would  be  to  enhance,  evaporation  of  NH^  from  the 
water  into  the  air. 

7.2.3  Air  Sparging 

Once  the  LNH^  spill  has  been  contained,  it  may  be  allowed  to  evaporate 
from  the  water  over  a  period  of  time,  or  a  more  active  course  may  be 
taken.  If  it  is  important  to  remove  the  contained  ammonia  from  the  water 
rapidly  and  there  is  no  significant  hazard  associated  with  transferring 
the  to  the  atmosphere,  air  sparging  may  be  used- 

Ammonia  can  be  stripped  from  water  through  the  use  of  large  quantities 
of  air.  Since  the  NH^OH  layer  is  only  a  few  inches  deep  and  the  power 
required  for  air  ejection  is  directly  proportional  to  the  depth  of  air 
ejection,  the  air  bubble  pipes  should  be  suspended  only  a  few  feet  below 
the  surface.  For  a  four  foot  pipe  depth,  2  psig  air  is  required,  and 
the  horsepower  requirement  is: 

HP  *=  0.024  Q 

where  Q  is  air  volume  in  cfm  (cubic  feet  per  minute) ,  and  compressor 
efficiency  and  pipe  and  fitting  losses  have  been  accounted  for.  For 
example,  if  the  NH^OH  layer  were  4  inches  deep  and  475  feet  in  radius 
(3,000  ton  spilled  instantly),  a  10b  HP  internal-combustion-engine  driving 
a  compressor  unit  would  pass  one  cubic  foot  of  air  through  each  cubic 


foot  of  N»i^OH  every  56  minutes.  This  type  of  self-contained 

compressor  unit  is  common  in  construction  work. 

7.2.4  Sorption 

In  some  eases  it  is  important  to  remove  the  ammonia  from  the  water 
without  releasing  it  into  the  atmosphere.  One  way  of  accomplishing  this 
is  by  absorption  or  adsorption  with  solids.  Absorption  refers  to  the 
containment  of  the  ammonia  in  the  bulk  of  the  aolid;  adsorption  refers 
to  the  adherence  of  the  ammonia  on  the  surface  of  the  solid.  Both  methods 
are  included  in  the  term  sorption. 

At  least  half  a  dozen  different  sorbent  booms  have  been  developed 
for  oil  spill  removal;  similar  principles  can  be  applied  to  ammonia.  The 
advantage  of  configuring  the  sorbent  in  the  form  of  the  boom  is  two-fold: 
it  may  be  towed  over  the  spill  repeatedly,  which  effectively  circumvents 
the  problem  of  dispensing  the  sorbent;  and  it  may  be  placed  on  board  ship 
or  towed  to  a  clean-up  area  after  use,  which  reduces  the  difficulty  of 
collecting  the  used  sorbent.  Since  the  depth  of  the  NH^OH  layer  is  only 
a  few  inches,  the  sorbent  boom  will  penetrate  the  entire  layer  and  practi¬ 
cally  all  of  the  NH^OH  should  come  into  contact  with  the  sorbent  material. 

When  the  sorbent  boom  is  saturated  with  ammonia,  it  must  be  removed 
to  a  clean-up  area  where  the  ammonia  can  be  removed  and  the  boom  be  readied 
for  further  use  as  required.  Ammonia  removal  may  be  accomplished  by 
washing  the  boom  with  water,  squeezing  the  sorbent  material  to  remove  • 
NH^OH,  allowing  the  NHg  to  evaporate  into  the  atmosphere,  or  heating  the 
sorbent  material  to  promote  evaporation  of  ammonia.  If  the  sorbent 
material  is  low  cost  and  readily  available,  it  may  be  disposed  of  in 
sealed  containers  or  by  chemical  treatment. 

The  sorbent  material  in  sorbent  booms  is  generally  contained  in 
porous  bags,  cloth  or  netting.  The  sorbent  materials  must  have  a  preference 
for  absorption  or  adsorption  of  KR^  or  NH^OH  rather  than  water. 

The  practical  difficulties  in  the  use  of  a  sorbent  for  large  spills 
lies  in  the  quantity  of  sorbent,  or  number  of  sorbent  boom  deployments 
required.  Sorbent  capacity  is  generally  stated  in  terms  of  a  retention 
ratio,  or  ratio  of  weight  of  taken  up  to  weight  of  sorbent.  With  a 


retention  ratio  of  2,  a  3,000  ton  spill  which  dissolved  2,000  tons  in  the 
water  would  require  1000  tons  of  sorbent  or  500  passes  with  a  2  ton  sorbent 
boon.  For  this  reason,  the  use  of  a  sorbent  may  only  be  practical  for 
small  spills.  It  should  be  noted,  however,  that  the  sorbent  is  generally 
inactive  chemically,  and  should  not  present  a  hazard  in  storage,  handling, 
or  transport. 

7.2.5  Chemical  Neutralization 

The  NH/+0H  layer  may  be  chemically  neutralized  by  the  addition  of  an 
appropriate  chemical  to  the  water.  However,  there  are  major  difficulties 
associated  with  this  type  of  corrective  action,  and  these  are  discussed  . 
below. 

Quantity.  The  neutralizing  agent  must  be  used  on  a  mole  for 
mole  basis  with  the  ammonia  to  completely  neutralize  it.  For  a  3,000 
ton  ammonia  spill  of  which  only  2,000  tons  dissolves  in  the  water,  4,300 
tons  of  HC1  would  be  required.  Neutralizing  chemicals  are  not  normally 
available  in  these  quantities. 

Dispensing  method.  The  neutralizing  chemical  must  be  dispensed 
in  a  uniform  and  controlled  manner  throughout  the  NH^OH  layer.  If  this  is 
not  accomplished,  additional  neutralizing  agent  may  be  required  and  con¬ 
centrations  of  excess  neutralizing  agent  may  occur  and  cause  additional 
adverse  conditions  in  the  body  of  water.  It  is  difficult  to  dispense 
the  agent  uniformly  dispensed  over  the  water  surface,  it  must  still  mix 
or  diffuse  through  the  NH^OH  layer  to  achieve  optimum  effect;  uniform 
mixing  may  be  Inhibited  by  density  differences  between  the  neutralizing 
agent  and  NH^OH  or  by  the  kinetics  of  the  chemical  reaction. 

Heat  generation.  Most  of  the  neutralizing  agents  that  could  be 
used  with  NH^OH  would  cause  exothermic  reactions,  and  create  a  significant 
temperature  rise  in  the  surface  layer  of  water.  Although  some  of  this 
heat  may  be  removed  by  evaporation,  convective  air  cooling  and  radiation, 
some  of  the  heat  will  be  removed  by  water  convection.  The  adverse  effects 
of  the  increase  in  water  temperature  must  be  taken  into  account. 


Reaction  products.  The  most  likely  product  of  a  neutralising 
reaction  is  a  salt,  i.e.:  NH^Cl  if  HC1  is  used.  Most  of  the  ammonium 
salts  that  would  result  from  use  of  a  neutralizing  agent  are  soluble 
in  water,  and  would  be  difficult  to  remove  from  a  large  body  of  water. 

The  adverse  effects  of  the  dissolved  salts  must  be  considered  in  the  use 
of  a  neutralizing  agent. 

7.2.6  Safety 

Most  neutralizing  agents  are  hazardous  chemicals  which  must  be 
stored,  transported  and  handled  with  caution.  The  circumstances 
surrounding  an  LNH^  spill  allow  little  time  for  adequate  practice  of 
safety  procedures;  the  use  of  a  neutralizing  agent  can  represent  a 
new  hazard  comparable  in  potential  adverse  impact  to  the  original 
LNH3  spill. 

It  appears,  from  the  above  discussed  difficulties,  that  the  use  of 
a  neutralizing  agent  is  not  the  preferred  corrective  action  for  an  LNH^ 
spill,  except  perhaps  for  very  small,  confined  spills. 

7.3  AMMONIA  VAPORS 

Depending  on  the  particular  situation,  it  may  be  dv jirable  to  disperse 
the  ammonia  vapor  into  the  atmosphere,  or  conversely,  to  trap  the  ammonia 
vapor  by  water  (or  other  means)  and  remove  it  from  the  atmosphere.  For 
a  3,000  ton  spill  of  LKH^  on  water  during  a  one  hour  period,  about  1,000 
tons  of  ammonia  will  be  generated  during  a  period  of  120  seconds.  Subse¬ 
quent  to  the  complete  evaporation  ot  the  liquid  ammonia  pool,  provided 
no  corrective  action  is  taken  with  respect  to  the  NH^OH  in  the  water, 
ammonia  will  continue  to  evaporate  from  the  water  surface  to  the  extent 
of  1,000  tons  during  the  nexjt,  few  minutes,  for  quiescent  water  conditions. 

7.3.1  Dispersal 

If  a  strong  wind  is  blowing,  the  ammonia  vapor  will  be  dispersed 
into  the  atmosphere,  leaving  little  alternative  as  to  any  possible 
correction  actions.  This  is  true  regardless  of  wind  conditions  since 
the  cold  ammonia  vapor  is  always  lighter  than  air.  Little  can  be  done 
to  prevent  rapid  dispersal  of  the  vapor  cloud. 
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7.3.2  Water  Fog 


If  atmospheric  conditions  are  such  that  the  ammonia  vapor  is  not 
rapidly  and  naturally  dispensed,  it  may  be  desirable  to  attempt  to 
remove  it  from  the  atmosphere.  Probably  the  tost  affective  way  of 
achieving  this  goal  is  by  the  use  of  a  water  fog  or  water  spray.  Equip¬ 
ment  for  generating  substantial  quantities  of  a  dense  water  fog  could 
be  developed  utilizing  the  principle  of  the  cooling  (to  near  ambient 
temperature)  of  a  warmer,  water  saturated  air  mass,  and  using  mechanical 
blowers  or  fans  to  direct  the  fog  through  the  ammonia  vapor.  The 
result  would  be  the  formation  of  a  NH^OH  fog  which,  as  it  cooled,  would 
condense  and  fall  into  the  water  surface  where  further  corrective  action 
could  be  taken. 

If  specialized  water  fog  generating  equipment  is  not  available, 
fire  hoses  could  be  utilized  (with  spray  nozzles,  if  possible)  to  play 
streams  of  spray  or  water  through  the  ammonia  vapor  cloud.  This  would  not 
be  as  effective  as  a  water  fog,  and  might  create  the  additional  problem 
of  forming  a  liquid  ammonia  aerosol  by  splashing  the  LNH^  (if  the  liquid 
ammonia  pool  has  not  yet  completely  evaporated) ,  which  is  more  difficult 
to  remove  with  the  water  spray  than  is  the  ammonia  vapor. 

7,3.3  Sorption  and  Chemical  Neutralization 

The  use  of  solids  to  absorb  or  adsorb  the  ammonia  vapor  may  be 
■unreasonably  difficult.  The  only  possibility  would  be  the  use  of  a 
pulverized  or  dust-like  sorbent  which  could  be  projected  through  the 
ammonia  vapor  cloud  by  a  mechanical  blower.  However,  large  quantities  of 
sorbent  would  be  required  (on  the  order  of  500  to  1,000  tons  for  a 
3,000  ton  spili) ,  and  the  problems  of  availability,  transport,  and 
subsequent  clean-up  of  the  sorbent  from  the  water  appear  to  be  major 
objections  to  this  type  of  corrective  action. 

The  use  of  a  chemical  neutralizing  agent  is  even  more  unlikely 
considering  the  difficulties  discussed  above  in  Section  7.2.5.  The 
hazards  associated  with  spraying  large  quanticies  of  a  neutralizing 
chemical  through  the  ammonia  cloud  appear  to  be  greater  than  the 
benefit  that  might  be  gained  by  neutralizing  the  ammonia  vapor. 
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APPENDIX  A 


This  appendix  contains  the  thermodynamic  properties  of  anhydrous 

ammonia.  Table  A-l  (taken  from  reference  33).  gives  the  important 

property  values.  Tables  A-2,  A-3  and  A-4  (from  reference  19)  give 

respectively  the  boiling  points  and  the  liquid  and  vapor  enthalpies, 

all  as  functions  of  pressure.  Figure  A-l  is  the  enthalpy  concentra- 

(1S) 

tion  diagram  for  an  ammonia-water  system. 
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Table  A~1  Important  Properties  of  Anhydrous  Ammoftia 


Chemical  formula . 

Molecular  Weight _ _ 1. 

Boiling  Temperature  at  Atmospheric  Pressme,  F 

Freezing  Temperature  at  Atmospheric  Pressure,  F _ 

Critical  Temperature,  F - ... _ _ _ _ .. _ 

Critical  Pressure,  psia  _ _ _ _ _ _ _ _ 

Critical  Density,  lb  per  cu  ft 

Density  of  Liquid  at  -28.0  F,  lb  pe$  cu  ft , _ _ _ 

Specific  Volume  of  Saturated  Vapor  at  -28.0  F,  cu  ft  per  lb. 

Specific  Heat  of  Liquid  at  86  F,  Btu  per  (lb)  (F)... . 

Specific  Heat  Ratio  (cj)/cv)  of  Vapor  at  86  F  and 

One  Atmosphere  Pressure  _ _ _ _ 

.Thermal  Conductivity,  (Btu)  (ft)  pet  (sq  ft)  (hr)  (F): 

* 

Saturated  Liquid  at  32  F  ... 

Saturated  Liquid  at  86  F  . . 

Vapor  at  One  Atmosphere  Pressure  at  32  F  . 

.  Vapor  at  One  Atmosphere  Pressure  s‘  86  F  . . 

Yiscosity,  Centipoises: 

Saturated  Liquid  at  5  F  _ _ _ 

Saturated  Liquid  at  86  F  . . . 

Vapor  at  One  Atmosphere  Pressure  at  5  F  . . 

Vapor  at  One  Atmosphere  Pressure  at  SS  F  . 

Relative  Dielectric  Strength  of  Vapor  at  84  F 

and  One  Atmosphere  Pressure  (Nitrogen  =1)  . . 

Color  . . . . . 

Odor  and  Detection  . . . . 


nh3 

17.03 


Flammability  Limits  (Pet  cent  by  volume) . 

Toxicity,  Underwriters'  Laboratories  Classification  .... 


0.29 

ase  (Avg.) 

0.0128 

0.0145 

0.250 

0.207 

0.0085 

0.0102 

0.82 

Clear  and  water  white 
Irritating  odor,  readily 
detectable  by  smell. 
Bun.ing  sulphur  used  for 
locating  leaks. 

16  to  25 
Group  2 
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RIATTVE  OrtHAW. 


FIG  A-t  ENTHALPY-CONCENTRATION  DIAGRAM  FOR  AMMONIA-WATER  MIXTURES 


Reproduced  from 
basf  available  copy. 
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Table  A-2.-B0IUNG  POINTS  OF  AMMONIA-WATER  MIXTURES  AS  A  FUNCTION  OF 

SATURATION  PRESSURE 


1 

101.7 

49.5 

12  5 

-16.5 

-40.5 

s 

126.1 

72-5 

35-0 

*•5 

-21.5 

4 

153.0 

99.0 

61.0 

29-5 

0.5 

6 

170.1 

116.0 

77-5 

*5-5 

16.0 

3 

182-9 

128,5 

89.0 

if:§ 

26.5 

ID 

1S3-2 

138. 0 

90-0 

34.5 

12 

202.0 

146.5 

107-0 

73.0 

42.0 

15 

213-0 

158.0 

H7.5 

830 

20 

228.0 

173.0 

132.0 

97.0 

§.5 

25 

240.1 

184.5 

1*3.5 

108.0 

fc? 

30 

250.3 

195-5 

153.0 

117.0 

35 

259-3 

204-5 

162.0 

125-5 

92.0 

40 

267.3 

213.0 

169-5 

132-5 

98.5 

50 

281.0 

225.0 

182.5 

145.0 

m.o 

6o 

292-7 

238-5 

248.0 

193-5 

156.0 

121.5 

s 

302.9 

20?.  5 

165.5 

130.5 

312.0 

253.0 

267.0 

J»12.5 

173.5 

138.5 

so 

320.3 

22C.  5 

181.5 

146.0 

loo 

327-8 

270-0 

229-5 

133.5 

153.0 

no 

33*- 8 

282.5 

236.0 

195-0 

158-5 

120 

341.3 

289.0 

242.0 

201.0 

164.0 

nc 

3*7-3 

295.0 

243.0 

206.5 

169.0 

174.0 

l*o 

353-0 

35$.* 

301.0 

25*- 0 

212.0 

150 

306-5 

259-5 

217.0 

179-0 

160 

353-5 

352.0 

254.5 

222.0 

183-5 

3S 

360.4 

317.0 

270. 0 

226.5 

187.5 

373-1 

322.0 

274.5 

230-5 

192-0 

ISO 

200 

Hl:i 

327-0 

33L-5 

ISio 

23*.  5 
238-5 

196- C 
200.0 

210 

385-9 

335-5 

238.0 

242.5 

203-5 

220 

389-9 

339-5 

291-5 

245.0 

207-5 

2J0 

393-7 

343.0 

295-0 

249-5 

211.0 

240 

397-4 

3*6.5 

290.5 

253. 0 

214.0 

250 

401.0 

350-0 

301.5 

256-5 

217-5 

260 

404.4 

353-0 

304.5 

259-5 

220-5 

270 

280  • 

407-3 

356.0 

JC7.5 

262.5 

224.0 

4U.1 

359-0 

310-3 

266.0 

227.0 

290 

414.2 

362.0 

313-0 

268. 5 

230.0 

300 

417-3 

364.5 

316.0 

271-5 

232-5 

ss 

420.4 

367.5 

319-0 

274.5 

255-5 

423-3 

370.0 

321-5 

277-5 

250.0 

?jS.o 

330 

*26.2 

372.5 

324.0 

2*1.0 

>40 

429-0 

37J-5 

326.5 

282.5 

243.5 

350 

431-7 

378.0 

329.0 

285.0 

£46.0 

360 

434.4 

380.5 

331.5 

267.5 

2*6-5 

*57;  o 

383-0 

334.0 

290.0 

251.0 

439-6 

336-5 

292-5 

253-5 

SS 

442.1 

339.0 

295-0 

255-5 

444.6 

390.5 

341.5 

297-5 

230. c 

410 

447.0 

392-5 

3*3-5 

299-5 

260.5 

4;o 

440.4 

395-0 

>46.0 

302.0 

262.5 

4  JO 

4W 

451.7 

397-5 

3*6.5 

304.5 

264.5 

454.0 

399-5 

402.0 

351-0 

306.5 

267-0 

*50 

456.3 

353-0 

309-0 

269-0 

460 

458.5 

404.5 

355-5 

311-0 

271.0 

S 

460.7 

407.0 

357-5 

313-0 

27>.C 

462.8 

400.0 

315-5 

275-5 

490 

464.3 

411.5 

317-5 

277-5 

500  , 

467-0 

414.0 

364.0 

319-5 

279.5 

-64.0 

-83-5 

-95.0 

:H:§ 

-102-5 

-105.0 

-45.5 

-65.0 

-76.5 

-85.5 

-S0.2 

-24.5 

-4*.  5 

-55.O 

-62. 5 

-66-5 

-69-6 

-10.0 

-3-1-0 

-43.O 

*50.5 

-55.O 

-57.6 

O.C 

-21.5 

-33-5 

-41.0 

-46-0 

-48.5 

8.0 

-13.5 

-26.0 

-53.5 

-35.5 

-*1.3 

15.  C 

-6.5 

-19.0 

-27.0 

-32-5 

-35-2 

24-5 

2.0 

-10.5 

-16.5 

-24.0 

-27.  J 

36.5 

14.5 

3U5 

-7-5 

-33-0 

-16.0 

46.5 

2*-5 

11.0 

1-5 

-4.5 

-8.0 

55.0 

32-5 

18-5 

9-5 

3-5 

-0.6 

52-5 

40.0 

25.5 

16.5 

If  0 

5-9 

69.O 

46.5 

32.0 

22-5 

.1  .0 

11.7 

80.5 

5S.0 

43.0 

33.0 

vJ.O 

21-7 

90.5 

67.5 

52.5 

42.0 

35-0 

30:2 

99-5 

76.0 

60.0 

*9-5 

VJ.O 

37.7 

IO7.5 

83-0 

67.5 

57.0 

49.0 

44.4 

II4.5 

90.0 

S:i 

§3.5 

55-5 

50.5 

56.0 

121-5 

97.0 

69-5 

6L5 

126.5 

102.5 

86.0 

74.5 

37-5 

61.2 

132.5 

I08.5 

91.0 

S.1 

89.5 

71-5 

66.  C 

137.5 

142.5 

i&c5 

96.0 
ICO.  5 

76.5 

81.0 

R:1 

147.0 

122.5 

3C5.0 

93.5 

S5.0 

78.8 

151.5 

127-0 

109-5 

97.5 

89.0 

82.6 

1=6.0 

131.  C 

113.5 

101.5 

93-0 

§6-3 

160.5 

135-0 

117-0 

105.0 

56.5 

89.8 

164.0 

139-0 

121.0 

103.5 

99-5 

93.1 

I68.0 

142-5 

124.5 

132.0 

105-0 

96.3 

172.0 

146.0 

128.0 

S.1 

105.5 

99-* 

175-5 

1*9-5 

131.5 

109.5 

102.4 

170.5 

153-0 

134.5 

121.5 

112.5 

105-3 

I08.I 

182.0 

156-0 

137.5 

124-5 

H5.5 

"185.0 

159-0 

1*1.0 

127-5 

118.5 

no. e 

138.0 

162.0 

1*3-5 

130-5 

221.0 

113-* 

191-0 

164-5 

1*6.5 

133-5 

123-5 

116.0 

193-5 

167-5 

140.5 

136-5 

126.0 

118.4 

196.5 

170.5 

152.0 

15S:? 

128.5 

120.  S 

199-0 

172,5 

15*. 5 

131-0 

123-2 

202.0 

175-0 

157-0 

144.0 

133-5 

125- C 

204.5 

207.0 

177-5 

180.0 

159-5 

162.0 

146.  * 
143-5 

135-5 

138.0 

127-1 

129-2 

209-5 

132.5 

13*. 0 

151.5 

140.0 

'  131-2 

212.0 

185-0 

166. 5  ' 

153.0 

142.0 

132-5 

214.5 

187.0 

16S-0 

155-0 

144.0 

134-9 

217.0 

139-5 

171-0 

157-C 

146  0 

137-1 

219.0 

191-5 

173-0 

159-5 

14S.0 

139-2 

221-5 

193-5 

175-5 

161-5 

150.0 

141.4 

223-5 

I96.O 

177-5 

163-5 

152-0 

1*3-5 

Z26.0 

193.0 

179-5 

ISI.5 

165.5 

154.0 

1*5-6 

228.0 

200.0 

167.0 

155- 5 

1*7-5 

230.5 

202-0 

183-5 

159-0 

157.5 

149-4 

232-5 

204.0 

185-5 

171-0 

159-5 

151-3 

234.5 

206.5 

187-5 

173.0 

161-5 

15M 

236-5 

208.5 

I60.5 

175-0 

I63.O 

15*.  5 

338. 5 

210.0 

191-5 

177.0 

178.5 

130. 5 

16*.5 

156.7 

2*0.5 

212.0 

193.0 

166.5 

1=3.4 

242.5 

214.0 

195-0 

168.5 

160.0 

2*4-5 

216.0 

197-0 

182.0 

170. 0 

161.6 

216 
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|  BUOYAKCY  0?  USSSOA  VAPOR  -  AIR  HIaTP§£  WITS  j 

|  AEROSOL  FRACTION  AS  A  PARAMETER  j 

i 

] 

Introduction  j 

I 

Whether  a  gaa  cloud  risen  in  air  (and  if  so,  its  rate  of  rise)  j 

is  determined  by  the  relative  densities  of  the  gas  cloud  and  : 

the  surrounding  air.  The  information  about  the  buoyancy  of  the  cloud  J 

| 

I  is  important  in  determining  its  dispersion  iu  the  atmosphere.  ■ 


I 

1 
j 

2 
J 

f 

I 

I 

i 

I 

i 

I 
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Pure  ammonia  vapor  is  lighter  than  ambient  air,  even  at  the  low 
saturation  temperature  (-28  °P) .  However ,  when  the  assaonia  vapor 
is  liberated  from  a  boiling  pool  of  liquid  araonia  on  water,  or  ou  Land, 
liquid  droplets  of  LHH-  are  thrown  into  the  vapor  phase  in  the  form  of 
aerosols.  The  density  of  such  an  aerosol  laden  ammonia  vapor  i3  oc-re 
than  that  of  pure  vapor.  Also  because  of  the  affinity  of  anaemia 
for  water  and  its  reaction  with  water,  several  phenomena  take  place 
when  moist  air  mixes  with  aerosol  laden  srsonia  vapor.  Ths  air  is 
cooled,  because  of  the  vapor's  low  temperature,  resulting  in  the 
possibility  of  water  vapor  condensation.  Ammonia  say  go  into  solution 
with  water,  liberating  heat.  A  part  of  this  heat  goes  to  heat  the 
mixture  of  the  gase3  and  another  part  to  evaporate  the  liquid  droplets 
of  ammonia.  The  final  equilibrium  condition  depends  on  the  relative 
amount e  of  aerosol,  humidity,  the  air  temperature,  and  the  mass  ratio 
of  vapor  to  air. 

To  obtain  the  final  fixture  condition  for  KH^ t  vapor  mixing  with 
moist  air  for  any  given  initial  conditions,  is  an  extremely  difficult 
and  involved  process.  However,  to  obtain  some  basic  understanding 
of  the  buoyant  nature  of  the  gas  with  different  aerosol  fractions,  s 
simplified  problem  is  worked  out£  wherein  only  dry  air  is  considered. 
Finally,  equations  are  developed  tc  describe  the  mixing  process  with 
moist  air  and  an  approach  to  the  solution  Is  discussed. 
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1)  Mixing  of  Saturated  Ammonia  Vapor  with  Bone  Dry  Air 

To  obtain  the  gaseous  xaixture  density,  the  mass  balance  equation 
and  energy  balance  equations  are  written  and  si Ived  for  the  final 
equilibrium  mixture  temperature.  Knowing  this  temperature,  the  mixture 
density  is  obtained.  In  deriving  the  equations,  the  following 
assumptions  are  made: 

a  All  equations  are  derived  on  the  basis  of  one  unit  mass  of 
bone  dry  air; 

•  The  gases  are  considered  to  be  perfect  gases; 

e  The  final  mixture  is  in  a  pure  gaseous  state  (no  droplets); 

e  Total  pressure  is  one  atmosphere. 

Let  F  **  Mass  of  cold  ammonia  vapor  (including  the  aerosol)  mixing 
with  unit  mass  of  air; 

F  11  Fraction  of  the  mass  of  ammonia  vapor  which  is  in  the  form 

u 

of  aerosol; 

I  «=  Enthalpy  per  unit  mass  of  any  substance  (subscripted) ; 

T^  3  Air  temperature 

T^  **  Initial  vapor  temperature  (saturation  temperature) 

^mix  =  temperature  °f  the  a*r  vapor  mixture. 

a)  Mass  Conservation 

Mass  of  air  +  mass  of  vapor  =  mass  of  final  mixture 

b)  Energy  Conservation 

Since  the  total  pressure  is  a  constant,  the  energy  conservation 
equation  becomes  the  enthalpy  conservation  equation,  i.e.. 


VV  +  F  1fl  'nf  <Ti>  +  CEi)J  -  <W 


+  F  I 


T(sat) 


vap 


(T  ,  ) 
mix 


(31) 
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The  air-vapor  mixture  temperature  is  obtained  from  the  above 
equation. 


c)  Equation  of  State 

R  R 

pup  m  p  £  2' 

F  a  p  ‘ i  mix  p  .  mix 
a  mix 


where 

’'‘mix 


(1  +  F)  p  p 
“  _ _ *  a  v 

(Uv  +  F  U  J 
v  a 


Ry  ~  universal  gas  constant 

p  *>  atmospheric  pressure 


p  "  molecular  weight 

X£  the  gases  can  be  assumed  to  have  constant  specific  heats, 
then  equation  B1  can  be  rearranged  to  give 


mix 


[c  T  +  F  (c  T.  -  F.  X)] 
— gj£— * _ ?,v i L 

IF  c  +  c  ] 

P.V  p,aJ 


(B2) 


(23) 


(24) 


where  c’s  are  the  constant  pressure  specific  heats  and  X  is  the  average 
latent  heat  of  vaporization  of  liquid  ammonia. 

Also  from  equations  B2  and  B3  we  get 

Fractional  Density  Deviation  »  A  «•  f  ~ -  1  1  (B5) 

V  *  >  (X  +  —  F)  J 

P 

V 

It  is  noted  that  physically  the  mixture  temperature  cannot  oe 
lass  than  the  saturation  temperature  of  ammonia  at  atmospheric  pressure. 

This  condition  imposes  an  upperbound  on  the  acceptable  value  of  the  mass 
ratio  of  vapor  to  air  (F)  for  a  given  fractional  aerosol  mass  in  vapor. 

By  substituting  Tmix  “  T,  in  equation  B4  and  rearranging  we  get 
for  the  critical  F, 

1 
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m 


) 

* 

! 


F 

c 


c 

Jh£ 


V 


and  F  <  F 
—  c 


(B6) 


Substituting  B6  in  B5  and  using  the  numerical  values  for  the 
molecular  weights  of  air  and  ammonia 


(1  +  Fc) 

(1  +  1.6985  F  ) 
c 


il 


(B7) 


Figure  B-l  shows  the  fractional  density  deviation  as  a  function 
of  vapor  air  ratio  (F)  with  the  aerosol  fraction  (F^)  as  a  parameter. 

The  critical  F  value  and  the  Critical  A  values  are  indicated  for  each 
FT .  The  buoyant  and  nonbuoyant  regions  (i.e.,  A  negative  and  A  positive, 

Lf 

respectively)  are  indicated  as  also  the  forbidden  region  consistent 
with  the  mixture  temperature  condition. 


Figure  B-2  indicates  the  mixture  temperature  for  different  F  and 
F_  values. 

•U 

2)  Mixing  of  Saturated  Ammonia  Vapor  with  Hoist  Air 
In  formulating  this  problem  we  assume  that: 

•  The  mixing  of  air  and  vapor  is  adiabatic. 


®  At  equilibrium  condition  after  mixing,  the  mixture  contains 
air,  ammonium  hydroxide  solution,  water  vapor  and  ammonia 
vapor.  The  last  two  components  are  in  equilibrium  with  the 
ammonium  hydroxide  solution. 


9 

« 

9 

e 


The  ammonium  hydroxide  solution  is  in  the  form  of  aerosols  and 
does  not  precipitate  out. 

The  volume  occupied  by  NR/ OH  drops  is  very  small. 

*4 

The  total  pressure  is  constant. 

A  unit  mass  of  dry  air  basis  is  used. 
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m 


2  3  4  5 

p  =  Mass  of  NH3  Vapor  at  -23°F 
Mass  of  Air 


FRACTIONAL  DENSITY  DEVIATION  OF  AJViMONIA  VAPGR- 
AiR  MIXTURE  WITH  AEROSOL  FRACTION  AS  A  PARAMETER 
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Air  Temperature  «  70°F 
Relative  Humidity  =  0% 


J _ _ _ - _ 1 _ _ _ _ L _ _ ...I . . i  , 

1.0  2.0  3.0  4.0  5.0 

F  =  Mass  of  NH3  Vapor  at  -28°F 
Mas  of  Air 


FIGURE  B— 2  TSMFERATURt  OF  AMMONIA  VAPOR-AIR  MIXTURE 
WITH  AEROSOL  FRACTION  AS  A  PARAMETER 


Let 


*  mass  of  gases  other  than  air  in  the  vapor  phase  (only  water 
and  ammonia) 

*  ctass  of  NH^OH  liquid  droplets  formed 
x  *  mass  fraction  of  water  in  the  liquid 

y  *  mass  fraction  of  water  vapor  in  the  water  vapor  ammonia  vapor 
mixture  (air  excluded) 

w^  »  initial  mass  of  water  vapor  in  air 

«  partial  pressure  of  air  in  the  mixture  of  gases 

V  •*  volume  of  the  mixture 


a)  Ammonia  Hags  Balance* 

(1  -  x)^  -J-  (1  -  y)Mv  -  F 

b)  Water  Mass  Conservation 
x  +  y  Mv  = 

c)  Phase  Equilibrium  Relation 

y  ■  f(x-  W  P  "  pa> 

d)  Partial  Pressure  of  Air 


p  V  -  —  T  , 
ra  u  mix 
£ 

e)  Volume  Relation 


where 


P..  = 


pv(s’  Tmix’ 


P  -  Pa) 
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f)  Relation  Between  Hass  of  Vepor  end  Mass  of  Liquid 


—  •  f(x,  T  ,  ,  p  -  p  ) 
'  *  mix*  *  *V 


Initial  enthalpy  of  air  +  Initial  enthalpy  of  saturated  ammonia 
vapor  ■  Final  enthalpy  of  mixture 

The  above  form  a  set  of  seven  coupled,  equations  for  the  seven 
unknowns:  M^,  pa,  x,  y  and  V.  The  solution  of  these 

equations  for  given  initial  values  of  F,  w.  and  T  is  extremely 
difficult.  However,  in  principle  the  solution  can  be  obtained  provided 
the  functional  relationships  connecting  the  various  parameters  for 
an  ammonia-water  system  are  known  completely. 
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HAZARD  ASSESSMENT  PROCEDURES  AHD 
APPLICATION  TO  SEVERAL  SPILL  SIZES 


Introduction 

The  procedures  for  calculating  the  hazards  caused  by  a  spill  of 
liquid  anhydrous  aaaaonla  (LNH^)  are  presented  in  this  appendix. 

The  procedures  ara  Illustrated  with  epecific  numerical  examples  for 
the  spill  of  3*000  tons  of  LNH^  on  water.  The  examples  Include  the 
calculation  of  vapor  concentration  at  ground  level  in  section  1,  asd 
in  section  2  the  calculation  of  ammonium  hydroxide  concentration  in 
water.  At  the  end  of  each  section,  plots  are  given  which  facilitate 
quick  hazard  aseessaent  for  several  spill  quantities  and  environmental 
conditions . 

Section  1:  Vapor  Dispersion 

The  procedure  for  calculating  the  ground  level  vapor  concentration 
can  be  summarized  by  the  following  steps: 

1.  Calculate  the  maximum  radius  of  spread  (Figure  5-7) 


2.  Calculate  the  mass  and  volume  of  saturated  vapor  liberated 

using  the  value  of  the  partition  function*  the  mass  of 

liquid  spilled,  and  the  density  of  vapor. 

3.  Calculate  the  fractional  density  deviation  of  the  vapor. 

4.  Calculate  the  buoyancy  flux  (section  6. 7. 2,1)  and  then  the 

buoyancy  length.  Knowing  this  and  the  downwind  distance* 
tha  vapor  plume  rise  can  be  obtained. 

5.  Knowing  the  value  of  the  potential  temperature  gradient 

in  the  atmosphere,  calculate  the  stratification  parameters 
(se2  section  6. 7. 2.1).  From  this  the  ceiling  height  s 


to  -Itch  the  cloud  rises  can  be  obtained.  Note  that 
cloud  reaches  a  ceiling  height  only  in  stable  atmospheres 
(atmospheres  S  and  F) .  Also,  in  the  case  of  fumigation 
type  of  atmospheres,  the  atmosphere  should  be  considered 
to  be  unstable  up  to  the  ceiling  height.  This  latter 
ceiling  height  should  be  obtained  from  weather  dsta 
(see  section  6. 5. 2. 2  ii) . 

6.  For  the  height  of  rise  of  cloud,  use  the  smaller  of  the 
values  obtained  in  steps  4  and  5. 

7.  Obtain  the  effective  downwind  distance  by  adding  to 
the  actual  distance  10  times  the  maximum  radius  of  spread. 

8.  Obtain  the  dispersion  parameter  values  (c's)  for  the 
appropriate  atmospheric  type,  at  the  effective  distance, 
from  Figures  6-i7a  and  6-17b. 

9.  Using  the  appropriate  formula  depending,  on  the  type  of 
atmosphere,  obtain  the  ground  level  concentrations. 

The  procedure  is  illustrated  below  with  a  specific  numerical 
example. 


Data 


Density  of  liquid  ammonia  at  -28°F 

42.56 

lbs/ft3 

Density  of  saturated  vapor  at  -28°F 

0.05556 

Iba/tt* 

Ambient  temperature 

68.0 

°F 

Density  cf  air  at  68°F 

0.07488 

lbs /ft3 

Mass  of  spilled  instantaneously 

«£ 

3000 

A 

tons 

■=>  6  x  10 

lbs 

Volume  of  liquid 

V£ 

1.056  x  106 

gallons 

Partition  ratio  (assumed) 

P 

0.6 

Liquid  regression  rate  (see  page  126) 

• 

y 

3.89  x  1C"3 

ft  io 
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Cal cul a tlgoB 


Maximum  radius  of  spread 
(Figure  5-7) 


2.5  (Sai) 


G.375 


453.75 


Mass  of  vapor  liberated 
Volume  of  vapor  liberated 


V 

v 


Fractional  density  deviation  A.  “ 

of  LNH^  1 

Fractional  density  deviation  A„  * 

of  vapor 


<i-p)ma 


(1-  “ ) 
pw 

(l-~) 

pair 


2.4  x  10  lbs 

4.32  x  107  ft3 


0.318 


«  0.259 


Maximum  vapor  liberation  rate  m 


Time  for  complete  evaporation  t 
of  LNHj  (equation  5-5)  t 


^max  p£  7(1"p)  = 


"  Vi 

. 8  y2  Ai. 


1.713  x  10  lbs/s 


123.4 


a)  Dispersion  of  Vaoor  Under  Unstable  Atmosphere  (c) 


Wind  velocity  (assumed)  U  *  20  knots  *  33.8  ft/s 

Downwind  distance  at  which  ground  x  «  5  km  «  2.7  nautical  miles 
level  concentration  is  needed  , 


Buoyancy  flux  Fp  (see  section  6. 7. 2.1)  “ 


4.32  s  10  x  32.2  x  0.259 
v  x  123.4 

9.29  x  105  ft4/s 


Buoyancy  length  Z  «*  — jj- 
r  IT 


Hence 

Effective  source  distance  “ 

x  «  x  +  10  R 
e  max 

Dispersion  parameters  for  x  ,  in 
atmosphere  c  (see  Figures  6-373 
and  6-lTo) 


**  24.0 


Height  of  rise  of  center  line  of  z  *  S.Sx^ 
vapor  cloud  (equation  6-56  with 
x  in  feet) 


s  *  3.3  x  (16400) J  2130  ft 


«  20936 


z  -  1050 
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Ground  level  concentration 
(set  equation  6-46b) 

2  x  2.4  s  106  e~1/,2U0' 

i0)2 

jo 1 

(2tr)1,5  18372  x  1050 

Concentration  in  ppa 

-  1.1  x  10-5 

-  1.1  x  10"5  x  (22.825  x 

10°) 

lb/ft3 

=*  251 

ppm 

b)  Dispersion  in  Stable  Atmosphere  (?) 


Potential  temperature  gradient 
in  the  atmosphere  (assumed) 

Wind  velocity 

Atmospheric  temperature 

Buoyancy  length  (1  ) 

P 

Height  of  rise  of  center  line  of 
vapor  cloud  in  a  distance  of  x 
feet  (see  equation  6-56) 

Brunt  Vaisalla  frequency 
(see  page  193) 

Stratification  parameter 
(see  page  193) 

Hence  the  maximum  height 
reached  by  the  vapor  cloud 
(equation  6-61)  given  by  the 
plume  theory* 

Downwind  dietance  at  which 
this  maximum  height  is  raachad 

Downwind  distance  at  which 
ground  level  concentration  is 
needed  (assumed) 

Effective  distance 

Dispersion  parameters  for 
atmosphere  F  at  the  distance 
x  (from  Figures  6-17a  and 
6^1?b) 


5  °K/km 


l5 

mph 

U 

17.3 

ft/s 

288 

°K 

!e 

2389 

ft 

U3 

z 

15.3  x2/3 

ft 

u 

c 

1.304  x  10“2 

7.3 

rad/s 

-  0.234 

5 

2383  x  1.304  x  10~2 

Zc 

2383  x  (.234) 2/3 

/aqq  \**5 

-  989  ft 

X* 

(sy 

-  520  ft 

nautical 

y 

110 

miles 

J* 

\ 60,800 

ft 

X 

x  +  10  R 

-  65,337  ft 

e 

max 

0 

1640 

ft 

7 

a 

197 

ft 

As  already  discussed  in  section  5.5  (page  128),  the  behavior  of  a  vapor 
cloud  liberated  by  a  massive  spill  is  likely  to  be  more  accurately 
described  by  the  puff  theory.  However,  because  of  the  lack  of  data 
for  massive  spills  of  LNH^  and  in  order  to  conform  to  the  findings  of 
the  present  test  program,  we  use  the  more  conservative  value  given  by 
the  plume  theory;  that  is,  zw  is  assumed  to  be  989  ft  (instead  of  sm 
■  1710  ft  given  by  puff  theory) . 
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sssssssasseaass. 


2  s  2.4  x  106  e~1/2 


tration  (equation  4-46b) 

(2ir)3/2  x  1S402  x  197 

. 

-20*0625 

e 

lbs /ft3 

C2 

0.044 

ppm 

Dispersion  Under  Stable  (F) :  Fumigation  Condition 

Assuming  a  ceiling  height 
(this  has  to  be  obtained 
from  weather  data) 

h  1,000 

ft 

Downwind  distance  -  same 
as  in  case  b) 

x  10 

nautical 

miles 

Dispersion  coefficient  at 
effective  distance  (x  )  for 
fumigation  condition 
(atmosphere  C)  from  Figures 

6-17a  and-6-17b) 

c  4,900 

y 

ft 

Hence  ground  level  concen¬ 
tration  (see  equation  6-47)  = 

2.4  x  106 

2rr  x  49002  x  1000 

a 

1.59  x  1G"5 

lbs/ft3 

a 

361.0 

ppm 

The  ground  level  concentrations  for  various  atmospheric  conditions 
are  plotted  as  functions  of  downwind  distances  in  Figures  C-l,  C-2,  C-3, 
and  C-4  respectively  for  100,  500,  1000,  and  3000  ton  spills.  The 
curves  for  atmosphere  F  indicate  the  fumigation  condition,  which  is 
the  worst  condition  as  can  be  seen  from  the  examples  in  b)  and  c) . 

In  calculating  these  stable  atmosphere  curves,  the  ceiling  height 
for  each  value  of  spill  was  evaluated  from  the  plume  theory  assuming 
a  temperature  gradient  of  5°K/km.  It  is  recalled  here  that  under 
fumigation  conditions,  the  ceiling  height  depends  cn  the  weather  and 
as  such  has  to  be  obtained  from  weather  data.  Therefore,  the  curves 
given  for  atmosphere  F  are  to  be  used  with  care  and  where  necessary 
have  to  be  recalculated  using  the  proper  ceiling  height  data. 
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Spill  Quantity:  100  Tons 


Downwind  Distance  (Nautical  Miles) 


Downwind  Distance  (Nautical  Milos) 


I 


Downwind  Distance  (Nautical  Miles) 


Section  2:  Water  Dispersion 


Introduction 

The  ammonium  hydroxide  (HE,  OH)  rosined  by  the  dissolution  of  LNH. 

•+  J 

in  water  is  dispersed  by  molecular  diffusion  in  quiescent  water  and  by 
turbulent  mixing  in  flowing  streams.  In  the  case  of  flowing  rivers 
and  streams,  the  downstream  NH^OH  concentration  is  calculated  by  using 
the  water  dispersion  models. 

In  this  section,  numerical  examples  are  given,  indicating  the 
calculating  procedure  to  obtain  the  NH^OH  concentration.  A  3,000-ton 
surface  spill  of  LNHL  is  considered. 

In  the  examples  shown  below,  boch  the  near  field  and  far  field 
calculations  are  indicated.  The  near  field  approximation  give6  NH^OH 
concentrations  in  the  immediate  vicinity  of  the  spill.  The  far  field 
model  is  used  for  distances  beyond  which  there  is  no  appreciable 
vertical  mixing.  Eecause  of  the  nature  of  the  assumption  made  in 
the  derivation  of  these  two  models,  there  is  bound  to  be  some 
discontinuity  in  the  values  of  NH^OH  concentrations  predicted  by  the 
models  In  the  common  regions  of  applicability.  Hence,  proper  judgment 
has  to  be  exercised  in  interpreting  the  concentration  values  obtained 
in  these  regions. 

Calculations 

In  the  calculations  shown  below,  the  following  assumptions  are 

made: 

a  Even  though  ammonium  hydroxide  is  lighter  than  water,  we  assume 
that  this  does  not  have  any  significant  effect  on  the  dilution 
process . 

»  No  evaporation  of  the  ammonium  hydroxide  takes  place  during 
dilution. 

e  The  spill  xs  instantaneous. 
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Spill  end  environmental  conditions: 
Quantity  of  Spill 

Density  of  LNH^ 

Partition  Ratio 
Location  of  Spill 


River  width 
River  depth 

Roughness  factor  for  riverbed 

Mean  velocity  in  the  river 

Downstream  distance  at  which 
maximum  concentration  is  to  be 
known  for  near  field  approxima¬ 
tion 

Downstream  distance  at  which  maximum 
concentration  is  to  be  known  for 
far  field  approximation 

Total  quantity  of  100%  ammonium 
hydrcxide  generated  in  a 
surface  spill 

Hydraulic  depth  of  stream  R^ 

Shear  velocity  (see  Table  6-3)  u* 


ra  3,000  tons  »  6  x  10  lbs 
a  IQ**  gallons 
“  42.56  lbs/ft^ 

0.6 

“  In  the  middle  of  a 
nontidal  river 

-  1,000  ft 

“  50  ft 


«  U  « 3  ft/sec 

»  0.5  nautical  miles 
r  3,000  ft 


«  10  nautical  miles 
s  60,000  ft 

°  0.6  x  6  x  10^  “  3.6  x  10^  lbs 


1000  X  50  ,  r  ,  r-  r*. 

“  1000  +  2*50  ■  45-45  ft  . 
*»  3.3  x  .03  x  3  x  (45.15) 
-  0.16  ft/s 


The  concentration  at  any  point  downstream  of  a  surface  spill  is 

a  maximum  on  the  water  surface.  For  the  near  field  analysis,  the 

source  is  considered  to  be  an  area  source  because  of  the  rapid  spread 

of  LNH^  on  water.  For  underwater  releases,  at  sufficient  depths  all 

of  the  LNBg  should  be  assumed  to  dissolve  in  water.  The  location  of 

the  maximum  NH.OH  concentration  downstream  for  underwater  release  is 
4 

not  necessarily  at  the  level  of  the  release.  This  is  because  of  the 
buoyant  nature  of  NH^OH.  However,  for  large  downstream  distances,  the 
location  of  release  hardly  influences  the  concentration. 


Near  Field  Approximation 

Time  at  which  concentration  is  maximum 
at  x  *  3000  ft 

Estimated  turbulent  diffusion  coefficient  e 
(see  Table  6-3)  for  narrow  river 


x 

U 


3000 

3 


100  sec 


Maximum  radius  of  spill ,  R 


max 


Equivalent  thickness  of  100%  concentration 
NH4OH  layer  (see  Figure  6-19  and  equation 
6-49) 

Maximum  NH^OH  concentration  (see 
equation  6-52) 


«  0.067  x  0.16  x  45.45 
*=  0.49  ft^/sec 

«  2.5  (gallons  spilled) 
S  454  ft 

3.6  x  106 

ct  -  — — *  ■■  .■■■—  —  —  —  —  — -  - 

42.56  xitx  (454) 2 

-  0.144  lb/ft3 
*=  2300  ppm  B  2300  mg /I 


0.375 


Far  Field  Approximation 

It  can  be  seen  that  the  near  field  approximation  gives  a  constant 
concentration  for  sufficiently  large  time.  This  occurs  because  the 
derivation  of  the  equation  neglects  longitudinal  dispersion  -  the  only 
mechanism  by  which  dilution  will  take  place  once  the  vertical  distri¬ 
bution  of  NH^OH  in  the  water  becomes  uniform. 


6- 


Therefcre,  for  concentration  predictions  at  long  times ,  equation 
is  used  to  allow  for  longitudinal  dispersion. 


Longitudinal  dispersion  coefficient 
(see  equation  6-55) 


Time  at  which  concentration  ie 
maximum  at  x  •=■  61,000  ft 

Maximum  NH,0H  concentration 
(see  equation  6-53) 


o  23.3  V 

=  23.3  x  .16  x  45.45 
«=  169.4  ft^/sec 

x  60,000  „„  ../I,. 

»  —  =  — 20,000  sec 

u  ^ 


3.6  x  10 


1000  x  50/4ir  x  169.4  x  20,000 
«  .011  lb/ft3 
“  177  ppm  =  177  mg/1 
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Discussion 


To  facilitate  a  rapid  hazard  assessment  of  the  water  pollution 
from  an  LNH^  spill  on  water t  Figures  C-5  through  C-8  are  drawn.  These 
Figures  indicate  the  maximum  NH^OH  concentration  as  a  function  of  the 
downstream  distance  in  nontidal  rivers.  The  stre"»  width  is  treated 
as  a  variable  in  addition  to  the  spill  quantity.  The  stream  depth  is 
kept  a  constant  at  50  feet. 

The  plots  were  generated  using  the  two  models  illustrated  earlier. 
However,  in  the  overlapping  regions,  continuity  between  the  two 
models  is  provided  using  subjective  judgment  (to  give  conservative 
concentration  estimates) . 
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Water  Surface  NH4OH  Concentration  (PPM) 


.1  1  10  100 

Downstream  Distance  (Nautical  Miles) 


FIGURE  C-7 

Water  Dispersion  Hazard  Assessment  Plot  for 
Surface  Spi.ll  Quantity  of  1000  Tons 


